bzoj3963[WF2011]MachineWorks cdq分治+斜率优化dp
3963: [WF2011]MachineWorks
Time Limit: 30 Sec Memory Limit: 256 MB
Submit: 270 Solved: 80
[Submit][Status][Discuss]
Description
Input
Output
Sample Input
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
转移并不满足单调啊,所以需要用cdq分治来维护单调的转移
也可以写splay
这篇博客写得不错http://www.cnblogs.com/zj75211/p/8148800.html
#include<bits/stdc++.h>
#define N 100050
#define ll long long
#define inf 2147483647
using namespace std;
int n,m,D,s[N];ll f[N];
struct date{int d,p,r,g,id;}q[N],a[N];
//struct point{int x;ll y;}p[N],b[N];
bool cmp1(date a,date b){return a.d<b.d;}
ll X(int i){
return q[i].r+f[q[i].id]-q[i].p-(ll)q[i].g*(q[i].d+);
}
int Y(int i){
return q[i].g;
} double G(int i,int j){
// if(Y(i)==Y(j))return (double)(X(i)>X(j)?inf:-inf);
return (double)(X(i)-X(j))/(double)(Y(i)-Y(j));
}
void solve(int l,int r){
if(l>r)return;
if(l==r){
f[l]=max(f[l],f[l-]);
return;
}
int mid=(l+r)>>;
int p1=l,p2=mid+;
//for(int i=l;i<=r;i++)q[i]=a[i];
solve(l,mid);int tp=,h=;
for(int i=l;i<=mid;i++){
if(q[i].p>f[q[i].id])continue;
if(Y(i)==Y(s[tp])&&tp){
if(X(i)>X(s[tp]))tp--;
else continue;
}
while(tp>&&G(s[tp],s[tp-])<=G(i,s[tp]))tp--;
s[++tp]=i;
}
for(int i=mid+;i<=r;i++){
while(h<tp&&G(s[h],s[h+])>=-q[i].d)h++;
f[i]=max(f[i],f[i-]);
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-));
}
solve(mid+,r);
p1=l;p2=mid+;
for(int i=l;i<=r;i++){
if(p1<=mid&&q[p1].g<=q[p2].g)a[i]=q[p1++];
else if(p2<=r)a[i]=q[p2++];
else a[i]=q[p1++];
}
for(int i=l;i<=r;i++)q[i]=a[i];
}
int main(){
int cas=;
while(scanf("%d%d%d",&n,&m,&D)!=EOF){
cas++;
if(!n&&!m&&!D)break;
memset(f,-,sizeof(f));
f[]=f[]=m;int d,p,r,g;
for(int i=;i<=n;i++)
scanf("%d%d%d%d",&d,&p,&r,&g),
q[i]=(date){d,p,r,g,i};
q[++n]=(date){D+,,,,n};
sort(q+,q++n,cmp1);
for(int i=;i<=n;i++)q[i].id=i;
solve(,n);
printf("Case %d: %lld\n",cas,f[n]);
}
return ;
}
bzoj3963[WF2011]MachineWorks cdq分治+斜率优化dp的更多相关文章
- BZOJ 3963 HDU3842 [WF2011]MachineWorks cdq分治 斜率优化 dp
http://acm.hdu.edu.cn/showproblem.php?pid=3842 写的check函数里写的<但是应该是<=,调了一下午,我是个zz. 就是普通的斜率优化因为有两 ...
- BZOJ 3963: [WF2011]MachineWorks [CDQ分治 斜率优化DP]
传送门 当然了WF的题uva hdu上也有 你的公司获得了一个厂房N天的使用权和一笔启动资金,你打算在这N天里租借机器进行生产来获得收益.可以租借的机器有M台.每台机器有四个参数D,P,R,G.你可以 ...
- 【BZOJ3963】[WF2011]MachineWorks cdq分治+斜率优化
[BZOJ3963][WF2011]MachineWorks Description 你是任意性复杂机器公司(Arbitrarily Complex Machines, ACM)的经理,公司使用更加先 ...
- 【uoj#244】[UER #7]短路 CDQ分治+斜率优化dp
题目描述 给出 $(2n+1)\times (2n+1)$ 个点,点 $(i,j)$ 的权值为 $a[max(|i-n-1|,|j-n-1|)]$ ,找一条从 $(1,1)$ 走到 $(2n+1,2n ...
- BZOJ1492:[NOI2007]货币兑换 (CDQ分治+斜率优化DP | splay动态维护凸包)
BZOJ1492:[NOI2007]货币兑换 题目传送门 [问题描述] 小Y最近在一家金券交易所工作.该金券交易所只发行交易两种金券:A纪念券(以下简称A券)和B纪念券(以下简称B券).每个持有金券的 ...
- BZOJ 1492: [NOI2007]货币兑换Cash [CDQ分治 斜率优化DP]
传送门 题意:不想写... 扔链接就跑 好吧我回来了 首先发现每次兑换一定是全部兑换,因为你兑换说明有利可图,是为了后面的某一天两种卷的汇率差别明显而兑换 那么一定拿全利啊,一定比多天的组合好 $f[ ...
- bzoj3672/luogu2305 购票 (运用点分治思想的树上cdq分治+斜率优化dp)
我们都做过一道题(?)货币兑换,是用cdq分治来解决不单调的斜率优化 现在它放到了树上.. 总之先写下来dp方程,$f[i]=min\{f[j]+(dis[i]-dis[j])*p[i]+q[i]\} ...
- bzoj1492[NOI2007]货币兑换Cash cdq分治+斜率优化dp
1492: [NOI2007]货币兑换Cash Time Limit: 5 Sec Memory Limit: 64 MBSubmit: 5541 Solved: 2228[Submit][Sta ...
- [BZOJ1492] [NOI2007] 货币兑换Cash(cdq分治+斜率优化)
[BZOJ1492] [NOI2007] 货币兑换Cash(cdq分治+斜率优化) 题面 分析 dp方程推导 显然,必然存在一种最优的买卖方案满足:每次买进操作使用完所有的人民币:每次卖出操作卖出所有 ...
随机推荐
- java连接jdbc Establishing SSL connection without server's identity verification is not recommended. According to MySQL 5.5.45+, 5.6.26+ and 5.7.6+ requirements SSL connection must be established by defa
conn = DriverManager.getConnection("jdbc:mysql://localhost:3306/jsp_db","root",& ...
- 浅谈 ThreadLocal
有时,你希望将每个线程数据(如用户ID)与线程关联起来.尽管可以使用局部变量来完成此任务,但只能在本地变量存在时才这样做.也可以使用一个实例属性来保存这些数据,但是这样就必须处理线程同步问题.幸运的是 ...
- VMware虚拟机误删除vmdk文件后如何恢复?
故障描述: Dell R710系列服务器(用于VMware虚拟主机),Dell MD 3200系列存储(用于存放虚拟机文件),VMware ESXi 5.5版本,因意外断电,导致某台虚拟机不能正常启动 ...
- c++ 中lambda
C++ 11中的Lambda表达式用于定义并创建匿名的函数对象,以简化编程工作. 1.Lambda表达式完整的声明格式如下: [capture list] (params list) mutable ...
- 【learning】多项式相关(求逆、开根、除法、取模)
(首先要%miskcoo,这位dalao写的博客(这里)实在是太强啦qwq大部分多项式相关的知识都是从这位dalao博客里面学的,下面这篇东西是自己对其博客学习后的一些总结和想法,大部分是按照其博客里 ...
- 算法题丨3Sum
描述 Given an array S of n integers, are there elements a, b, c in S such that a + b + c = 0? Find all ...
- myeclipse的导航器
在myeclipse的导航器下面可以看到编译后的文件目录结构 如何打开导航器试图呢? 窗口->显示视图->导航器 windows->show view->Navigator 这 ...
- 阿里云API网关(6)用户指南(开放 API )
网关指南: https://help.aliyun.com/document_detail/29487.html?spm=5176.doc48835.6.550.23Oqbl 网关控制台: https ...
- hadoop原理
MapReduce工作原理图文详解 前言: 前段时间我们云计算团队一起学习了hadoop相关的知识,大家都积极地做了.学了很多东西,收获颇丰.可是开学后,大家都忙各自的事情,云计算方面的动静都不太 ...
- windows计划任务
前段时间写了一个小工具,实现两个数据库的数据同步. 需求:要求每天的某个时间自动同步数据 功能写好之后,发现了windows一个自带的功能不错,可以实现我的需要,不用写定时器触发了 控制面板-> ...