codeforces578C. Weakness and Poorness

最新推荐文章于 2020-02-18 12:40:25 发布

原创最新推荐文章于 2020-02-18 12:40:25 发布 · 476 阅读

0 ·

CC 4.0 BY-SA版权

文章标签：

#数学

数学专栏收录该内容

6 篇文章

订阅专栏

本文介绍了一种利用三分法优化复杂子段和问题的解决方案，通过设定初始区间并逐步缩小范围来找到最优解。算法实现包括数学运算、函数调用和循环结构，适用于需要高效求解最大子段和场景。

摘要生成于 C知道，由 DeepSeek-R1 满血版支持，前往体验 >

传送门：http://codeforces.com/problemset/problem/578/c

思路：设f(x)为取x时的最大子段和，f(x)是先减后增的，于是可以用三分法求最值

先确定初始区间[l,r]，mid1=(l+r)/2,mid2=(mid1+r)/2

O(n)求出f(mid1)和f(mid2)

若f(mid1)>f(mid2)则令l=mid1

否则令r=mid2

直到精度达到结束

#include<cmath>
#include<cstdio>
#include<cstring>
#include<algorithm>
typedef long double ld;
const ld eps=1e-12;
const int maxn=200010;
using namespace std;
double a[maxn];ld b[maxn];int n;

ld getans(){
	ld res=0.0,ans=-1.0;
	for (int i=1;i<=n;i++){
		res+=b[i];
		if (res<0) res=0;
		ans=max(ans,res);
	}
	return ans;
}

ld f(ld x){
	ld ans=-1.0;
	for (int i=1;i<=n;i++) b[i]=a[i]-x;
	ans=getans();
	for (int i=1;i<=n;i++) b[i]*=-1.0;
	ans=max(ans,getans());
	return ans;
}

int main(){
	scanf("%d",&n);
	for (int i=1;i<=n;i++) scanf("%lf",&a[i]);
	ld l=-10000.0,r=10000.0,mid1,mid2,res1=-1,res2=-1;
	while (l+eps<=r){
		mid1=(l+r)/2,mid2=(mid1+r)/2;
		res1=f(mid1),res2=f(mid2);
		if (res1>res2+eps) l=mid1;
		else if (res2>res1+eps) r=mid2;
		else break;
	}
	printf("%.15lf\n",(double)res1);
	return 0;
}