算法竞赛宝典 分治算法 一元三次方程

一元三次方程求解

总时间限制: 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

//利用精度逐个枚举

#include<bits/stdc++.h>
using namespace std;
double a,b,c,d;
double f(double x)
{
	return ((a*x+b)*x+c)*x+d; 
}
int main()
{
	cin>>a>>b>>c>>d;
	for(int i=-10000; i<=10000; i++)
	{
		double x=(double)i/100;
		if(f(x)>=-0.00001 && f(x)<=0.00001)
		{
			printf("%.2lf ",x);
		}
	}
	return 0;
}

//二分算法首当其冲，也是利用精度，还用利用题目的根与根之差的绝对值>=1 的条件对根的f(x1)*f(x2)进行优化处理

//对二分算法有了进一步的认知，循环推出的条件又多种多样，根据题目的意思去做出判断，特别是对精度问题的处理

这个二分算法是对穷举算法的优化
使得一些没有可能有答案的区间进行精度确认，进而直接跳过
减少时间的消耗

#include<algorithm>
#include<iostream>
using namespace std;
int n;
double ans[3];
double a, b, c, d;

double f(double x)
{
	return ((a*x + b)*x + c)*x + d;
}
int main()
{
	while (cin >> a >> b >> c >> d)
	{
		double x1, x2,xx;
		n = 0;
		for (int i = -100; i <= 100; i++)
		{
			x1 = i;
			x2 = i + 1;                //因为二分的时候x1，x2会变，所以用i固定
			if (f(x1) == 0)
				ans[n++] = x1;
			else if (f(x1)*f(x2) < 0)
			{
				while (x2 - x1 >= 0.001)        //注意这一的终止边界条件，
				{
					xx = (x1 + x2) / 2;
					if ((f(x1)*f(xx)) <= 0)            //二分的细节处理
						x2 = xx;
					else
						x1 = xx;
				}
				ans[n++] = x1;
			}
		}
		printf("%.2lf %.2lf %.2lf\n", ans[0], ans[1], ans[2]);
	}
	return 0;
}