不确定用二分写行不行。担心会不会卡精度,于是就设置一下精度下限就过了。
1038 一元三次方程求解
2001年NOIP全国联赛提高组
时间限制: 1 s
空间限制: 128000 KB
题目等级 : 白银 Silver
查看运行结果
题目描述 Description
有形如:ax3+bx2+cx+d=0 这样的一个一元三次方程。给出该方程中各项的系数(a,b,c,d 均为实数),并约定该方程存在三个不同实根(根的范围在-100至100之间),且根与根之差的绝对值>=1。要求由小到大依次在同一行输出这三个实根(根与根之间留有空格),并精确到小数点后2位。
提示:记方程f(x)=0,若存在2个数x1和x2,且x1<x2,f(x1)*f(x2)<0,则在(x1,x2)之间一定有一个 根。
输入描述 Input Description
一个三次方程的各项系数
输出描述 Output Description
三个解
样例输入 Sample Input
1 -5 -4 20
样例输出 Sample Output
-2.00 2.00 5.00
----------------------------------------------------------------------------
#include<cstdio>
#include<cstring>
#include<iostream>
#include<algorithm>
using namespace std;
double ans[4];
double a,b,c,n;
int cur=0;
void pd(double l,double r){
if(r-l<0.0001) return;
double m=(l+r)/2.0;
double s=l*l*l*a+l*l*b+l*c+n;
double t=r*r*r*a+r*r*b+r*c+n;
if(s*t<0){
if(r-l<0.001){
ans[++cur]=l;
return ;
}
}
pd(l,m);
pd(m,r);
return;
}
int main(){
scanf("%lf%lf%lf%lf",&a,&b,&c,&n);
pd(-100.0,100.0);
printf("%.2lf %.2lf %.2lf\n",ans[1],ans[2],ans[3]);
return 0;
}
----------------------------------------------------------------------------