题解 P6087 【[JSOI2015]送礼物】

分析

先考虑如果没有 \(L\) 和 \(R\) 的长度限制,我们发现最优解中对于区间 \(l\sim r\),\(A_l\) 和 \(A_r\) 应该分别为 \(A_l,A_{l+1},\cdots,A_r\) 中的最大值和最小值。

考虑加上 \(L\) 和 \(R\) 的长度限制,若原来的区间为 \(l\sim r\),则会分如下 3 种情况进行讨论:

  1. 如果 \(r-l+1< L\),那么我们可以考虑把它的长度补成 \(L\)。

  2. 如果 \(L\le r-l+1\le R\),那么这直接可以作为答案。

  3. 如果 \(R<r-l+1\),那么它坏了。

综上所述,只有情况 \(1,2\) 是需要考虑的。

对于第一种情况

区间长度固定为 \(L\),分 \(A_r\) 为最大值和最小值跑一下单调队列就好了。

对于第二中情况

我们再回来看那个式子,发现是个分数规划,非常常规地转换成判断性问题:

\[\begin{aligned} \frac{M(l,r)-m(l,r)}{r-l+K}&\ge mid\\ M(l,r)-m(l,r)&\ge mid\times(r-l+K)\\ M(l,r)-m(l,r)-mid\times r+mid\times l&\ge mid\times K \end{aligned} \]

不妨另 \(A_r\) 为区间 \(l\sim r\) 的最大值,\(A_l\) 为区间 \(l\sim r\) 的最小值。

\[\begin{aligned} A_r-A_l-mid\times r+mid\times l&\ge mid\times K\\ (A_r-mid\times r)-(A_l-mid\times l)&\ge mid\times K \end{aligned} \]

用单调队列实现一下就可以了。

\(A_r\) 为区间 \(l\sim r\) 的最小值,\(A_l\) 为区间 \(l\sim r\) 的最大值同理。

实现细节

就是极为有用的小 trick——弱化最大值,最小值。

你会发现用上述方法实现非常困难,因为你要强制取最大或取最小。

但是稍作思考,你会发现,如果不取最大值最小值答案一定不会更优,所以写代码的时候完全不需要考虑。

#include<bits/stdc++.h>
#define log(a) cerr<<"\033[32m[DEBUG] "<<#a<<'='<<(a)<<" @ line "<<__LINE__<<"\033[0m"<<endl
#define LL long long
#define SZ(x) ((int)x.size()-1)
#define ms(a,b) memset(a,b,sizeof a)
#define F(i,a,b) for(int i=(a);i<=(b);++i)
#define DF(i,a,b) for(int i=(a);i>=(b);--i)
using namespace std;
inline int read(){char ch=getchar(); int w=1,c=0;
	for(;!isdigit(ch);ch=getchar()) if (ch=='-') w=-1;
	for(;isdigit(ch);ch=getchar()) c=(c<<1)+(c<<3)+(ch^48);
	return w*c;
}
const int N=5e4+10;
const double eps=1e-6;
int a[N],n,k,L,R,q[N],s,t;
double b[N],c[N];
inline double situation1(){
	double ans=0;
	s=1,t=0;
	F(i,1,n){
		while(s<=t&&i-q[s]>=L)s++;
		if(s<=t)ans=max(ans,(a[i]-a[q[s]])/((double)(L+k-1)));
		while(s<=t&&a[q[t]]>a[i])t--;
		q[++t]=i;
	}
	s=1,t=0;
	F(i,1,n){
		while(s<=t&&i-q[s]>=L)s++;
		if(s<=t)ans=max(ans,(a[q[s]]-a[i])/((double)(L+k-1)));
		while(s<=t&&a[q[t]]<a[i])t--;
		q[++t]=i;
	}return ans;
}
inline bool check(double x){
	F(i,1,n)b[i]=a[i]-x*i,c[i]=a[i]+x*i;
	s=1,t=0;
	F(i,1,n){
		while(s<=t&&i-q[s]+1>R)s++;
		if(s<=t&&b[i]-b[q[s]]>=x*k)return true;
		if(i>=L){
			while(s<=t&&b[q[t]]>b[i-L+1])t--;
			q[++t]=i-L+1;
		}
	}
	s=1,t=0;
	F(i,1,n){
		while(s<=t&&i-q[s]+1>R)s++;
		if(s<=t&&c[q[s]]-c[i]>=x*k)return true;
		if(i>=L){
			while(s<=t&&c[q[t]]<c[i-L+1])t--;
			q[++t]=i-L+1;
		}
	}
	return false;
}
inline double situation2(){
	double l=0,r=10001;
	while(l+eps<r){
		double mid=(l+r)/2.0;
		if(check(mid))l=mid;
		else r=mid;
	}return l;
}
signed main(){
	int _=read();
	while(_--){
		n=read(),k=read(),L=read(),R=read();
		F(i,1,n)a[i]=read();
		cout<<fixed<<setprecision(4)<<max(situation1(),situation2())<<endl;
	}
	return 0;
}
上一篇:「JSOI2015」串分割


下一篇:mvc源码解析2-1:mvchandler的执行