Function:凸包,单调栈,题意转化,单峰函数三分,离线处理
很难啊啊啊!!!
bzoj5380原题,应该可以粘题面。
问题转换:
有一个n列1e9行的矩阵,每一列上都写着相同的数字Ai。
你从位置(x,y)出发每一步可以向左上方或左方走一步,最后走到第一行。
要求最小化路径上的总权值。
首先题意转化就让我挂了。。。
然后题解里一个显然的结论让我又挂了一回:最优决策是先往左上走几步,在往上一直走。
证明比较简单。因为如果你往上走了几步后再往左上走,那么一定不如先往左上走,不然就不是最优决策。
于是可以枚举转弯位置。(27%暴力)
#include<cstdio>
#include<iostream>
#define int long long
int n,q,a[],sum[],x,y,ans;
main(){
scanf("%lld",&n);
for(int i=;i<=n;++i)scanf("%lld",&a[i]),sum[i]=sum[i-]+a[i];
scanf("%lld",&q);
while(q--){
scanf("%lld%lld",&x,&y);
ans=x*a[y];
for(int i=;i<y;++i)if(y-i<x)ans=std::min(ans,sum[y]-sum[i]+a[i]*(x+i-y));
printf("%lld\n",ans);
}
}
可以发现,对于出发点(x,i),它在第j列转弯的话,费用是一个以x为自变量的一次函数。
我们同时把所有的直线画出来,然后对x时的直线们取最低点即可。
怎么维护这些直线呢?
首先可以证明,对于某一列上的任意x,其可能的最优解直线的A一定是递增的。
那么我们就可以离线所有询问,依次处理每一列就好,始终维护
感性理解的话就是你斜着往前走一定会找到一个A更小的再往上走,如果你跨过了更小的A走到更大的A上再往上走当然不优。。。
具体证明的话没时间写啦。。。推荐blog:https://blog.csdn.net/Rose_max/article/details/82250809
这样的话我们维护一个单调栈,栈中的直线按A递增,加入新的直线时弹栈把更大的A弹掉就好了。
但仅仅这样还不够,并没有把所有没用的直线弹掉。
考虑栈里的两条直线:
栈顶:y=x
倒数第二条:y=0.5x+0.5
现在要加入直线:y=2x-4
这三条直线满足A单调递增,那么可以把直线直接加入了?
并不是,在坐标系上画出这三条直线,我们会发现没有任何一个最小值在直线y=x上
在点(3,2)以后y=0.5x+0.5最小,以前是y=2x-4
这就是维护凸壳的事情了。。。我们通过判断新加入的直线与栈顶的后两条直线的交点位置来判定栈顶的直线还有没有用。
如果栈顶直线的交点横坐标更大,那么栈顶直线就没有用了。
及时弹栈之后我们就有了那个凸壳,现在剩下的问题就是对于已知x怎么求它在凸壳上的值了。
可以证明,从最优的那条直线往栈两边扫,费用都是单调递增的。
所以费用关于斜率是个单谷函数。三分求出所在直线即可。
#include<cstdio>
#include<algorithm>
using namespace std;
#define int long long
struct qs{int x,y,id;friend bool operator<(qs a,qs b){return a.y<b.y;}}q[];
int n,a[],Q,b[],sum[],ans[],sta[];
int cal(int p,int x){return a[p]*x+b[p];}
double jiao_dian(int p,int q){return (1.0*b[q]-b[p])/(a[p]-a[q]);}
signed main(){
scanf("%lld",&n);
for(int i=;i<=n;++i)scanf("%lld",&a[i]),sum[i]=sum[i-]+a[i];
scanf("%lld",&Q);
for(int i=;i<=Q;++i)scanf("%lld%lld",&q[i].x,&q[i].y),q[i].id=i;
sort(q+,q++Q);
for(int i=,alq=,top=;i<=n;++i){
b[i]=-sum[i-]+a[i]*(i-);//含义是:计算时的定值的sum做差,然后多走了i-1步(本来可以走到前面,但是在这里选择了往上走)
while(top&&a[sta[top]]>=a[i])top--;//栈内斜率一定递增否则不优
while(top>&&jiao_dian(sta[top-],i)<jiao_dian(sta[top],i))top--;//如果栈顶的直线在任意的x上都不是最又决策那么就pop掉
sta[++top]=i;
while(alq<=Q&&q[alq].y==i){
int X=q[alq].x,Y=q[alq].y,l=lower_bound(sta+,sta+top+,Y-X+)-sta,r=top-;//找到可行的栈内区间
while(l<r)if(cal(sta[l+r>>],X-Y)>cal(sta[(l+r>>)+],X-Y))l=(l+r>>)+;else r=l+r>>;//三分单谷,X-Y即为至少要走的步数,剩下的步数在b里面
while(l<top&&cal(sta[l],X-Y)>=cal(sta[l+],X-Y))l++;//后延直到最优解
ans[q[alq].id]=cal(sta[l],X-Y)+sum[i];alq++;//要加上斜着走的费用
}
}
for(int i=;i<=Q;++i)printf("%lld\n",ans[i]);
}
Function:凸包,单调栈,题意转化,单峰函数三分,离线处理的更多相关文章
- [CSP-S模拟测试]:导弹袭击(数学+凸包+单调栈)
题目背景 $Guess$准备向敌军阵地发起进攻了!$Guess$的武器是自动制导导弹.然而在机房是不允许游戏的,所以班长$XZY$对游戏界面进行了降维打击,结果... 题目描述 众所周知,环境因素对导 ...
- 「洛谷5300」「GXOI/GZOI2019」与或和【单调栈+二进制转化】
题目链接 [洛谷传送门] 题解 按位处理. 把每一位对应的图都处理出来 然后单调栈处理一下就好了. \(and\)操作处理全\(1\). \(or\)操作处理全\(0\). 代码 #include & ...
- HDU 5875 H - Function 用单调栈水过了
http://acm.hdu.edu.cn/showproblem.php?pid=5875 单调栈,预处理to[i]表示第一个比a[i]小的数字,一直跳就可以. 这题是数据水而已. 这里学习下单调栈 ...
- Poj 3250 单调栈
1.Poj 3250 Bad Hair Day 2.链接:http://poj.org/problem?id=3250 3.总结:单调栈 题意:n头牛,当i>j,j在i的右边并且i与j之间的所 ...
- 2019牛客暑期多校训练营(第一场) - A - Equivalent Prefixes - 单调栈
A - Equivalent Prefixes - 单调栈 题意:给定两个n个元素的数组a,b,它们的前p个元素构成的数组是"等价"的,求p的最大值."等价"的 ...
- [Agc005D/At2060] Minimum Sum - 单调栈
鉴于早上那题让我怀疑单调栈白学,特意来复习下单调栈 题意 考虑按照每个元素对答案的贡献来统计,那么我们只需要找到每个元素左边右边第一个比它小的就可 这题给的又是排列,简直不能再良心 #include ...
- 2016 大连网赛---Function(单调栈)
题目链接 http://acm.split.hdu.edu.cn/showproblem.php?pid=5875 Problem Description The shorter, the simpl ...
- CF535E Tavas and Pashmaks 单调栈、凸包
传送门 题意:有一场比赛,$N$个人参加.每个人有两种参数$a,b$,如果存在正实数$A,B$使得$\frac{A}{a_i} + \frac{B}{b_i}$在$i=x$处取得最大值(可以有多个最大 ...
- HDU 5875 Function (线段树+gcd / 单调栈)
题意:给你一串数a再给你一些区间(lef,rig),求出a[lef]%a[lef+1]...%a[rig] 题解:我们可以发现数字a对数字b取模时:如果a<b,则等于原数,否则a会变小至少一半. ...
随机推荐
- IDEA 学习笔记之 Java项目开发
Java项目开发: 新建模块: 添加JDK: 导入本地Jars: 从远程Maven仓库下载: 创建package: 新建类/接口/枚举等: 字体太小,改字体: Duplicate Scheme 修改编 ...
- Tomcat线程参数maxThreads、acceptCount
一.配置Tomcat/conf/server.xml修改配置 <Connector port="8080" protocol="org.apache.coyote. ...
- 从零开始的 phpstorm+wamp 组合下的debug环境搭建(纯小白向)
本文主要是为了帮自己记住每次重装系统后需要干点啥,如果能帮到你,烦请给个好评 环境说明: 1. windows10 64bit 2. wampservers 3.0.6(x86) apache2.4. ...
- Redis未授权访问写Webshell和公私钥认证获取root权限
0x01 什么是Redis未授权访问漏洞 Redis 默认情况下,会绑定在 0.0.0.0:,如果没有进行采用相关的策略,比如添加防火墙规则避免其他非信任来源 ip 访问等,这样将会将 Redis 服 ...
- Webshell免杀绕过waf
转自圈子404师傅 0x01 前言# 尽最大努力在一文中让大家掌握一些有用的WEBSHELL免杀技巧 0x02 目录# 关于eval 于 assert 字符串变形 定义函数绕过 回调函数 回调函数变形 ...
- 《锋利的jQuery》学习总结
通过对<锋利的jQuery>(第二版)一书的学习,发现此书讲解通俗易懂,是学习jQuery的一本很好的指导书,特作如下总结.此书主要讲解了jQuery的常用操作,包括认识jQuery,jQ ...
- JavaScript函数总结—越努力,越幸运!
JavaScript 函数总结 JavaScript为web的编程脚本语言. JavaScript由三部分组成:emc(语法) dom(文档对象模型) bom(浏览器对象模型). [函数的定义] 1. ...
- 解决seajs ie8 对象不支持charAt 属性。
在使用 seajs做项目,今天偶然发现在ie9以下的ie版本会 报出 对象不支持charAt 属性.刚开始还以为是自己写的js部分出了问题,经过几个小时的奋战.最终找到了其根源.在sea-debug. ...
- [JOYOI1307] 联络员
题目限制 时间限制 内存限制 评测方式 题目来源 1000ms 131072KiB 标准比较器 Local 题目描述 Tyvj已经一岁了,网站也由最初的几个用户增加到了上万个用户,随着Tyvj网站的逐 ...
- maven子项目导出成jar包及运行
第一步:选这idea右侧栏的maven projects 第二步:选中需要打包成jar包的项目下的lifecycle 第三步:选中package 第四步:点击开始导出 第五步:使用winRAR打开ja ...