BZOJ3316: JC loves Mkk

最新推荐文章于 2020-01-13 20:23:43 发布

forezxl

最新推荐文章于 2020-01-13 20:23:43 发布

阅读量199

点赞数 1

分类专栏： BZOJ 数据结构---单调队列其他---二分/三分蒟蒻zxl的Blog专栏文章标签： BZOJ 二分单调队列

本文链接：https://blog.csdn.net/a1799342217/article/details/82749472

版权

蒟蒻zxl的Blog专栏同时被 3 个专栏收录

460 篇文章 6 订阅

订阅专栏

BZOJ

281 篇文章 1 订阅

订阅专栏

其他---二分/三分

14 篇文章 0 订阅

订阅专栏

二分单调队列

题目传送门

套路断环为链。二分答案mid(用long double)，记一个前缀和 $s [i] = s [i - 1] + w [i] - m i d$ ，那么我们只要判断是否有一段长度为偶数且 $s[r]-s[l-1]\geq0$ 即可。

用单调队列维护 $s [i]$ 的最小值。保证长度为偶数的话用两个分奇偶维护。

二分完后我们可以得到分母，乘上答案四舍五入就可以得到分子。

这道题充满了套路。。。

代码：

#include<cctype>
#include<cstdio>
#include<cstring>
#include<algorithm>
#define N 200005
#define F inline
#define eps 1e-6
using namespace std;
typedef long long LL;
typedef long double DB;
int n,L,R,mx,w[N],q[2][N];
DB s[N]; LL ans1,ans2;
F char readc(){
	static char buf[100000],*l=buf,*r=buf;
	if (l==r) r=(l=buf)+fread(buf,1,100000,stdin);
	return l==r?EOF:*l++;
}
F int _read(){
	int x=0; char ch=readc();
	while (!isdigit(ch)) ch=readc();
	while (isdigit(ch)) x=(x<<3)+(x<<1)+(ch^48),ch=readc();
	return x;
}
F LL gcd(LL a,LL b){ return !b?a:gcd(b,a%b); }
bool pd(DB mid){
	for (int i=1;i<=n;i++) s[i]=s[i-1]+w[i]-mid;
	int l[2]={1,1},r[2]={0,0};
	for (int i=L;i<=n;i++){
		int p=i-L,x=i&1;
		while (l[x]<=r[x]&&s[p]<s[q[x][r[x]]]) r[x]--;
		q[x][++r[x]]=p; if (i-q[x][l[x]]>R) l[x]++;
		if (s[i]-s[q[x][l[x]]]>0) return ans2=i-q[x][l[x]],true;
	}
	return false;
}
int main(){
	n=_read(),L=_read(),R=_read(),L+=L&1,R-=R&1;
	for (int i=1;i<=n;i++) w[i]=w[i+n]=_read(),mx=max(mx,w[i]);
	DB l=0,r=mx,mid; n<<=1;
	while (r-l>eps) pd(mid=(l+r)/2)?l=mid:r=mid;
	ans1=(LL)((l+r)/2*ans2+0.5);
	LL t=gcd(ans1,ans2); ans1/=t,ans2/=t;
	return ans2==1?printf("%lld",ans1):printf("%lld/%lld",ans1,ans2),0;
}