bzoj3316: JC loves Mkk（单调队列+分数规划）

Description

Input

第1行，包含三个整数。n，L，R。
第2行n个数，代表a[1..n]。

Output

仅1行，表示询问答案。
如果答案是整数，就输出整数；否则，输出既约分数“P/Q”来表示。

 

Sample Input

5 3 4
3 1 2 4 5

Sample Output

7/2

HINT
1≤L≤R≤n≤10^5，0≤ai≤10^9，保证问题有解,数据随机生成

 

 

首先这是一个分数规划，于是我们得二分，设答案为mid，那么原数列变成a[i]-mid，然后就是要找一段使得区间和大于0

前缀和可以先预处理，然后找到满足s[j]<s[i]且i<j的j，发现满足条件的j中s[j]越小越好，于是用单调队列维护

然后得保证选的数的个数是偶数，于是开两个单调队列，分别维护位置为奇数和偶数的

//minamoto
#include<iostream>
#include<cstdio>
#define ll long long
using namespace std;
#define getc() (p1==p2&&(p2=(p1=buf)+fread(buf,1,1<<21,stdin),p1==p2)?EOF:*p1++)
char buf[1<<21],*p1=buf,*p2=buf;
template<class T>inline bool cmax(T&a,const T&b){return a<b?a=b,1:0;}
template<class T>inline bool cmin(T&a,const T&b){return a>b?a=b,1:0;}
int read(){
    #define num ch-'0'
    char ch;bool flag=0;int res;
    while(!isdigit(ch=getc()))
    (ch=='-')&&(flag=true);
    for(res=num;isdigit(ch=getc());res=res*10+num);
    (flag)&&(res=-res);
    #undef num
    return res;
}
const int N=1e5+5;
int n,m,L,R,h1,h2,t1,t2;ll ans1,ans2,g,A[N<<1],S[N<<1];
double v[N<<1],s[N<<1];int q1[N<<1],q2[N<<1];
ll gcd(ll a,ll b){return b?gcd(b,a%b):a;}
bool check(double x){
    for(int i=1;i<=m;++i) v[i]=A[i]-x,s[i]=s[i-1]+v[i];
    h1=h2=t1=1,t2=0,q1[1]=0;
    for(int i=L;i<=m;++i){
        while(h1<=t1&&q1[h1]<i-R) ++h1;
        while(h2<=t2&&q2[h2]<i-R) ++h2;
        if(!(i&1)&&h1<=t1&&s[q1[h1]]<=s[i]){
            ans1=S[i]-S[q1[h1]],ans2=i-q1[h1],g=gcd(ans1,ans2),ans1/=g,ans2/=g;return 1;
        }
        if((i&1)&&h2<=t2&&s[q2[h2]]<=s[i]){
            ans1=S[i]-S[q2[h2]],ans2=i-q2[h2],g=gcd(ans1,ans2),ans1/=g,ans2/=g;return 1;
        }
        if(!((i-L+1)&1)){
            while(h1<=t1&&s[q1[t1]]>=s[i-L+1]) --t1;
            q1[++t1]=i-L+1;
        }else{
            while(h2<=t2&&s[q2[t2]]>=s[i-L+1]) --t2;
            q2[++t2]=i-L+1;
        }
    }
    return 0;
}
int main(){
//    freopen("testdata.in","r",stdin);
    n=read(),L=read(),R=read(),m=n<<1;
    double l=1<<30,r=0;
    for(int i=1;i<=n;++i) A[i]=A[i+n]=read(),cmin(l,(double)A[i]),cmax(r,(double)A[i]);
    for(int i=1;i<=m;i++) S[i]=S[i-1]+A[i];
    for(int i=1;i<=50;++i){
        double mid=(l+r)/2;
        check(mid)?l=mid:r=mid;
    }
    printf("%lld/%lld",ans1,ans2);
    return 0;
}