洛谷1024 一元三次方程求解

原题地址

https://www.luogu.org/problem/show?pid=1024

二分答案

解题思路

据说暴力可过，只能说数据水了怎么搞都能过……反正这题可以拿来练二分，只是不完全符合二分性质，每次判断中点是否可行之后左右区间都需要再继续搜，直到三个解都被找到或者区间长度<=0.001（必须是10^-3不是10^-2）

参考代码

#include
#include
#include
#include
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;
}