bzoj 3963: [WF2011]MachineWorks (DP+CDQ分治)

题目描述

传送门

题解

将所有的机器按照d从小到大排序。
f[i] 表示买入第i台机器后的最大收益,则 f[i]=max{f[j]+(d[i]d[j]1)g[j]+r[j]}p[i] ,其中 j<i
我们将式子写开,得到 f[i]=f[j]+d[i]g[j]d[j]g[j]g[j]+r[j]p[i]
移项,得到 f[j]d[j]g[j]g[j]+r[j]=f[i]p[i]d[i]g[j]
那么 f[i]p[i] 可以看成是经过点 (g[j],f[j]d[j]g[j]g[j]+r[j]) 的斜率为 d[i] 的直线的最大截距.
我们对下标进行CDQ分治,将[l,mid]得到的点维护成一个斜率单减的凸壳,[mid+1,r]中的直线按照斜率递减的顺序依次进行匹配,计算答案即可。

代码

#include<iostream>
#include<cstdio>
#include<cstring>
#include<algorithm>
#include<cmath>
#define N 100033
#define eps 1e-9
#define LL long long
using namespace std;
const double inf=1e19;
int n,C,D,top,st[N];
LL f[N];
struct data{
    LL d,r,p,g,id;
    double k;
}a[N],np[N];
struct point{
    LL x,y; bool pd;
    bool operator <(const point &b) const {
      return x<b.x||(x==b.x&&y>b.y);
    } 
}p[N],np1[N];
int cmp(data a,data b)
{
    return a.d<b.d;
}
double getk(int i,int j)
{
    if (p[i].x-p[j].x==0) return -inf;
    return (double)(p[i].y-p[j].y)*1.0/(double)(p[i].x-p[j].x)*1.0;
}
void divide(int l,int r)
{
    if (l==r) {
        f[l]=max(C-a[l].p,f[l]); p[l].pd=0;
        if (f[l]<0) f[l]=-1,p[l].pd=1;
        p[l].x=a[l].g;
        p[l].y=f[l]-a[l].d*a[l].g-a[l].g+a[l].r;
        return;
    }
    int mid=(l+r)/2;
    int l1=l; int l2=mid+1;
    for (int i=l;i<=r;i++)
     if (a[i].id<=mid) np[l1++]=a[i];
     else np[l2++]=a[i];
    for (int i=l;i<=r;i++) a[i]=np[i];
    divide(l,mid);
    int top=0;
    for (int i=l;i<=mid;i++) {
        if (p[i].pd) continue;
        while (top>=2&&getk(i,st[top])+eps>getk(st[top],st[top-1])) top--;
        st[++top]=i;
    }
    if (top) {
        int j=1;
        for (int i=mid+1;i<=r;i++) {
            while (j<top&&getk(st[j],st[j+1])>-a[i].d) 
             j++;
            int k=st[j];
            f[a[i].id]=max(f[a[i].id],p[k].y+p[k].x*a[i].d-a[i].p);
        }
    }
    divide(mid+1,r);
    l1=l; l2=mid+1;
    for (int i=l;i<=r;i++)
     if ((p[l1]<p[l2]||l2>r)&&l1<=mid) np1[i]=p[l1++];
     else np1[i]=p[l2++];
    for (int i=l;i<=r;i++)
      p[i]=np1[i];
}
int main()
{
    freopen("a.in","r",stdin);
    freopen("my.out","w",stdout);
    int T=0;
    while (true) {
        scanf("%d%d%d",&n,&C,&D); T++;
        if (n==0&&C==0&&D==0) break;
        for (int i=1;i<=n;i++) scanf("%I64d%I64d%I64d%I64d",&a[i].d,&a[i].p,&a[i].r,&a[i].g);
        ++n; a[n].d=D+1; a[n].p=0; a[n].r=0; a[n].g=0;
        for (int i=1;i<=n;i++) f[i]=-1;
        sort(a+1,a+n+1,cmp);
        for (int i=1;i<=n;i++) a[i].id=i;
        divide(1,n);
        printf("Case %d: %lld\n",T,f[n]);
    }
}
评论 1
添加红包

请填写红包祝福语或标题

红包个数最小为10个

红包金额最低5元

当前余额3.43前往充值 >
需支付:10.00
成就一亿技术人!
领取后你会自动成为博主和红包主的粉丝 规则
hope_wisdom
发出的红包
实付
使用余额支付
点击重新获取
扫码支付
钱包余额 0

抵扣说明:

1.余额是钱包充值的虚拟货币,按照1:1的比例进行支付金额的抵扣。
2.余额无法直接购买下载,可以购买VIP、付费专栏及课程。

余额充值