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

来源:互联网 发布:理财软件哪款靠谱 编辑:程序博客网 时间:2024/05/17 08:34

基本和锯木厂一样。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 longint 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;}
阅读全文
0 0