BZOJ4476 JSOI2015送礼物（分数规划+单调队列）

　　看到这个式子当然先二分答案。得max-min-(j-i+k)ans>=0。

　　显然max-min相同的情况下所选区间长度越短越好，所以max和min都应该取在边界。那么实际上我们根本不用管端点是否真的是max或min，因为即使不是将他们计入也不会对最终答案造成影响。不妨设右端点是max，则要最大化a_j-a_i-(j-i)ans=(a_j-jans)-(a_i-ians)，单调队列维护即可。左端点是max同理。

　　为了防止不存在长度在l~r的这样的区间，先对长度l的区间单调队列暴力跑一次。

#include<iostream>
#include<cstdio>
#include<cmath>
#include<cstdlib>
#include<cstring>
#include<algorithm>
using namespace std;
int read()
{
    int x=,f=;char c=getchar();
    while (c<''||c>'') {if (c=='-') f=-;c=getchar();}
    while (c>=''&&c<='') x=(x<<)+(x<<)+(c^),c=getchar();
    return x*f;
}
#define N 50010
const double eps=1E-;
int T,n,k,l,r,a[N],p[N],q[N];
bool check(double ans)
{
    int head=,tail=,head1=,tail1=;
    for (int i=;i<=n;i++)
    {
        while (head<=tail&&p[head]<=i-l) head++;
        while (head<=tail&&a[i]>=a[p[tail]]) tail--;
        p[++tail]=i;
        while (head1<=tail1&&q[head1]<=i-l) head1++;
        while (head1<=tail1&&a[i]<=a[q[tail1]]) tail1--;
        q[++tail1]=i;
        if (i>=l&&a[p[head]]-a[q[head1]]>=(l+k-)*ans) return ;
    }
    head=,tail=;
    for (int i=;i<=n;i++)
    {
        while (head<=tail&&q[head]<=i-r) head++;
        if (i>=l)
        {
            while (head<=tail&&a[q[tail]]-q[tail]*ans>=a[i-l+]-(i-l+)*ans) tail--;
            q[++tail]=i-l+;
            if ((a[i]-i*ans)-(a[q[head]]-q[head]*ans)>=k*ans) return ;
        }
    }
    head=,tail=;
    for (int i=;i<=n;i++)
    {
        while (head<=tail&&q[head]<=i-r) head++;
        if (i>=l)
        {
            while (head<=tail&&a[q[tail]]+q[tail]*ans<=a[i-l+]+(i-l+)*ans) tail--;
            q[++tail]=i-l+;
            if ((a[q[head]]+q[head]*ans)-(a[i]+i*ans)>=k*ans) return ;
        }
    }
    return ;
}
int main()
{
#ifndef ONLINE_JUDGE
    freopen("bzoj4476.in","r",stdin);
    freopen("bzoj4476.out","w",stdout);
    const char LL[]="%I64d\n";
#else
    const char LL[]="%lld\n";
#endif
    T=read();
    while (T--)
    {
        n=read(),k=read(),l=read(),r=read();
        for (int i=;i<=n;i++) a[i]=read();
        double L=,R=,ans;
        while (L+eps<=R)
        {
            double mid=(L+R)/;
            if (check(mid)) ans=mid,L=mid+eps;
            else R=mid-eps;
        }
        printf("%.4f\n",ans);
    }
    return ;
}