输入
3 4 2 2 3 1 2 4 1 3 2
输出
4.6666666667
题目大意:给你两串序列a1~an,b1~bm,根据这两串序列可以构成一个n*m规模的矩阵,问题是让我们在这个矩阵中找到一个规模大于等于x*y规模的,并且矩阵中所有元素的平均值最大的矩阵,并输出这个最大平均值。
题目解析:因为题目中已经告诉我们Wi,j=ai+bj,所以如果想要求出平均值最大的矩阵,让这个矩阵中对应的a序列和b序列的平均值分别最大即可,这样我们就把一个二维问题转化成了一个一维问题,所以现在的问题就是在两个序列中找到平均值最大的一段的平均值是多少,因为我们发现序列中最大值是1e5所以,我们可以采用二分平均值,然后去序列中判断是否存在大于等于指定长度的一段即可,最后矩阵的平均值就是两段序列最大平均值的加和,做题的过程中二分平均值是我没想到的,而且代码实现起来也不是太容易,有许多细节需要考虑,通过这道题也充分认识到了自己的不足之处。
上代码:
#include <iostream>
#include <cstring>
#include <cstdio>
#include <algorithm>
#include <cmath>
const int N=1e5+10;
using namespace std;
int a[N],b[N];
double temp[N];
int n,m,x,y;
bool checka(double t)
{
for(int i=1;i<=n;i++)
temp[i]=temp[i-1]+a[i]-t;
double p=0x3f3f3f3f;
for(int i=x;i<=n;i++)
{
p=fmin(p,temp[i-x]);
if(temp[i]-p>=0)
return true;
}
return false;
}
bool checkb(double t)
{
for(int i=1;i<=m;i++)
temp[i]=temp[i-1]+b[i]-t;
double p=0x3f3f3f3f;
for(int i=y;i<=m;i++)
{
p=fmin(p,temp[i-y]);
if(temp[i]-p>=0)
return true;
}
return false;
}
int main()
{
cin>>n>>m>>x>>y;
for(int i=1;i<=n;i++)
cin>>a[i];
for(int i=1;i<=m;i++)
cin>>b[i];
double l=0,r=100000;
double ans=0;
while(r-l>=1e-10)
{
double mid=(l+r)/2;
if(checka(mid))
l=mid;
else
r=mid;
}
ans+=(l+r)/2;
l=0;
r=100000;
while(r-l>=1e-10)
{
double mid=(l+r)/2;
if(checkb(mid))
l=mid;
else
r=mid;
}
printf("%.10lf\n",ans+(l+r)/2);
return 0;
}