记录编号 |
86169 |
评测结果 |
AAAAAAAAAA |
题目名称 |
[UVa 1476] 误差曲线 |
最终得分 |
100 |
用户昵称 |
cstdio |
是否通过 |
通过 |
代码语言 |
C++ |
运行时间 |
1.818 s |
提交时间 |
2014-01-21 20:35:40 |
内存使用 |
0.54 MiB |
显示代码纯文本
#include<iostream>
#include<cstdio>
#include<algorithm>
#include<cstring>
#include<cmath>
#include<vector>
#include<queue>
#include<iomanip>
#include<queue>
#include<set>
#include<map>
using namespace std;
const double eps=1e-12;
const int SIZEN=10001;
int N;
double a[SIZEN]={0},b[SIZEN]={0},c[SIZEN]={0};
double F(double x){
double ans=a[1]*x*x+b[1]*x+c[1];
for(int i=2;i<=N;i++) ans=max(ans,a[i]*x*x+b[i]*x+c[i]);
return ans;
}
double extremum(double (*f)(double),double a,double b,bool type){//求单峰函数f在[a,b]上的极值,type是极值类型
//type=1是下凸函数最小值,type=0是上凸函数最大值
if(b-a<eps) return f((a+b)/2);
double m1,m2;
m1=a+(b-a)/3,m2=b-(b-a)/3;
if(f(m1)>f(m2)){
if(type) return extremum(f,m1,b,type);
else return extremum(f,a,m2,type);
}
else{
if(type) return extremum(f,a,m2,type);
else return extremum(f,m1,b,type);
}
}
void work(void){
scanf("%d",&N);
for(int i=1;i<=N;i++) scanf("%lf%lf%lf",&a[i],&b[i],&c[i]);
printf("%.4lf\n",extremum(F,0,1000,1));
}
int main(){
freopen("errorcurves.in","r",stdin);
freopen("errorcurves.out","w",stdout);
int T;
scanf("%d",&T);
while(T--) work();
return 0;
}