[BZOJ4476] [JSOI2015] 送礼物 (01分数规划+ST表)

题面

给出n,k,l,r和序列a,要求从a中选一段连续的区间[i,j]出来,使得M(i,j)-m(i,j)/(j-i+k)最大,且[i,j]长度在区间[l,r]内。

M(i,j)表示[i,j]中的最大值,m(i,j)表示[i,j]中的最小值。

Ai< =10^8,N,K< = 50,000

分析

看到\(\frac{max(i,j)-min(i,j)}{j-i+k}\)的形式,我们敏感的想到01分数规划。考虑最终答案我们取的最大值和最小值在序列中的位置。

1.如果最大值和最小值在序列中的距离<L,那么用长度为L的区间包含它最优。因为如果区间更长的话,分子不变,分母会变大,不优。这种情况区间长度已经确定,可以不用01分数规划,直接用单调队列\(O(n)\)求即可。类似滑动窗口问题。

2.如果最大值和最小值在序列中的距离>=L,那么它们在区间的左右端点最优。同样如果区间更长的话,分子不变,分母会变大,不优。

按照01分数规划的套路,二分答案mid,判断当前答案是比mid大还是比mid小。如果比mid大,则有

\(|a_i-a_j| \geq mid(j-i+k)\)

(1)$ (a[i]+mid \times i)-(a[j]+mid \times j) \geq mid\times k $,

(区间左端点为最大值)只要暴力枚举i,再在i右侧满足区间长度条件的范围[i+L-1,i+R-1]内求一个\(a[j]+mid \times j\)最小的j即可

(2)$ (a_j-mid \times j)-(a_i-mid \times i) \geq mid \times k ,$

(区间右端点为最大值)只要暴力枚举j,再在j左侧满足区间长度条件的范围[i-L+1,i-R+1]内求一个\(a[j]+mid \times j\)最小的j即可

如果[i,j]内有比a[i]更大的a[k]怎么办?枚举到k的时候就会计入答案,且k为左端点的答案肯定比i优,会覆盖掉i为左端点的不合法情况

用ST表维护最小值即可。

时间复杂度\(O(n \log ^2 n)\),因为每次检查mid都要重新建一次ST表

代码

#include<iostream>
#include<cstdio>
#include<cstring>
#include<algorithm>
#include<queue>
#include<cmath>
#include<deque>
#define INF 1e9
#define eps 1e-6
#define maxn 50000
#define maxlogn 20
using namespace std;
inline void qread(int &x){
x=0;
int sign=1;
char c=getchar();
while(c<'0'||c>'9'){
if(c=='-') sign=-1;
c=getchar();
}
while(c>='0'&&c<='9'){
x=x*10+c-'0';
c=getchar();
}
x=x*sign;
}
int t;
int n,K,L,R;
int a[maxn+5]; //最大值最小值在两端的情况
//最大值最小值在区间两端的时候显然更优
// 01分数规划的套路,二分答案mid
//|a[i]-a[j]| >=mid(j-i+k)
//(1) (a[i]+mid*i)-(a[j]+mid*j)>=mid*k , 只要暴力枚举i,再在i右侧满足区间长度条件的范围内求一个值最小的j
//(2) (a[j]-mid*j)-(a[i]+mid*i)>=mid*k , 只要暴力枚举j,再在j左侧满足区间长度条件的范围内求一个值最小的j
//如果[i,j]内有比a[i]更大的a[k]怎么办?枚举到k的时候就会计入答案,且k为左端点的答案肯定比i优,会覆盖掉i为左端点的不合法情况
int log_2[maxn+5];
double st1[maxlogn+5][maxn+5];//维护a[i]+mid*i的最小值
double st2[maxlogn+5][maxn+5];//维护a[i]-mid*i的最小值
inline double max(double a,double b){
return a>b?a:b;
}
inline double min(double a,double b){
return a<b?a:b;
}
double query1(int l,int r){
int k=log_2[r-l+1];
return min(st1[k][l],st1[k][r-(1<<k)+1]);
}
double query2(int l,int r){
int k=log_2[r-l+1];
return min(st2[k][l],st2[k][r-(1<<k)+1]);
}
bool check(double mid){
for(int i=1;i<=n;i++) st1[0][i]=a[i]+mid*i;
for(int i=1;i<=n;i++) st2[0][i]=a[i]-mid*i;
for(int j=1;(1<<j)<=n;j++){
for(int i=1;i+(1<<j)-1<=n;i++){
st1[j][i]=min(st1[j-1][i],st1[j-1][i+(1<<(j-1))]);
st2[j][i]=min(st2[j-1][i],st2[j-1][i+(1<<(j-1))]);
}
}
double ans=-INF;
for(int i=1;i<=n;i++){
if(i+L-1<=n) ans=max(ans,a[i]+mid*i-query1(i+L-1,min(i+R-1,n)));
}
for(int i=1;i<=n;i++){//这里的i对应的是上面的j
if(i-L+1>=1) ans=max(ans,a[i]-mid*i-query2(max(i-R+1,1),i-L+1));
}
return ans>=mid*K;
} //区间长度为L的时候,区间最大值不一定在左右端点,也可能成为最大答案
//直接长度为L的滑块窗口求解
double sp_solve(){
deque<int>qmax;
deque<int>qmin;
double ans=-INF;
for(int i=1;i<=n;i++){
while(!qmax.empty()&&qmax.front()<i-L+1) qmax.pop_front();
while(!qmin.empty()&&qmin.front()<i-L+1) qmin.pop_front();
while(!qmax.empty()&&a[qmax.back()]<a[i]) qmax.pop_back();
qmax.push_back(i);
while(!qmin.empty()&&a[qmin.back()]>a[i]) qmin.pop_back();
qmin.push_back(i);
ans=max(ans,1.0*(a[qmax.front()]-a[qmin.front()])/(L-1+K));
}
return ans;
}
double bin_search(){
double l=0,r=1000;
double mid;
double res=0;
// printf("db %d\n",check(0+eps)?1:0);
while(r-l>=eps){
mid=(l+r)/2;
if(check(mid)){
res=mid;
l=mid+eps;
}
else r=mid-eps;
}
return res;
}
int main(){
// freopen("1.in","r",stdin);
log_2[0]=-1;
for(int i=1;i<=maxn;i++) log_2[i]=log_2[i>>1]+1;
qread(t);
while(t--){
qread(n);
qread(K);
qread(L);
qread(R);
for(int i=1;i<=n;i++) qread(a[i]);
double ans1=sp_solve();
// printf("ans1=%.4f\n",ans1);
double ans2=bin_search();
double ans=max(ans1,ans2);
printf("%.4f\n",ans);
} }

[BZOJ4476] [JSOI2015] 送礼物 (01分数规划+ST表)的更多相关文章

  1. P6087 [JSOI2015]送礼物 01分数规划+单调队列+ST表

    P6087 [JSOI2015]送礼物 01分数规划+单调队列+ST表 题目背景 \(JYY\) 和 \(CX\) 的结婚纪念日即将到来,\(JYY\) 来到萌萌开的礼品店选购纪念礼物. 萌萌的礼品店 ...

  2. BZOJ4476 JSOI2015送礼物(分数规划+单调队列)

    看到这个式子当然先二分答案.得max-min-(j-i+k)ans>=0. 显然max-min相同的情况下所选区间长度越短越好,所以max和min都应该取在边界.那么实际上我们根本不用管端点是否 ...

  3. BZOJ_4476_[Jsoi2015]送礼物_01分数规划+单调队列

    BZOJ_4476_[Jsoi2015]送礼物_01分数规划+单调队列 Description JYY和CX的结婚纪念日即将到来,JYY来到萌萌开的礼品店选购纪念礼物. 萌萌的礼品店很神奇,所有出售的 ...

  4. [BZOJ4476][JSOI2015]送礼物[分数规划+单调队列]

    题意 题目链接 分析 分数规划之后可以得到式子:\(max-min-r*mid+l*mid\geq k*mid\) . 贪心选择,肯定区间的端点是极小或者极大值.特殊处理区间长度 \(\leq L\) ...

  5. bzoj4476 [Jsoi2015]送礼物

    化简式子 $M>=m+ans*(r-l+k)$ 发现$M,m$确定时,总区间长度越小越好,于是假定右端点为最小值$M+ans*l>=m+ans*r+ans*k$, 右面都确定了,但最大值仍 ...

  6. 【BZOJ4476】[Jsoi2015]送礼物 分数规划+RMQ

    [BZOJ4476][Jsoi2015]送礼物 Description JYY和CX的结婚纪念日即将到来,JYY来到萌萌开的礼品店选购纪念礼物.萌萌的礼品店很神奇,所有出售的礼物都按照特定的顺序都排成 ...

  7. POJ3621Sightseeing Cows[01分数规划 spfa(dfs)负环 ]

    Sightseeing Cows Time Limit: 1000MS   Memory Limit: 65536K Total Submissions: 9703   Accepted: 3299 ...

  8. ZOJ 2676 Network Wars ★(最小割算法介绍 && 01分数规划)

    [题意]给出一个带权无向图,求割集,且割集的平均边权最小. [分析] 先尝试着用更一般的形式重新叙述本问题.设向量w表示边的权值,令向量c=(1, 1, 1, --, 1)表示选边的代价,于是原问题等 ...

  9. P2877 [USACO07JAN]牛校Cow School(01分数规划+决策单调性分治)

    P2877 [USACO07JAN]牛校Cow School 01分数规划是啥(转) 决策单调性分治,可以解决(不限于)一些你知道要用斜率优化却不会写的问题 怎么证明?可以暴力打表 我们用$ask(l ...

随机推荐

  1. 使用pycharm创建git项目的过程

    首先建立远程仓库,然后将远程仓库克隆到本地 然后在pycharm中以该目录创建项目(如果遇到说目录非空,不用管它,Location直接粘贴古来,不然找不到路径) 如果构建好项目说无效的SDK,那么选择 ...

  2. VUE+Ionic,项目搭建&打包成APK

    安装Vue&创建Vue工程 1.安装Vue CLI: npm install -g vue-cli 2.创建新的Vue项目,创建最后一步会提醒是否使用npm install 自动安装,如果选择 ...

  3. win10如何设置软件开机启动

    想要实现应用程序在所有的用户登录系统后都能自动启动,就把该应用程序的快捷方式放到“系统启动文件夹”里C:\ProgramData\Microsoft\Windows\Start Menu\Progra ...

  4. 移动端300ms兼容问题(移动端经典问题)

    移动端300ms延迟原因 2007 年初.苹果公司在发布首款 iPhone 前夕,遇到一个问题:当时的网站都是为大屏幕设备所设计的.于是苹果的工程师们做了一些约定,应对 iPhone 这种小屏幕浏览桌 ...

  5. HTML5基础知识汇总(一)

    一.HTML的开发工具和使用的浏览器 开发工具:记事本等文本编辑器,Atom.VisualStudioCode( VSCode).Brackets.Sublime text和Hbuider. 浏览器: ...

  6. 【bzoj1588】[HNOI2002]营业额统计

    题目描述: 营业额统计 Tiger最近被公司升任为营业部经理,他上任后接受公司交给的第一项任务便是统计并分析公司成立以来的营业情况. Tiger拿出了公司的账本,账本上记录了公司成立以来每天的营业额. ...

  7. Form 表单用 Ajax 提交数据并用 jQuery Validate 验证

    表单填写需要验证可用插件 jQuery Validate 提交数据使用 Ajax 可操控性得到提到 注意:请自行引入 jQuery 和 jQuery Validate HTML 代码 <form ...

  8. 转: Github上关于iOS的各种开源项目集合

    https://blog.csdn.net/jiashaoying/article/details/79079500 下拉刷新 EGOTableViewPullRefresh - 最早的下拉刷新控件. ...

  9. PHP代码修改后提交,无法立即生效

    今天遇到一个坑爹的问题,就是我修改了PHP代码,提交到网站根目录以后,刷新浏览器居然无法立即生效,差不多得一分钟的样子再去刷新,才会看到更改后的效果. 出现这个问题的原因,是因为手头的项目需要较高的P ...

  10. Atcoder arc079 D Decrease (Contestant ver.) (逆推)

    D - Decrease (Contestant ver.) Time limit : 2sec / Memory limit : 256MB Score : 600 points Problem S ...