BZOJ.1096.[ZJOI2007]仓库建设(DP 斜率优化)

题目链接

\(Solution\)

　　注: 刚学 可能过程有疏漏，仅供参考（其实我不是很懂这里是否需要证决策单调）。

　　抽象下问题，即将一个序列分成任意多段，设\(f[i]\)表示以\(i\)作为一个右端点时\([1,i]\)的最小值，则\(f[r]=f[l]+(l,r]这一段的价值\)，即\[f[r]=f[l]+cost[r]+\sum_{i=l+1}^{r-1}P[i]*(X[r]-X[i])\].

　　令\[S[x]=\sum_{i=1}^xX[i]\ ,\ W[x]=\sum_{i=1}^xP[i]*X[i]\]

　　化简式子得

\[f[l]+X[r]*(S[r-1]-S[l])-(W[r-1]-W[l])+cost[r]\]

　　假设对于i来说，之前有\(j\)比\(k\)更优，即

\[f[j]+X[i]*(S[i-1]*S[j])-(W[i-1]-W[j])+cost[i]>f[k]+X[i]*(S[i-1]*S[k])-(W[i-1]-W[k])+cost[i]\]

　　移项

　　\[\frac{f[j]+W[j]-(f[k]+W[k])}{S[j]-S[k]}<X[i]\]

　　所以如果满足上式，则\(j\)比\(k\)更优。

#include<cstdio>
#include<cctype>
#define gc() getchar()
typedef long long LL;
const int N=1e6+5;

int n,q
;
LL dis
,cost
,num,S
,W
,f
;

inline int read()
{
int now=0,f=1;register char c=gc();
for(;!isdigit(c);c=gc()) if(c=='-') f=-1;
for(;isdigit(c);now=now*10+c-'0',c=gc());
return now*f;
}
inline LL X(int j,int k) {return S[j]-S[k];}
inline LL Y(int j,int k) {return f[j]-f[k]+W[j]-W[k];}
inline LL Calc(int i,int p) {return f[p]+dis[i]*(S[i-1]-S[p])-W[i-1]+W[p]+cost[i];}

int main()
{
n=read();
for(int a,i=1;i<=n;++i)
{
dis[i]=read(), num=read(), cost[i]=read();
S[i]=S[i-1]+num, W[i]=W[i-1]+dis[i]*num;
}
int h=1,t=1; q[1]=0;
for(int i=1;i<=n;++i)
{
while(h<t && Y(q[h+1],q[h])<=dis[i]*X(q[h+1],q[h])) ++h;
f[i]=Calc(i,q[h]);
while(h<t && Y(i,q[t])*X(q[t],q[t-1])<Y(q[t],q[t-1])*X(i,q[t])) --t;
q[++t]=i;
}
printf("%lld",f
);

return 0;
}