三分法的模板,一直提交,一直WA
结果,今天把精度判断从l-r>=t
改成 fabs(F(l)-F(r))>=t
,过了三个点,再将 \(t\) 改成 \(10^{-7}\) 就 \(AC\) 了。。。
这给我一个十分重要的启示:精度问题十分重要!!!
#include <iostream>
#include <cstdio>
#include <cmath>
using namespace std;
const int MAXN = 10000+5;
int a[MAXN],b[MAXN],c[MAXN],n;
inline 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;
}
int main() {
int T;
cin >> T;
while(T--) {
cin >> n;
for(int i=1; i<=n; i++) cin >> a[i] >> b[i] >> c[i];
double l=0,r=1000,t=1e-7;
while(fabs(F(l)-F(r))>t) {
double p=(2*l+r)/3,q=(l+2*r)/3;
if(F(p) > F(q)) {
l=p;
} else r=q;
}
printf("%.4lf\n",F(l));
}
return 0;
}
原文地址:https://www.cnblogs.com/zhinv/p/9816676.html
时间: 2024-11-08 22:35:30