记录编号 236514 评测结果 WWWWWWWWWW
题目名称 [UVa 1476] 误差曲线 最终得分 0
用户昵称 Gravatarliu_runda 是否通过 未通过
代码语言 C++ 运行时间 1.619 s
提交时间 2016-03-14 15:34:51 内存使用 0.36 MiB
显示代码纯文本
#include<cstdio>
#include<cassert>
const int maxn=10005;
const double gold=0.6180339;
int a[maxn],b[maxn],c[maxn],n;
double max(double a,double b){
	return a>b?a:b;
}
double f(double x){
	double p=x*x;
	double ans=p*a[1]+x*b[1]+c[1];
	for(int i=2;i<=n;++i){
		ans=max(ans,p*a[i]+x*b[i]+c[i]);
	}
	return ans;
}
int read(){
	int x;char ch;bool minus=false;
	while(ch=getchar(),(ch<'0'||ch>'9')&&ch!='-');
	if(ch=='-'){
		minus=true;
		ch=getchar();
	}
	x=ch-48;
	while(ch=getchar(),ch<='9'&&ch>='0')x=x*10+ch-48;
	return minus?-x:x;
}

double getmax(){
	double l=0,r=1000,pl=381.96601125,pr=618.03398874989;
	double fl=f(pl),fr=f(pr);
	while(r-l>0.0000001){
		
		if(fl>fr){
			l=pl;pl=pr;fl=fr;
			pr=l+(r-l)*gold;fr=f(pr);
		}else{
			r=pr;pr=pl;fr=fl;
			pl=r-(r-l)*gold;fl=f(pl);
		}
	}
	return f((pl+pr)/2);
}
int main(){
	freopen("errorcurves.in","r",stdin);
	freopen("errorcurves.out","w",stdout);
	int t=read();
	while(t--){
		n=read();
		for(int i=1;i<=n;++i){
			a[i]=read();
			b[i]=read();
			c[i]=read();
		}
		printf("%.1lf\n",getmax());
	}
	fclose(stdin);fclose(stdout);
	return 0;
}