1. 程式人生 > 實用技巧 >洛谷 P1024 一元三次方程求解(二分)

洛谷 P1024 一元三次方程求解(二分)

題目連結:https://www.luogu.com.cn/problem/P1024

首先在(-100,100)上暴力列舉,如果有f(i)*f(j)<=0,那麼答案一定在[i,j]上,在這裡記錄下f在[i,j]上的單調性,進行二分。二分結束後繼續在(ans,100)上暴力列舉,重複是上述操作。

AC程式碼:

 1 #include<cstdio>
 2 #include<iostream>
 3 using namespace std;
 4 double a,b,c,d;
 5 double f(double x){
 6     return a*x*x*x+b*x*x+c*x+d;
7 } 8 double erfen(double l,double r,int flag){ 9 while(l+0.0001<r){ 10 double mid=(l+r)/2; 11 if(f(mid)>=0){ 12 if(flag) r=mid; 13 else l=mid; 14 } 15 else{ 16 if(flag) l=mid; 17 else r=mid; 18 } 19 }
20 return l; 21 } 22 int main(){ 23 scanf("%lf%lf%lf%lf",&a,&b,&c,&d); 24 double low=-99.99; 25 for(double i=low;i<100;i++){ 26 if(f(i)*f(low)<=0){ 27 int flag=(f(i)>f(low)); 28 double ans=erfen(low,i,flag); 29 printf("
%.2f ",ans); 30 low=i; 31 } 32 } 33 return 0; 34 }
AC程式碼