记录编号 |
236514 |
评测结果 |
WWWWWWWWWW |
题目名称 |
[UVa 1476] 误差曲线 |
最终得分 |
0 |
用户昵称 |
liu_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;
}