三分的模版题吧,不过得注意函数是什么,别搞错了。三分的思想就是每次压缩区间吧,所以修改的是最外边的值。一般用来求极值。
代码:
#include<iostream>
#include<cstdio>
#include<vector>
#include<string>
#include<queue>
#include<cstring>
#define maxn 10005
#define INF 0xfffffff
#define mem(a,b) memset(a,b,sizeof(a))
#define FOR(i,s,t) for(int i=s;i<=t;i++)
#define ull unsigned long long
#define ll long long
using namespace std;
inline int rd(int &x)
{
char c=getchar();
while(!isdigit(c))c=getchar();
x=0;
while(isdigit(c))
{
x=x*10+c-'0';
c=getchar();
}
return x;
}
int buf[10];
inline void ot(int i)
{
int p=0;
if(i==0) p++;
else while(i)
{
buf[p++]=i%10;
i/=10;
}
for(int j=p-1; j>=0; j--) putchar('0'+buf[j]);
}
int a[maxn],b[maxn],c[maxn],n;
double f(double x)
{
double ret=-INF;
for(int i=0;i<n;i++)
{
ret=max(ret,a[i]*x*x+b[i]*x+c[i]);
}
return ret;
}
int main()
{
int t;
scanf("%d",&t);
while(t--)
{
scanf("%d",&n);
for(int i=0; i<n; i++)
{
scanf("%d%d%d",&a[i],&b[i],&c[i]);
}
double l=0,r=1000,m1,m2;
for(int ii=0; ii<100; ii++)
{
m1=l+(r-l)/3;
m2=r-(r-l)/3;
if(f(m1)<f(m2))
{
r=m2;
}
else
{
l=m1;
}
}
printf("%.4lf\n",f(l));
}
return 0;
}