bzoj1096: [ZJOI2007]仓库建设 斜率优化DP

纠结的公式推导。。。g[i]表示到第1个到第i个仓库的货物总和，sum [i]表示把i以前的所有货物都送到i要多少费用，当然还要加上i点建工厂的费用。

#include <iostream>
#include <cstring>
#include <algorithm>
#include <cstdio>
#include <cmath>
using namespace std;
#define ll long long
#define maxn 1100000
struct node{ll x,p,c;}a[maxn];
ll f[maxn],q[maxn],h,t,n,sum[maxn],g[maxn];
bool cmp(node aa,node bb){return aa.x<bb.x;}
double getk(ll k,ll j)
{
    double tmp=(f[k]-f[j]-sum[k]+sum[j]+g[k]*a[k].x-g[j]*a[j].x)/(g[k]-g[j]);
    return tmp;
}
int main()
{
    scanf("%lld",&n);
    for(int i=1;i<=n;i++)
    {
        scanf("%lld%lld%lld",&a[i].x,&a[i].p,&a[i].c);
    }
    sort(a+1,a+1+n,cmp);
    for(int i=1;i<=n;i++)
    {
        sum[i]=sum[i-1]+g[i-1]*(a[i].x-a[i-1].x);
        g[i]=g[i-1]+a[i].p;
    }
    h=t=1;
    for(int i=1;i<=n;i++)
    {
        while(h<t&&getk(q[h+1],q[h])<double(a[i].x)) h++;
        f[i]=f[q[h]]-sum[q[h]]+sum[i]-g[q[h]]*(a[i].x-a[q[h]].x)+a[i].c;
        while(h<t&&getk(q[t],q[t-1])>getk(i,q[t])) t--;
        q[++t]=i;
    }
    printf("%lld\n",f[n]);
    return 0;
}