[BZOJ]4476: [Jsoi2015]送礼物分数规划+RMQ

最新推荐文章于 2020-12-04 05:03:31 发布

200815147

最新推荐文章于 2020-12-04 05:03:31 发布

阅读量314

点赞数 3

分类专栏： 01分数规划

本文链接：https://blog.csdn.net/baidu_36797646/article/details/78692848

版权

01分数规划专栏收录该内容

7 篇文章 0 订阅

订阅专栏

Description

JYY和CX的结婚纪念日即将到来，JYY来到萌萌开的礼品店选购纪念礼物。
萌萌的礼品店很神奇，所有出售的礼物都按照特定的顺序都排成一列，而且相邻的礼物之间有一种神秘的美感。于是，JYY决定从中挑选连续的一些礼物，但究
竟选哪些呢？
【问题描述】
假设礼品店一共有N件礼物排成一列，每件礼物都有它的美观度。排在第i
1< =i< =N个位置的礼物美观度为正整数Ai,。JYY决定选出其中连续的一段，
即编号为礼物i,i+1,…,j-1,j的礼物。选出这些礼物的美观程度定义为
(M（i,j)-m(i,j))/(j-i+k)
其中M(i,j)表示max{Ai,Ai+1….Aj}，m(i,j)表示min{Ai,Ai+1….Aj}，K为给定的正整数。
由于不能显得太小气，所以JYY所选礼物的件数最少为L件；同时，选得太
多也不好拿，因此礼物最多选R件。JYY应该如何选择，才能得到最大的美观程
度？由于礼物实在太多挑花眼，JYY打算把这个问题交给会编程的你。

题解：

唉，连分数规划都不会做了……对于长度为L的序列，直接RMQ求出答案；否则必定是序列的左端点为最小值，右端点为最大值（所以我们要把序列反过来再做一次）。那么假设我们二分出的答案为 $v$ ，那么有

a [ r ] - a [ l ] r - l + 1 > = v

$\frac{a[r]-a[l]}{r-l+1}>=v$ 移项得

(a [r] - r \times v) - (a [l] - l \times v) > = k \times v

$(a[r]-r\times v)-(a[l]-l\times v)>=k\times v$ 然后用单调队列维护一下就好了。

代码：

#include<bits/stdc++.h>
using namespace std;
#define LL long long
#define pa pair<int,int>
const double eps=1e-7;
const int Maxn=50010;
const int inf=2147483647;
int read()
{
    int x=0,f=1;char ch=getchar();
    while(ch<'0'||ch>'9'){if(ch=='-')f=-1;ch=getchar();}
    while(ch>='0'&&ch<='9')x=(x<<3)+(x<<1)+(ch^48),ch=getchar();
    return x*f;
}
int n,k,L,R,a[Maxn],mx[17][Maxn],mn[17][Maxn],Log[Maxn];
double b[Maxn];
int get(int l,int r,int op)//0 max 1 min
{
    int o=Log[r-l+1];
    return((!op)?(max(mx[o][l],mx[o][r-(1<<o)+1])):(min(mn[o][l],mn[o][r-(1<<o)+1])));
}
int q[Maxn];
bool check(double v)
{
    for(int i=1;i<=n;i++)b[i]=(double)(a[i])-(double)(i)*v;
    int l=1,r=0;double o=(double)(k)*v;
    for(int i=1;i<=n;i++)
    {
        while(l<r&&i-q[l]+1>R)l++;
        if(l<=r){if(i-q[l]+1>L&&b[i]-b[q[l]]>=o)return true;}
        while(l<=r&&b[q[r]]>=b[i])r--;
        q[++r]=i;
    }
    for(int i=1;i<=n;i++)b[i]=(double)(a[n-i+1])-(double)(i)*v;
    l=1,r=0;
    for(int i=1;i<=n;i++)
    {
        while(l<r&&i-q[l]+1>R)l++;
        if(l<=r){if(i-q[l]+1>L&&b[i]-b[q[l]]>=o)return true;}
        while(l<=r&&b[q[r]]>=b[i])r--;
        q[++r]=i;
    }
    return false;
}
int main()
{
    Log[1]=0;for(int i=2;i<=50000;i++)Log[i]=Log[i>>1]+1;
    int T=read();
    while(T--)
    {
        n=read(),k=read(),L=read(),R=read();
        for(int i=1;i<=n;i++)a[i]=read(),mx[0][i]=mn[0][i]=a[i];
        for(int i=1;(1<<i)<=n;i++)
        for(int j=1;j+(1<<i)-1<=n;j++)
        mx[i][j]=max(mx[i-1][j],mx[i-1][j+(1<<(i-1))]),
        mn[i][j]=min(mn[i-1][j],mn[i-1][j+(1<<(i-1))]);
        double l=0,r=1000,ans=0;
        for(int i=1;i+L-1<=n;i++)ans=max(ans,(double)(get(i,i+L-1,0)-get(i,i+L-1,1))/(double)(L-1+k));
        while(r-l>=eps)
        {
            double mid=(l+r)*0.5;
            if(check(mid))l=mid+eps;
            else r=mid-eps;
        }
        printf("%.4lf\n",max(l-eps,ans));
    }
}