2011 WorldFinal Machine Works

CDQ分治+斜率优化dp

首先考虑dp，有dp[i]表示到第i天所获得的价值最大是多少(机器卖出)，所以最后的答案就是dp[n+1]

然后推一推式子：dp[i]=max(dp[j]+(d[i]-d[j]-1)*g[j]+r[j]-p[j])

设当j>k时从dp[j]转移会比从dp[k]转移更优，那么有：

dp[j]+d[i]*g[j]-(d[j]+1)*g[j]+r[j]-p[j]>=dp[k]+d[i]*g[k]-(d[k]+1)*g[k]+r[k]-p[k]

设f[t]=dp[t]-(d[t]+1)*g[t]+r[t]-p[t]

所以有：f[j]+d[i]*g[j]>=f[k]+d[i]*g[k]

                      f[j]-f[k]>=d[i]*(g[k]-g[j])

设g[k]<=g[j]

所以有： (f[j]-f[k])/(g[j]-g[k])>=-d[i]

所以要使用斜率优化要满足d[i]单调递增((-d[i])单调递减)，g[i]单调递减

所以我们CDQ分治一下，将第一维d[i]排序，第二维g[i]分治一下就ok了

#include <iostream>
#include <cstdio>
#include <cstdlib>
#include <cstring>
#include <algorithm>
#include <cmath>
using namespace std;
const int maxn=100010;
int n,ca;
long long c,dd,dp[maxn];
struct da{long long d,p,r,g;}a[maxn];
struct po{long long ff,gg;}pp[maxn],now[maxn];
bool cmp(da aa,da bb){return aa.d<bb.d;}
bool cmp1(po aa,po bb){if (aa.gg==bb.gg) return aa.ff<bb.ff;return aa.gg<bb.gg;}
void cdq(int l,int r)
{
    int mid=(l+r)>>1; if (l==r) return;
    cdq(l,mid); int na=0,nb=0;
    for (int i=l;i<=mid;i++) 
     if (dp[i]>=a[i].p) 
      pp[++na]={(dp[i]-a[i].d*a[i].g-a[i].g+a[i].r-a[i].p),a[i].g};
    sort(pp+1,pp+na+1,cmp1);    
    for (int i=1;i<=na;i++)
    {
        while (nb>1 && (double)(now[nb].ff-now[nb-1].ff)*(now[nb-1].gg-pp[i].gg)>=(double)(pp[i].ff-now[nb-1].ff)*(now[nb-1].gg-now[nb].gg)) nb--;
        now[++nb]=pp[i];    
    }
    int cct=1;
    for (int i=mid+1;i<=r;i++)
    {
        while (cct<nb && (double)(now[cct+1].ff-now[cct].ff)>=(double)a[i].d*(now[cct].gg-now[cct+1].gg)) cct++;
        dp[i]=max(dp[i],now[cct].ff+now[cct].gg*a[i].d);
    }
    cdq(mid+1,r);
}
int main()
{
    while (scanf("%d%lld%lld",&n,&c,&dd))
    {
        if(!n && !c && !dd) break;
        for (int i=0;i<=n+1;i++) dp[i]=c;
        for (int i=1;i<=n;i++) 
         scanf("%lld%lld%lld%lld",&a[i].d,&a[i].p,&a[i].r,&a[i].g);
        a[++n].d=dd+1; a[n].p=a[n].r=a[n].g=0; sort(a+1,a+n+1,cmp);
        cdq(0,n); printf("Case %d: %lld\n",++ca,dp[n]);
    }
return 0;
}