洛谷OJ-P 1024-一元三次方程式の解(二分解答)


タイトルの説明
有形は、ax 3+bx 2+cx+d=0のような一元三次方程式である.この方程式の各係数(a,b,c,dはいずれも実数)を与え、この方程式には3つの異なる実根(根の範囲は-100から100の間)が存在し、根と根の差の絶対値>=1が規定されている.この3つの実ルート(ルートとルートの間にスペースが残っている)を小さいから大きいまで順次同じ行に出力し、小数点以下2桁まで正確にする必要があります.ヒント:方程式f(x)=0を記し、2つの数x 1とx 2が存在し、x 1
入力
1行、4個の実数A,B,C,D.
しゅつりょく
1行、3つの実根を小数点以下2桁に正確にします.
サンプル入力
1 -5 -4 20

サンプル出力
-2.00 2.00 5.00

テーマの考え方
問題はルートが-100から100の間にあり、ルートとルートの差の絶対値>=1である場合、ルートがこの2つの数字の間にあると判断された場合、各2つの隣接する数字を列挙する.では、この2つの数字の中で2点で答えを探せばいいです.
タイトルコード
#include  
#include 
#include 
#include 
#include 
#include 
#include 
#include 
#include 
#include 
#define LL long long  
using namespace std;
double ans1, ans2, ans3;
double a, b, c, d;
double l, r, mid;


double f1(double x){
	return x*x*x*a + x*x*b + x*c +d;	
}

double F(double x, double y){
	l = x; r = y;
	while(l < r){
		mid = (l+r) / 2.0;
		double ans1 = f1(l);
		double ans2 = f1(mid);
		if((ans1*ans2) <= 0)
			r = mid;		
		else
			l = mid + 0.00001;
	}
	return r;
}

int main(){

	while(scanf("%lf%lf%lf%lf",&a,&b,&c,&d) != EOF){
		for(double i = -100; i <= 100; i++){
			double ans1 = f1(i);
			double ans2 = f1(i+1);
			if(ans1 == 0)
				printf("%.2f ",i);
			
			if((ans1*ans2) < 0){
				ans1 = F(i, i+1.0);
				printf("%.2f ",ans1);
			}
		}
		printf("
"); } return 0; }