【Openjudge】一元三次方程求解_二分

7891:一元三次方程求解
總時間限制:  1000ms  記憶體限制:  65536kB 描述

有形如：ax3+bx2+cx+d=0  這樣的一個一元三次方程。

給出該方程中各項的係數(a，b，c，d  均為實數)，並約定該方程存在三個不同實根(根的範圍在-100至100之間)，且根與根之差的絕對值>=1。要求由小到大依次在同一行輸出這三個實根(根與根之間留有空格)，並精確到小數點後2位。 輸入 一行，包含四個實數a，b，c，d，相鄰兩個數之間用單個空格隔開。 輸出 一行，包含三個實數，為該方程的三個實根，按從小到大順序排列，相鄰兩個數之間用單個空格隔開，精確到小數點後2位。 範例輸入

1.0 -5.0 -4.0 20.0

範例輸出

-2.00 2.00 5.00

根據題目描述，求根的過程即為二分，只要f(l)*f(r)<=0,其間必有根。 由於一定存在三個根，且兩根的距離大於1，故從左往右，依次掃描、二分。
注意：因為此題要用到浮點，注意函數參數為Double型。 我就這麼Wa了很久 。

#include<iostream>#include<cstdio>#include<cstring>#include<cmath>#include<cstdlib>#include<queue>#include<stack>#include<algorithm>using namespace std;#define MAXN#define MAXM#define INF 0x3f3f3f3f#define LL long longdouble a,b,c,d;double f(double x){return a*x*x*x+b*x*x+c*x+d;}double dfs(double l,double r){if(r-l<=0.001)return l;double mid=(l+r)/2;if(f(l)*f(mid)<=0)return dfs(l,mid);else return dfs(mid,r);}int main(){scanf("%lf%lf%lf%lf",&a,&b,&c,&d);double ans[5];int cnt=0;for(int i=-100;i<100;++i){double l=i,r=i+1;if(f(l)==0)ans[++cnt]=l;else if(f(l)*f(r)<0)ans[++cnt]=dfs(l,r);if(cnt>=3)break;}printf("%0.2lf %0.2lf %0.2lf\n",ans[1],ans[2],ans[3]);}