P1024 [NOIP2001 提高组] 一元三次方程求解

题目描述

有形如：a x^3 + b x^2 + c x + d = 0 这样的一个一元三次方程。给出该方程中各项的系数（a,b,c,d均为实数），并约定该方程存在三个不同实根（根的范围在 -100至100 之间），且根与根之差的绝对值 ≥1。要求由小到大依次在同一行输出这三个实根(根与根之间留有空格)，并精确到小数点后 2 位。

提示：记方程 f(x) = 0f(x)=0，若存在 2 个数 x1​ 和 x2​，且 x1​<x2​,f(x1​)×f(x2​)<0，则在(x1​,x2​) 之间一定有一个根。

输入格式

一行，44 个实数 a, b, c, d。

输出格式

一行，3 个实根，从小到大输出，并精确到小数点后 2位。

输入输出样例

输入 #1复制

1 -5 -4 20

输出 #1复制

-2.00 2.00 5.00

说明/提示

【题目来源】

NOIP 2001 提高组第一题

 解题方法：二分。

Code:

#include <iostream>
#include <stdio.h>

using namespace std;

const int N = 5000;

double a, b, c, d, l, r, mid, x1, x2;

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

signed main()
{
	cin >> a >> b >> c >> d;
	for(int i = -100; i < 100; i ++ )
	{
		l = i, r = i + 1;
		x1 = g(l), x2 = g(r);
		if(x1 == 0)  printf("%.2lf ", l);
		else
		{
			if(x1 * x2 < 0)
			{
				while(r - l > 1e-8)
				{
					mid = (l + r) / 2;
					if(g(mid) * g(r) >= 0)  r = mid;
					else  l = mid;
				}
				printf("%.2lf ", l);
			}
		}
	}
}