[SDOI2012]任务安排 BZOJ2726 斜率优化+二分查找
网上的题解...状态就没有一个和我一样的...这让我有些无从下手...
分析:
我们考虑,正常的斜率优化满足x(i)单调递增,k(i)单调递增,那么我们就可以只用维护一个单调队列满足对于当前的x(i)有最小值即可,因为x(i)满足单调递增。这样的话,我们就可以维护一个单调队列让队首元首最小。而这道题,可以发现有部分数据满足x(i)单调递增,那么直接裸上就可以,但是由于时间有负数,所以x(i)并不满足单调性。但是由于k(i)仍然满足单调性,因此,我们依然可以发现更新x(i)的时候满足斜率优化的性质,也就是我们已经可以将已经完全被覆盖的直线忽略。(即:我们依然可以维护一个大凹包)
考虑n^2的DP方程:f[i]表示当前的最后一段区间在以i为终点时的最小花费。
转移:f[i]=min{f[j]+(sum1[i]-sum1[j]+m)*(sum2[n]-sum2[j])};
之后,我们发现在我们维护的队列中,对于相同的x(i),不同的j,更新的f[i]具有单峰性,那么导函数具有单调性,二分导函数,找到最大值,之后更新答案,剩下的就是斜率优化的常规操作了。
附上代码:
#include <cstdio>
#include <algorithm>
#include <cmath>
#include <cstring>
#include <cstdlib>
#include <queue>
#include <iostream>
using namespace std;
#define N 300005
#define ll long long
#define K(x) (-sum2[x])
#define B(x) (f[x]-sum1[x]*sum+sum1[x]*sum2[x]-m*sum2[x])
#define Y(x,y) (K(y)*sum1[x]+B(y))
int n,q[N];
ll f[N],sum1[N],sum2[N],sum,m;
bool cmp(int i,int j,int k)
{
ll t1=(K(j)-K(k))*(B(i)-B(k));
ll t2=(K(i)-K(k))*(B(j)-B(k));
return t1<=t2;
}
bool check(int p1,int p2,int i)
{
return /*f[p2]-sum2[p2]*(sum1[i]+m)*/Y(i,p2)</*f[p1]-(sum1[i]+m)*sum2[p1]*/Y(i,p1);
}
int main()
{
scanf("%d%lld",&n,&m);
for(int i=1,x,y;i<=n;i++)
{
scanf("%d%d",&x,&y);
sum1[i]=sum1[i-1]+x;
sum2[i]=sum2[i-1]+y;
}
sum=sum2[n];int t=0;memset(f,0x3f,sizeof(f));f[0]=0;
for(int i=1;i<=n;i++)
{
int l=0,r=t;
while(l<r)
{
int mid=(l+r)>>1;
if(check(q[mid],q[mid+1],i))l=mid+1;
else r=mid;
}
// printf("%d\n",q[l]);
f[i]=Y(i,q[l])+sum*sum1[i]+m*sum;
while(t&&cmp(q[t-1],q[t],i))t--;
q[++t]=i;
}
printf("%lld\n",f[n]);return 0;
}
[SDOI2012]任务安排 BZOJ2726 斜率优化+二分查找的更多相关文章
- BZOJ_2726_[SDOI2012]任务安排_斜率优化+二分
BZOJ_2726_[SDOI2012]任务安排_斜率优化+二分 Description 机器上有N个需要处理的任务,它们构成了一个序列.这些任务被标号为1到N,因此序列的排列为1,2,3...N.这 ...
- BZOJ2726 [SDOI2012]任务安排 【斜率优化 + cdq分治】
题目 机器上有N个需要处理的任务,它们构成了一个序列.这些任务被标号为1到N,因此序列的排列为1,2,3...N.这N个任务被分成若干批,每批包含相邻的若干任务.从时刻0开始,这些任务被分批加工,第i ...
- 2018.09.05 bzoj2726: [SDOI2012]任务安排(斜率优化dp+二分)
传送门 跟Ti" role="presentation" style="position: relative;">TiTi为正数的时候差不多. ...
- BZOJ.2726.[SDOI2012]任务安排(DP 斜率优化)
题目链接 数据范围在这:https://lydsy.com/JudgeOnline/wttl/thread.php?tid=613, 另外是\(n\leq3\times10^5\). 用\(t_i\) ...
- hdu2993之斜率dp+二分查找
MAX Average Problem Time Limit: 2000/1000 MS (Java/Others) Memory Limit: 32768/32768 K (Java/Othe ...
- BZOJ2726【SDOI2012】任务安排(斜率优化Dp+二分查找)
由题目条件显然可以得到状态 f[i][j] 表示以 i 为结尾且 i 后作为断点,共做了 j 次分组的最小代价. 因此转移变得很显然:f[i][j]=min{f[k][j-1]+(s×j+sumT[i ...
- [BZOJ2726][SDOI2012]任务安排(DP+凸壳二分)
2726: [SDOI2012]任务安排 Time Limit: 10 Sec Memory Limit: 128 MBSubmit: 1580 Solved: 466[Submit][Statu ...
- 洛谷P2365/5785 任务安排 题解 斜率优化DP
任务安排1(小数据):https://www.luogu.com.cn/problem/P2365 任务安排2(大数据):https://www.luogu.com.cn/problem/P5785 ...
- 2018.09.05 任务安排(斜率优化dp)
描述 这道题目说的是,给出了n项必须按照顺序完成的任务,每项任务有它需要占用机器的时间和价值.现在我们有一台机器可以使用,它每次可以完成一批任务,完成这批任务所需的时间为一个启动机器的时间S加上所有任 ...
随机推荐
- requireJS基本概念及使用流程(2)
上一篇我们一起研究了研究requireJS,这一篇我们来说一说requireJS具体的使用过程 其实很简单的,我总结了总结就是分为四步走 第一步:在页面中引入requireJS并且引入入口文件 第二步 ...
- SD从零开始57-58,第三方订单处理,跨公司销售
[原创] SD从零开始57 第三方订单处理流程 第三方订单处理的流程Processes for Third-Party Order Processing 客户的采购订单首先在你公司的一个销售组织作为一 ...
- 关于kernel-devel、kernel-header和kernel src的区别
A kernel-header package would contain 'header files' needed by some applications which would be buil ...
- Java语法基础(四)----循环结构语句
一.循环结构: 循环语句可以在满足循环条件的情况下,反复执行某一段代码,这段被重复执行的代码被称为循环体语句,当反复执行这个循环体时,需要在合适的时候把循环判断条件修改为false,从而结束循环,否则 ...
- Apache2启动错误Could not reliably determine the server's fully qualified domain name
错误情况: AH00558: apache2: Could not reliably determine the server's fully qualified domain name, using ...
- Safari中的Date
在js中处理Date时,发现Safari和其他浏览器的支持方式不一致 1.例如:2017-01-01 12:00:00 在其他浏览器中,使用这个格式的字符串进行new Date操作没有问题,但是在Sa ...
- 红帽7中firewall常用指令
1.端口管理 (1)列出DMZ区域开放的端口 ~]#firewall-cmd --zone=dmz --list-ports (2)8080端口加入dmz区 ~]#firewall-cmd --zon ...
- centos7 安装ldap
ldap首先我们要知道这个ldap的概念, LDAP是轻量目录访问协议(Lightweight Directory Access Protocol)的缩写 目录是一个为查询.浏览和搜索而优化的专业分布 ...
- 通过C/C++基于http下载文件
简介 Windows系统如何通过C/C++下载互联网上的文件呢?这里笔者给大家演示一个最简单的方法,利用Windows提供的urlmon库,可以快速实现文件下载的简单实例. 注:本文内容部分参考了&l ...
- 使用let声明变量的理解
先看阮大神的[ECMAScript 6 入门]中关于这一部分的描述 var a = []; for (let i = 0; i < 10; i++) { a[i] = function () { ...