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;
}