2018年信息学奥赛NOIP资料下载
题目描述
有形如:ax3+bx2+cx1+dx0=0ax
3
+bx
2
+cx
1
+dx
0
=0 这样的一个一元三次方程。给出该方程中各项的系数(a,b,c,da,b,c,d均为实数),并约定该方程存在三个不同实根(根的范围在-100−100至100100之间),且根与根之差的绝对值 \ge 1≥1。要求由小到大依次在同一行输出这三个实根(根与根之间留有空格),并精确到小数点后22位。
提示:记方程f(x)=0f(x)=0,若存在22个数x_1x
1
和x_2x
2
,且x_1<x_2x
1
<x
2
,f(x_1) \times f(x_2)<0f(x
1
)×f(x
2
)<0,则在(x_1,x_2)(x
1
,x
2
)之间一定有一个根。
输入输出格式
输入格式:
一行,44个实数A,B,C,DA,B,C,D。
输出格式:
一行,33个实根,并精确到小数点后22位。
输入输出样例
输入样例#1:
1 -5 -4 20
输出样例#1:
-2.00 2.00 5.0
#include<iostream>
#include<cstdio>
#include<cmath>
#include<algorithm>
using namespace std;
double a, b, c, d;
double f(double x)
{
return (a*x*x*x + b * x*x + c * x + d);
}
int main()
{
double x1, x2, xx;
cin >> a >> b >> c >> d;
for (int x = -100; x <= 100; x++)
{
x1 = x; x2 = x + 1;//确定根可能所在的区间
if (f(x1) == 0)printf("%.2f ", x1); //考虑根在区间端点上的情况
else if (f(x1)*f(x2)<0)//如果小于0,则确定根在区间[x1,x2]中
{
while (x2 - x1 >= 0.001) //二分法确定根的值,由于要精确到小数点后2位,所以这里要算到第三位
{
xx = (x1 + x2) / 2;
if ((f(x1)*f(xx)) <= 0)x2 = xx;
else x1 = xx;
}
printf("%.2f ", x1);
}
}
return 0;
}