2021牛客暑期多校训练营4 - J - Average（ 二分 + 最大连续字段和 ）

题目链接：点击进入

题目

题意

求矩阵的（ 宽 >= x , 高 >= y 的平均值最大的 ）子矩阵。

思路

要是看懂了上面的图片，就明白问题就转换成了分别求 a、b 数组的 有长度限制的（ a 是不小于 x , b 是不小于 y ） 连续子序列的最大平均值。这个问题不懂的可以点点这里。

代码

//#pragma GCC optimize(3)//O3
//#pragma GCC optimize(2)//O2
#include<iostream>
#include<string>
#include<map>
#include<set>
//#include<unordered_map>
#include<queue>
#include<cstdio>
#include<vector>
#include<cstring>
#include<stack>
#include<algorithm>
#include<iomanip>
#include<cmath>
#include<fstream>
#define X first
#define Y second
#define best 131
#define INF 0x3f3f3f3f3f3f3f3f
#define pii pair<int,int>
#define lowbit(x) x & -x
#define inf 0x3f3f3f3f
//#define int long long
//#define double long double
//#define rep(i,x,y) for(register int i = x; i <= y;++i)
using namespace std;
typedef long long ll;
typedef unsigned long long ull;
const double pai=acos(-1.0);
const int maxn=1e6+10;
const int mod=998244353;
const double eps=1e-6;
const int N=5e3+10;
 
//inline int read()
//{
//    int k = 0, f = 1 ;
//    char c = getchar() ;
//    while(!isdigit(c)){if(c == '-') f = -1 ;c = getchar() ;}
//    while(isdigit(c)) k = (k << 1) + (k << 3) + c - 48 ,c = getchar() ;
//    return k * f ;
//}
 
int n,m,x,y,a[maxn],b[maxn];
double sum[maxn];
 
bool check1(double aver)
{
    for(int i=1;i<=n;i++)
        sum[i]=sum[i-1]+a[i]-aver;
    double minn=1e10;
    double ans=-1e10;
    for(int i=x;i<=n;i++)
    {
        minn=min(minn,sum[i-x]);
        ans=max(ans,sum[i]-minn);
    }
    return ans>=0;
}
 
bool check2(double aver)
{
    for(int i=1;i<=m;i++)
        sum[i]=sum[i-1]+b[i]-aver;
    double minn=1e10;
    double ans=-1e10;
    for(int i=y;i<=m;i++)
    {
        minn=min(minn,sum[i-y]);
        ans=max(ans,sum[i]-minn);
    }
    return ans>=0;
}
 
int main()
{
//  ios::sync_with_stdio(false);
//  cin.tie(0);cout.tie(0);
    scanf("%d%d%d%d",&n,&m,&x,&y);
    for(int i=1;i<=n;i++) scanf("%d",&a[i]);
    for(int i=1;i<=m;i++) scanf("%d",&b[i]);
    double l=0,r=1e10,ans1=0,ans2=0;
    while((r-l)>=eps)
    {
        double mid=(l+r)/2.0;
        if(check1(mid))
            l=mid,ans1=mid;
        else
            r=mid;
    }
//  cout<<ans<<endl;
    l=0,r=1e10;
    while((r-l)>=eps)
    {
        double mid=(l+r)/2.0;
        if(check2(mid))
            l=mid,ans2=mid;
        else
            r=mid;
    }
    printf("%.10lf",ans1+ans2);
    return 0;
}