原题地址
https://www.luogu.org/problem/show?pid=1024
二分答案
解题思路
据说暴力可过,只能说数据水了怎么搞都能过……反正这题可以拿来练二分,只是不完全符合二分性质,每次判断中点是否可行之后左右区间都需要再继续搜,直到三个解都被找到或者区间长度<=0.001(必须是10^-3不是10^-2)
参考代码
#include<iostream>
#include<algorithm>
#include<cstdio>
#include<cmath>
using namespace std;
bool ff=0;
int cnt=0;
long double a,b,c,d,f1,f2,f0;
long double p[5];
long double f(long double x)
{
long double t0,t1,t2,t3;
t3=a*x*x*x;
t2=b*x*x;
t1=c*x;
t0=d;
return t0+t1+t2+t3;
}
void ef(long double l,long double r,int k)
{
if(r-l<1&&k==0||ff) return;
if(r-l<=0.001&&k==1)
{
cnt++;
p[cnt]=l;
if (cnt==3) ff=1;
return;
}
if(r-l<=0.001&&k==0) return;
long double mid=(l+r)/2;
f1=f(l);
f2=f(r);
f0=f(mid);
if (f0==0)
{
cnt++;
p[cnt]=mid;
if (cnt==3) ff=1;
if (f1>0)ef(l,mid,0);
else ef(l,mid,1);
if (f2>0)ef(mid,r,0);
else ef(mid,r,1);
}
if (f1*f0<0)ef(l,mid,1);
else ef(l,mid,0);
if (f2*f0<0)ef(mid,r,1);
else ef(mid,r,0);
}
int main()
{
a,b,c,d;
cin>>a>>b>>c>>d;
long doublel=-100.0,r=100.0;
ef(l,r,1);
sort(p+1,p+4);
for (int i=1;i<=3;i++)
printf("%.2llf%c",p[i],' ');
return 0;
}