luogu P6087 [JSOI2015]送礼物 二分 单调队列 决策单调性
LINK:送礼物
原本想了一个 \(nlog^2\)的做法 然后由于线段树常数过大 T到30.
以为这道题卡\(log^2\)没想到真的有神仙写\(log^2\)的过了 是我常数大了 抱歉。
能过的\(log^2\)的做法是看到了一个 决策单调性优化的dp 证明我不会。
不过由此得到的启示是 一些类似或者就是dp的题目 要多往决策单调性上想 大胆猜想 不管证明。
线段树的做法是二分之后 也是维护dp的决策最优性。发现换set很难更换 所以放弃治疗。
正解是一个log的做法:
考虑 M(i,j) 和 m(i,j) 之间的距离关系 容易发现两种 i-j<=L-2 i-j>=L-1.
前者 可以发现 区间大小永远为L 后者 可以发现i和j是区间的两端。
分类讨论 对于第一种情况 可以发现直接找到这样的l和r 然后由于区间大小固定 所以直接计算即可。
考虑后者 一个r可能存在很多的l 先二分一个mid.
可以得到 \(M(i,j)-m(i,j)-i\cdot mid+j\cdot mid>=k\cdot mid\)
考虑到最优解的l和r之间的距离是>=L-1的 且分属两端 那么强制认为 i为最大值 j为最小值 对于每个i找到一个j即可。
这样的好处是 可以发现一定可以扫到那个最优解 且 如果是j不是最小值 或者i不是最大值 那么刚才那个值只会更小不会更大.
至此两种情况讨论清楚了 所以复杂度为nlogn.注意精度要高一点不然会wa.
//#include<bits\stdc++.h>
#include<iostream>
#include<iomanip>
#include<cstdio>
#include<cstring>
#include<string>
#include<ctime>
#include<cmath>
#include<cctype>
#include<cstdlib>
#include<queue>
#include<deque>
#include<stack>
#include<vector>
#include<algorithm>
#include<utility>
#include<bitset>
#include<set>
#include<map>
#define ll long long
#define db double
#define INF 1000000000
#define ldb long double
#define pb push_back
#define put_(x) printf("%d ",x);
#define get(x) x=read()
#define gt(x) scanf("%d",&x)
#define gi(x) scanf("%lf",&x)
#define put(x) printf("%d\n",x)
#define putl(x) printf("%lld\n",x)
#define gc(a) scanf("%s",a+1)
#define rep(p,n,i) for(RE int i=p;i<=n;++i)
#define go(x) for(int i=lin[x],tn=ver[i];i;tn=ver[i=nex[i]])
#define fep(n,p,i) for(RE int i=n;i>=p;--i)
#define pii pair<int,int>
#define mk make_pair
#define RE register
#define P 1000000007
#define F first
#define gf(x) scanf("%lf",&x)
#define pf(x) ((x)*(x))
#define ull unsigned long long
#define ui unsigned
#define EPS 1e-7
#define sq sqrt
#define mod 998244353
#define l(p) t[p].l
#define r(p) t[p].r
#define sum(p) t[p].sum
#define tag(p) t[p].tag
#define zz p<<1
#define yy p<<1|1
using namespace std;
char buf[1<<15],*fs,*ft;
inline char getc()
{
return (fs==ft&&(ft=(fs=buf)+fread(buf,1,1<<15,stdin),fs==ft))?0:*fs++;
}
inline int read()
{
RE int x=0,f=1;RE char ch=getc();
while(ch<'0'||ch>'9'){if(ch=='-')f=-1;ch=getc();}
while(ch>='0'&&ch<='9'){x=x*10+ch-'0';ch=getc();}
return x*f;
}
const int MAXN=50010;
int T,cnt,cnt1;
int n,k,L,R,maxx,l,r;
int a[MAXN],Log[MAXN];
int f[MAXN][20],q[MAXN];
inline int cmp(int x,int y){return a[x]>a[y]?x:y;}
inline int ask(int l,int r)
{
int z=Log[r-l+1];
return cmp(f[l][z],f[r-(1<<z)+1][z]);
}
inline int check(db mid)
{
l=1;r=0;
rep(L,n,i)
{
while(l<=r&&i-q[l]>=R)++l;
while(l<=r&&(i-L+1)*mid-a[i-L+1]>=q[r]*mid-a[q[r]])--r;
q[++r]=i-L+1;
if(a[i]-i*mid-a[q[l]]+q[l]*mid>=k*mid)return 1;
}
l=1;r=0;
rep(L,n,i)
{
while(l<=r&&i-q[l]>=R)++l;
while(l<=r&&a[i-L+1]+(i-L+1)*mid>=a[q[r]]+q[r]*mid)--r;
q[++r]=i-L+1;
if(a[q[l]]+q[l]*mid-a[i]-i*mid>=k*mid)return 1;
}
return 0;
}
int main()
{
freopen("1.in","r",stdin);
get(T);
rep(2,50000,i)Log[i]=Log[i>>1]+1;
while(T--)
{
cnt=cnt1=maxx=0;
get(n);get(k);get(L);get(R);
rep(1,n,i)get(a[i]),f[i][0]=i;
rep(1,Log[n],j)
rep(1,n-(1<<j)+1,i)f[i][j]=cmp(f[i][j-1],f[i+(1<<(j-1))][j-1]);
db ans=0;l=1;r=0;
rep(1,n,i)
{
while(l<=r&&a[i]<=a[q[r]])--r;
while(l<=r&&i-q[l]>=L-1)++l;
q[++r]=i;
if(i>=L-1)ans=max(ans,(a[ask(i-L+2,i)]-a[q[l]])*1.0/(L-1+k));
}
db wl=0,wr=1000;
while(wl+EPS<wr)
{
db mid=(wl+wr)/2;
if(check(mid))wl=mid;
else wr=mid;
}
ans=max(ans,wl);
printf("%.4lf\n",ans);
}
return 0;
}
luogu P6087 [JSOI2015]送礼物 二分 单调队列 决策单调性的更多相关文章
- P6087 [JSOI2015]送礼物 01分数规划+单调队列+ST表
P6087 [JSOI2015]送礼物 01分数规划+单调队列+ST表 题目背景 \(JYY\) 和 \(CX\) 的结婚纪念日即将到来,\(JYY\) 来到萌萌开的礼品店选购纪念礼物. 萌萌的礼品店 ...
- BZOJ.1010.[HNOI2008]玩具装箱toy(DP 斜率优化/单调队列 决策单调性)
题目链接 斜率优化 不说了 网上很多 这的比较详细->Click Here or Here //1700kb 60ms #include<cstdio> #include<cc ...
- BZOJ_4476_[Jsoi2015]送礼物_01分数规划+单调队列
BZOJ_4476_[Jsoi2015]送礼物_01分数规划+单调队列 Description JYY和CX的结婚纪念日即将到来,JYY来到萌萌开的礼品店选购纪念礼物. 萌萌的礼品店很神奇,所有出售的 ...
- [BZOJ4476] [JSOI2015] 送礼物 (01分数规划+ST表)
[BZOJ4476] [JSOI2015] 送礼物 (01分数规划+ST表) 题面 给出n,k,l,r和序列a,要求从a中选一段连续的区间[i,j]出来,使得M(i,j)-m(i,j)/(j-i+k) ...
- 【BZOJ4476】[Jsoi2015]送礼物 分数规划+RMQ
[BZOJ4476][Jsoi2015]送礼物 Description JYY和CX的结婚纪念日即将到来,JYY来到萌萌开的礼品店选购纪念礼物.萌萌的礼品店很神奇,所有出售的礼物都按照特定的顺序都排成 ...
- 两种解法-树形dp+二分+单调队列(或RMQ)-hdu-4123-Bob’s Race
题目链接: http://acm.hdu.edu.cn/showproblem.php?pid=4123 题目大意: 给一棵树,n个节点,每条边有个权值,从每个点i出发有个不经过自己走过的点的最远距离 ...
- [BZOJ1044][HAOI2008]木棍分割 二分 + 单调队列优化dp + 滚动数组优化dp
Description 有n根木棍, 第i根木棍的长度为Li,n根木棍依次连结了一起, 总共有n-1个连接处. 现在允许你最多砍断m个连接处, 砍完后n根木棍被分成了很多段,要求满足总长度最大的一段长 ...
- NOIP模拟 最佳序列 - 二分 + 单调队列
题意: 各一个n(\(\le 20000\))的序列,定义纯洁序列为长度len满足\(L \le len \le R\)的序列,纯洁值为某一纯洁序列的平局值,输出所有纯洁序列中最大平均值. 分析: 二 ...
- [USACO12MAR]花盆 二分 单调队列
[USACO12MAR]花盆 二分 单调队列 存在一个长度为\(x\)的区间\([l,r]\),使得区间中最大值与最小值差至少为\(w\),求这个最小的\(x\) \(n\le 100000\),\( ...
随机推荐
- # Mysql常用函数总结(一)
Mysql常用函数总结(一) 博客已搬家,更多内容查看https://liangyongrui.github.io/ 遇到什么总结什么 DATE_SUB(date,INTERVAL expr type ...
- NanoPi arm架构下的程序 ./ 运行黑屏 Qt环境可运行
首先之所以QtCreator环境下可直接运行,但是在终端下 ./ 则不能运行(黑屏但是不报错),判断肯定不是程序或者是库的问题.于是猜想是环境问题,即终端环境与QtCreator环境不同. 然后就查看 ...
- 【转】Web端测试点整理
版权声明:本文为博主原创文章,遵循 CC 4.0 by-sa 版权协议,转载请附上原文出处链接和本声明.本文链接:https://blog.csdn.net/qq_24373725/article/d ...
- 控制流程之while循环, for循环
条件循环:while,语法如下 while 条件: # 循环体 # 如果条件为真,那么循环体则执行,执行完毕后再次循环,重新判断条件... # 如果条件为假,那么循环体不执行,循环终止死循环 基本使用 ...
- python 装饰器(一):装饰器基础(一)装饰器形式,何时执行
简介 装饰器是可调用的对象,其参数是另一个函数(被装饰的函数). 装饰器可能会处理被装饰的函数,然后把它返回,或者将其替换成另一个函数或可调用对象. 形式 假如有个名为 decorate 的装饰器: ...
- How to start MySQL on Linux
启动MySQL数据库 service mysql start 查看MySQL进程 ps -ef |grep mysql 查看MySQL端口号 cd /etc/init.d/ netstat -atnp ...
- uPDF-功能强大的PDF文件处理小工具
前几天因为工作原因,需要将一个PDF压缩一下. 网上找了半天,要么收费,要么就是转换的质量不太好.论坛也找到一些破解的软件,但是总有点不太合适,有些功能还挺复杂.也有些在线转换的,又考虑到自己较为隐私 ...
- sql多表语句
多条件查询条件判空 最优写法 3三表带条件查询
- 什么是控制反转(IoC)?什么是依赖注入(DI)?以及实现原理
IoC不是一种技术,只是一种思想,一个重要的面向对象编程的法则,它能指导我们如何设计出松耦合.更优良的程序.传统应用程序都是由我们在类内部主动创建依赖对象,从而导致类与类之间高耦合,难于测试:有了 ...
- Burp Suite Extender Module - 扩展模块
模块功能: 在扩展模块可以通过使用自定义代码,进行Burp 的自定义操作. 1. Burp Extensions页面 2. BApp Store中可以购买和安装别人写好的扩展功能 3. 在APIs界面 ...