bzoj1096 [ZJOI2007]仓库建设(斜率优化)

xiaoxiao2021-02-28  80

基本和锯木厂一样。dp[i]表示在i建一个锯木厂的最优方案。则 dp[i]=min{dp[j]+c[i]-c[j]-sw[j]*(sd[i]-sd[j]}+w[i]。 具体见锯木厂。

#include <cstdio> #include <cstring> #define ll long long int const N=1000000+10; int d[N],w[N],n,q[N],h=0,t=0; ll sw[N],sd[N],c[N],dp[N]; inline int read(){ int x=0,f=1;char ch=getchar(); while(ch<'0'||ch>'9'){if(ch=='-') f=-1;ch=getchar();} while(ch>='0'&&ch<='9') x=x*10+ch-'0',ch=getchar(); return x*f; } inline double slope(int k1,int k2){ return (dp[k1]-dp[k2]+c[k2]-c[k1]+sw[k1]*sd[k1]-sw[k2]*sd[k2])*1.0/(1.0*(sw[k1]-sw[k2])); } int main(){ // freopen("a.in","r",stdin); n=read();//注意读入的是sd。。。 for(int i=1;i<=n;++i){ sd[i]=read();int x=read();w[i]=read(); sw[i]=sw[i-1]+x; } for(int i=1;i<=n;++i) d[i]=sd[i]-sd[i-1]; for(int i=1;i<=n;++i) c[i]=c[i-1]+sw[i-1]*d[i]; for(int i=1;i<=n;++i){ while(h<t&&slope(q[h],q[h+1])<sd[i]) ++h; dp[i]=dp[q[h]]+c[i]-c[q[h]]-sw[q[h]]*(sd[i]-sd[q[h]])+w[i]; while(h<t&&slope(q[t],i)<slope(q[t-1],q[t])) --t; q[++t]=i; } printf("%lld",dp[n]); return 0; }
转载请注明原文地址: https://www.6miu.com/read-37236.html

最新回复(0)