题意:
每个数都减x,找到合适的x,使得最大的abs(连续子序列的和)->最小
解法:
三分,类开口朝上的二次函数
实例代码:
#include <cstdio>
#include <iostream>
#include <cmath>
#include <cstring>
#include <algorithm>
#include <stdlib.h>
#include <map>
#include <vector>
#define rep(i,n) for(i=0;i<n;i++)
#define cle(x) memset(x,0,sizeof(x))
#define ll long long
const int maxn=200000+5;
using namespace std;
ll a[maxn];
long double b[maxn];
long double ans;
int n;
long double doit (long double x){
for(int i=1;i<=n;i++)b[i]=a[i];
for(int i=1;i<=n;i++)b[i]-=(x*i);
long double mi=0;
long double ma=0;
long double re=0;
for(int i=1;i<=n;i++){
re=max(re,max(abs(b[i]-mi),abs(b[i]-ma)));
mi=min(mi,b[i]);
ma=max(ma,b[i]);
}
ans=min(ans,re);
return re;
}
int main()
{
#ifndef ONLINE_JUDGE
freopen("in.txt","r",stdin);
//freopen("out.txt","w",stdout);
#endif
while(cin>>n){
int i;
for(i=1;i<=n;i++)scanf("%I64d",&a[i]);
for(i=1;i<=n;i++)a[i]=a[i-1]+a[i];
long double l=-10000,r=10000;
int T=100;
ans=3000000000.0;
while(T--){
long double ml=l+(r-l)/3.0;
long double mr=r-(r-l)/3.0;
if(doit(ml)>doit(mr))l=ml;
else r=mr;
}
printf("%.12lf\n",(double)ans);
}
return 0;
}