bzoj3963[WF2011]MachineWorks cdq分治+斜率优化dp

Posted wsy

tags:

篇首语:本文由小常识网(cha138.com)小编为大家整理,主要介绍了bzoj3963[WF2011]MachineWorks cdq分治+斜率优化dp相关的知识,希望对你有一定的参考价值。

3963: [WF2011]MachineWorks

Time Limit: 30 Sec  Memory Limit: 256 MB
Submit: 270  Solved: 80
[Submit][Status][Discuss]

Description

你是任意性复杂机器公司(Arbitrarily Complex Machines, ACM)的经理,公司使用更加先进的机械设备生产先进的机器。原来的那一台生产机器已经坏了,所以你要去为公司买一台新的生产机器。你的任务是在转型期内尽可能得到更大的收益。在这段时间内,你要买卖机器,并且当机器被ACM公司拥有的时候,操控这些机器以获取利润。因为空间的限制,ACM公司在任何时候都只能最多拥有一台机器。
在转型期内,有若干台可能卖出的机器。作为先进机器的专家,对于每台机器Mi,你已经知道了其价格Pi和可以买入的日期Di。注意,如果不在第Di天买入机器Mi,那么别的人也会买走这一台机器,也就是说,以后你将没有机会购买这台机器了。如果ACM的钱低于一台机器的价格,那么你显然不可能买到这一台机器。
如果你在第Di天买入了机器Mi,那么ACM公司可以从第(Di)+1天开始使用这一台机器。每使用这台机器一天,就可以为公司创造出Gi美元的收益。
你可以决定要在买入之后的某一天,以一定的折扣价卖出这一台机器。收购市场对于每一台机器,都有一个折扣价Ri。你不能在卖出的那一天使用机器,但是你可以在卖出的那一天再买入一台新的。
在转型期结束后,ACM公司会卖掉当前所拥有的机器。你的任务就是最大化转型期间ACM公司可以得到的收入。
 

 

Input

输入包含若干组测试用例。每一组测试用例的第一行有3个正整数N,C和D。N是将会卖出的机器的台数(N<=10^5),C是在转型期开始时公司拥有的美元数量(C<=10^9),D是转型期持续的天数(D<=10^9)。
之后的N行每一行描述了一台机器的情况。每一行有4个正整数Di,Pi,Ri和Gi,分别表示这台机器卖出的时间,购买这台机器需要的美元数量,卖出这台机器的折扣价和使用这台机器可以得到的利润。这些数字满足1<=Di<=D,1<=Ri<Pi<=10^9且1<=Gi<=10^9.
最后一组测试用例后面的一行由3个0组成,表示输入数据。
 

 

Output

对于每一组测试用例,输出测试用例的编号,之后给出ACM公司在第D+1天结束后可以得到的最大数量的美元。
请依照下面给出的样例输出。
 

 

Sample Input

6 10 20
6 12 1 3
1 9 1 2
3 2 1 2
8 20 5 4
4 11 7 4
2 10 9 1
0 0 0

Sample Output

Case 1: 44

 

转移并不满足单调啊,所以需要用cdq分治来维护单调的转移
也可以写splay
这篇博客写得不错http://www.cnblogs.com/zj75211/p/8148800.html

 1 #include<bits/stdc++.h>
 2 #define N 100050
 3 #define ll long long
 4 #define inf 2147483647
 5 using namespace std;
 6 int n,m,D,s[N];ll f[N];
 7 struct date{int d,p,r,g,id;}q[N],a[N];
 8 //struct point{int x;ll y;}p[N],b[N];
 9 bool cmp1(date a,date b){return a.d<b.d;}
10 ll X(int i){
11     return q[i].r+f[q[i].id]-q[i].p-(ll)q[i].g*(q[i].d+1);
12 }
13 int Y(int i){
14     return q[i].g;
15 }
16 
17 double G(int i,int j){
18 //    if(Y(i)==Y(j))return (double)(X(i)>X(j)?inf:-inf);
19     return (double)(X(i)-X(j))/(double)(Y(i)-Y(j));
20 }
21 void solve(int l,int r){
22     if(l>r)return;
23     if(l==r){
24         f[l]=max(f[l],f[l-1]);
25         return;
26     }
27     int mid=(l+r)>>1;
28     int p1=l,p2=mid+1;
29     //for(int i=l;i<=r;i++)q[i]=a[i];
30     solve(l,mid);int tp=0,h=1;
31     for(int i=l;i<=mid;i++){
32         if(q[i].p>f[q[i].id])continue;
33         if(Y(i)==Y(s[tp])&&tp){
34             if(X(i)>X(s[tp]))tp--;
35             else continue;
36         }
37         while(tp>1&&G(s[tp],s[tp-1])<=G(i,s[tp]))tp--;
38         s[++tp]=i;
39     }
40     for(int i=mid+1;i<=r;i++){
41         while(h<tp&&G(s[h],s[h+1])>=-q[i].d)h++;
42         f[i]=max(f[i],f[i-1]);
43         if(h<=tp)f[i]=max(f[i],f[q[s[h]].id]-q[s[h]].p+q[s[h]].r+(ll)q[s[h]].g*(q[i].d-q[s[h]].d-1));
44     }
45     solve(mid+1,r);
46     p1=l;p2=mid+1;
47     for(int i=l;i<=r;i++){
48         if(p1<=mid&&q[p1].g<=q[p2].g)a[i]=q[p1++];
49         else if(p2<=r)a[i]=q[p2++];
50         else a[i]=q[p1++];
51     }
52     for(int i=l;i<=r;i++)q[i]=a[i];
53 }
54 int main(){
55     int cas=0;
56     while(scanf("%d%d%d",&n,&m,&D)!=EOF){
57         cas++;
58         if(!n&&!m&&!D)break;
59         memset(f,-1,sizeof(f));
60         f[1]=f[0]=m;int d,p,r,g;
61         for(int i=1;i<=n;i++)
62         scanf("%d%d%d%d",&d,&p,&r,&g),
63         q[i]=(date){d,p,r,g,i};
64         q[++n]=(date){D+1,0,0,0,n};
65         sort(q+1,q+1+n,cmp1);
66         for(int i=1;i<=n;i++)q[i].id=i;
67         solve(1,n);
68         printf("Case %d: %lld\\n",cas,f[n]);
69     }
70     return 0;
71 }

 

以上是关于bzoj3963[WF2011]MachineWorks cdq分治+斜率优化dp的主要内容,如果未能解决你的问题,请参考以下文章

●BZOJ 3963 [WF2011]MachineWorks

BZOJ 3963: [WF2011]MachineWorks [CDQ分治 斜率优化DP]

BZOJ3958[WF2011]Mummy Madness 二分+扫描线+线段树

BZOJ 2673[Wf2011]Chips Challenge

bzoj3963--cdq分治

BZOJ 4078: [Wf2014]Metal Processing Plant