牛客多校赛4

牛客多校赛4

Average

一条典中典的简单的题目
稍微观察一下可以发现这条题目可以转化为序列 a i a_i ai b i b_i bi的一段区间 [ l , r ] [l,r] [l,r]的最大均值。可以二分答案,再找是否有一个区间 [ l , r ] [l,r] [l,r]的均值比二分时的 m i d mid mid大。维护前缀和 − - 区间长度* m i d mid mid。先扫一遍确定 r r r,维护 [ 0 , l ] [0,l] [0,l]中的最小值,只要用一个值记录一下就好,一开始想了主席树有想到堆,太过离谱了

#include<bits/stdc++.h>
using namespace std;
const int N = 1e5+10;
double a[N],b[N],sum_a[N],sum_b[N];
int n,m,x,y;
bool calca(double mid){
	sum_a[0] = 0;
	for(int i = 1;i<=n;i++){
		sum_a[i] = sum_a[i-1]+a[i]-mid;
		if((i>=x)&&(sum_a[i]>=0))return true;
	}
	double mn =  sum_a[0];
	for(int i = x;i<=n;i++){
		mn = min(mn,sum_a[i-x]);
		if(sum_a[i]>=mn)return true;
	}
	return false;
}
bool calcb(double mid){
	sum_b[0] = 0;
	for(int i = 1;i<=m;i++){
		sum_b[i] = sum_b[i-1]+b[i]-mid;
		if((i>=y)&&(sum_b[i]>=0))return true;
	}
	double mn = sum_b[0];
	for(int i = y;i<=m;i++){
		mn = min(mn,sum_b[i-y]);
		if(sum_b[i]>=mn)return true;
	}
	return false;
}
int main(){
	scanf("%d%d%d%d",&n,&m,&x,&y);
	double l1=N,l2=N,r1=-1,r2=-1;
	for(int i = 1;i<=n;i++){
		scanf("%lf",&a[i]);
		l1 = min(l1,a[i]);
		r1 = max(r1,a[i]);
	}
	for(int i = 1;i<=m;i++){
		scanf("%lf",&b[i]);
		l2 = min(l2,b[i]);
		r2 = max(r2,b[i]);
	}
	double mid1,mid2;
	for(int i=0;i<100;i++){
    	mid1=(l1+r1)/2;
    	if(calca(mid1))l1=mid1;
    	else r1=mid1;
	}
	for(int i=0;i<100;i++){
    	mid2=(l2+r2)/2;
    	if(calcb(mid2))l2=mid2;
    	else r2=mid2;
	}
	printf("%.10lf",mid1+mid2);
	return 0;
}

未完待续

评论
添加红包

请填写红包祝福语或标题

红包个数最小为10个

红包金额最低5元

当前余额3.43前往充值 >
需支付:10.00
成就一亿技术人!
领取后你会自动成为博主和红包主的粉丝 规则
hope_wisdom
发出的红包
实付
使用余额支付
点击重新获取
扫码支付
钱包余额 0

抵扣说明:

1.余额是钱包充值的虚拟货币,按照1:1的比例进行支付金额的抵扣。
2.余额无法直接购买下载,可以购买VIP、付费专栏及课程。

余额充值