SSL 1015一元三次方程求解

原题在这里

解析

当已知区间(a,b)内有一个根时，用二分法求根，若区间(a,b)内有根，则必有f(a)·f(b)<0。重复执行如下的过程：
(1)若a+0.0001>b或f((a+b)/2)=0，则可确定根为(a+b)/2并退出过程；
(2)若f(a)* f((a+b)/2)<0，则由题目给出的定理可知根在区间(a,(a+b)/2)中，故对区间重复该过程；
(3)若f(a)* f((a+b)/2)>0 ，则必然有f((a+b)/2)* f(b)<0 ，根在((a+b)/2,b)中，对此区间重复该过程。
执行完毕，就可以得到精确到0.0001的根。

code

#include<iostream>
#include<cstdio>
using namespace std;
double a,b,c,d,l,r,mid;
double f(double x){
	return a*x*x*x+b*x*x+c*x+d;
}
int main(){
	scanf("%lf%lf%lf%lf",&a,&b,&c,&d);
	for(int i=-100;i<=100;i++){
		if(f(i)==0){								//假如根为整数直接输出
			printf("%.2lf ",double(i));
		}
		else if(f(i)*f(i+1)<0){
			l=i;r=i+1;								//在小数范围内进行二分查找
			while(1){
				mid=(l+r)/2;
				if(l+0.0001>r or f(mid)==0){
					printf("%.2lf ",mid);break;
				}else if(f(l)*f(mid)<0){
					r=mid;
				}
				else l=mid;
			}
		}
	}
	return 0;
}