●BZOJ 3963 [WF2011]MachineWorks
题链:
http://www.lydsy.com/JudgeOnline/problem.php?id=3963
题解:
斜率优化DP,CDQ分治。
先按时间排序。(规定以下内容的第i台机器的卖出时间D[i]大于第i-1台机器的卖出时间D[i-1])
定义DP[i]表示在在第i台机器可以交易的那天,卖出所有机器后能够得到的最大收益,(最后答案是DP[all_day+1])
转移显然:
$DP[i]=DP[i-1]$
$DP[i]=max(DP[j]+(D_i-D_j-1)*G_j-P_j+R_j) (j<i且DP[j]>=P[j])$
令$Y_j=DP[j]-(D_j+1)*G_j-P_j+R_j$,如果存在两个转移点k,j且G[k]<G[j],假设j点优于k点
那么 $Y_j-Y_k>-D_i(G_j-G_k)$
$\quad\quad\frac{Y_j-Y_k}{G_j-G_k}>-D_i$
那么得到结论,如果 G[k]<G[j],且Slope(j,k)>-D[i]的话,则j点优于k点。
同时如果存在三个转移来源点:k,j,i,满足G[k]<G[j]<G[i],
同时Slope(i,j)>Slope(j,k),则j点无效。
所以对于每个来源点二元组(G[j],Y[j]),只需要在平面上维护一个上凸壳即可。
但是G不单调,所以用CDQ分治。
对于分治的每一层l~r,先递归处理左边l~mid,然后把左边按G从小到大排序,并维护好上凸壳。
由于D单调递增,所以遍历一边凸壳以及右边mid+1~r进行贡献就好。
代码:
#include<cstdio>
#include<cstring>
#include<iostream>
#include<algorithm>
#define MAXN 100050
#define ll long long
using namespace std;
ll Y[MAXN],G[MAXN],D[MAXN],P[MAXN],R[MAXN],DP[MAXN];
int h[MAXN],N,C,T,S;
bool cmp(int i,int j){
return D[i]<D[j];
}
struct Moque{
int q[MAXN],l,r;
#define Slope(i,j) (1.0*(Y[i]-Y[j])/(G[i]-G[j]))
void Reset(){l=1; r=0;}
void Push(int i){
if(l<=r&&G[i]==G[q[r]])
{if(Y[i]>Y[q[r]]) r--; else return;}
while(l+1<=r&&Slope(i,q[r])>Slope(q[r],q[r-1])) r--;
q[++r]=i;
}
int Query(int i){
while(l+1<=r&&Slope(q[l+1],q[l])>-D[i]) l++;
return q[l];
}
}Q;
void solve(int l,int r){
static int tmp[MAXN],cl,cr,p; static ll MAXDP;
if(l==r) return (void)(Y[h[l]]=DP[h[l]]-(D[h[l]]+1)*G[h[l]]-P[h[l]]+R[h[l]]);
int mid=(l+r)>>1;
solve(l,mid);//之后,左边G单调
Q.Reset(); MAXDP=0;
for(int i=l;i<=mid;i++){
MAXDP=max(MAXDP,DP[h[i]]);
if(DP[h[i]]>=P[h[i]]) Q.Push(h[i]);//上凸壳斜率单减
}
//Di单增,-Di单减,正序枚举就好。
for(int i=mid+1,j;i<=r;i++){
j=Q.Query(h[i]);
DP[h[i]]=max(DP[h[i]],MAXDP);
DP[h[i]]=max(DP[h[i]],DP[j]+(D[h[i]]-D[j]-1)*G[j]-P[j]+R[j]);
}
solve(mid+1,r);
cl=l; cr=mid+1; p=l;
while(cl<=mid||cr<=r){
if(cl>mid) tmp[p]=h[cr],cr++;
else if(cr>r||G[h[cl]]<G[h[cr]]) tmp[p]=h[cl],cl++;
else tmp[p]=h[cr],cr++; p++;
}
for(int i=l;i<=r;i++) h[i]=tmp[i];
}
int main(){
int Cas=0;
while(1){
memset(DP,0,sizeof(DP));
scanf("%d%d%d",&N,&S,&T);
if(N==0&&DP[1]==0&&T==0) break;
memset(Y,0,sizeof(Y));
for(int i=1,a,b,c,d;i<=N;i++){
scanf("%d%d%d%d",&a,&b,&c,&d);
h[i]=i; D[i]=a; P[i]=b; R[i]=c; G[i]=d;
}
N++; D[N]=T+1; h[N]=N;
sort(h+1,h+N+1,cmp);
DP[h[1]]=S; solve(1,N);
printf("Case %d: %lld\n",++Cas,DP[N]);
}
return 0;
}
●BZOJ 3963 [WF2011]MachineWorks的更多相关文章
- HDU 3824/ BZOJ 3963 [WF2011]MachineWorks (斜率优化DP+CDQ分治维护凸包)
题面 BZOJ传送门(中文题面但是权限题) HDU传送门(英文题面) 分析 定义f[i]f[i]f[i]表示在iii时间(离散化之后)卖出手上的机器的最大收益.转移方程式比较好写f[i]=max{f[ ...
- BZOJ 3963: [WF2011]MachineWorks [CDQ分治 斜率优化DP]
传送门 当然了WF的题uva hdu上也有 你的公司获得了一个厂房N天的使用权和一笔启动资金,你打算在这N天里租借机器进行生产来获得收益.可以租借的机器有M台.每台机器有四个参数D,P,R,G.你可以 ...
- BZOJ 3963: [WF2011]MachineWorks 斜率优化 + splay动态维护凸包
Description 你是任意性复杂机器公司(Arbitrarily Complex Machines, ACM)的经理,公司使用更加先进的机械设备生产先进的机器.原来的那一台生产机器已经坏了,所以 ...
- 【BZOJ】3963: [WF2011]MachineWorks
[题意]给定n台在时间di可以买入的机器,pi买入,可在任意时间ri卖出,买入和卖出之间的持有时间每天产生gi金钱,任意时间至多持有一台机器.给定初始钱数c和总天数T,求最大收益.n<=10^5 ...
- Solution -「WF2011」「BZOJ #3963」MachineWorks
\(\mathcal{Description}\) Link. 给定你初始拥有的钱数 \(C\) 以及 \(N\) 台机器的属性,第 \(i\) 台有属性 \((d_i,p_i,r_i,g_i ...
- bzoj3963[WF2011]MachineWorks cdq分治+斜率优化dp
3963: [WF2011]MachineWorks Time Limit: 30 Sec Memory Limit: 256 MBSubmit: 270 Solved: 80[Submit][S ...
- BZOJ3963: [WF2011]MachineWorks 【CDQ+斜率优化DP】*
BZOJ3963: [WF2011]MachineWorks Description 你是任意性复杂机器公司(Arbitrarily Complex Machines, ACM)的经理,公司使用更加先 ...
- 【BZOJ3963】[WF2011]MachineWorks cdq分治+斜率优化
[BZOJ3963][WF2011]MachineWorks Description 你是任意性复杂机器公司(Arbitrarily Complex Machines, ACM)的经理,公司使用更加先 ...
- BZOJ 3963 HDU3842 [WF2011]MachineWorks cdq分治 斜率优化 dp
http://acm.hdu.edu.cn/showproblem.php?pid=3842 写的check函数里写的<但是应该是<=,调了一下午,我是个zz. 就是普通的斜率优化因为有两 ...
随机推荐
- MySQL的小Tips
交集和差集 MySQL中没有这两个运算,但是有并集运算,所以可以利用这个来间接实现. 差集: SELECT ID FROM ( SELECT DISTINCT A.AID AS ID FROM TAB ...
- OSI七层协议模型、TCP/IP四层模型学习笔记
1. OSI七层和TCP/IP四层的关系 1.1 OSI引入了服务.接口.协议.分层的概念,TCP/IP借鉴了OSI的这些概念建立TCP/IP模型. 1.2 OSI先有模型,后有协议,先有标准,后进行 ...
- OO第一次总结
第一次作业: 第一次作业的指导书发下来之后我按着上面的步骤一步一步的做了之后发现项目拉下来了,怎么开始码代码呢...然后在舍友的帮助下才知道怎么建包建类,然后对Java的语法又不是很了解,于是就先把C ...
- JAVA_SE基础——69.Date类
package cn.itcast.other; import java.text.ParseException; import java.text.SimpleDateFormat; import ...
- Python之旅_计算机基础入门
一.计算机基础 1.Python是编程语言 语言:一种事物与另一种事物沟通的介质. 编程语言:程序员与计算机沟通的介质. 什么是编程:程序员用编程语言把自己的逻辑思想下来,编程的结果就是一堆文件. 为 ...
- ajax和jquery使用技巧
1.使用ajax的方法的时候可以使用u方法来获取连接,这样更加安全:alert弹窗的时候需要单引号双引号火狐浏览器会报错!
- docker实践3
我的docker学习笔记3 $docker run ubuntu echo'hello world' $docker run -i -t ubuntu /bin/bash #ps -ef #exi ...
- Linux后台运行命令 nohup command > myout.file 2>&1
Linux命令后台运行 转自北国的雨,谢谢:http://www.cnblogs.com/lwm-1988/archive/2011/08/20/2147299.html 有两种方式:1. comma ...
- j2ee基础(1)servlet的生命周期
Servlet的生命周期 Servlet 生命周期规定了 Servlet 如何被加载.实例化.初始化. 处理客户端请求,以及何时结束服务. 该生命周期可以通过 javax.servlet.Servle ...
- Error loading MySQLdb module: No module named 'MySQLdb'----------- django成功连接mysql数据库的方法
在进行django学习过程中,尝试使用框架连接mysql数据库,启动服务器的时候经常遇到Error loading MySQLdb module: No module named 'MySQLdb' ...