luogu1024一元三次方程求解(NOIP2001提高组第1题)

xiaoxiao2021-02-28  41

luogu1024  一元三次方程求解(NOIP2001提高组第1题)

时空限制    1000ms/128MB

题目描述

有形如: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)之间一定有一个根。

输入输出格式

输入格式:

一行,4个实数A,B,C,D。

输出格式:

一行,三个实根,并精确到小数点后2位。

输入输出样例

输入样例#1: 1 -5 -4 20 输出样例#1: -2.00 2.00 5.00

代码

#include<cstdio> 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(){ scanf("%lf%lf%lf%lf",&a,&b,&c,&d); int cnt=0; for (double x=-100,x1,x2,xx; x<=100; x+=0.9999){ //枚举可能有答案的区间 x1=x; x2=x+0.9999; if (f(x1)==0) ++cnt==3?printf("%.2lf\n",x1):printf("%.2lf ",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; } ++cnt==3?printf("%.2lf\n",x1):printf("%.2lf ",x1); } } return 0; }
转载请注明原文地址: https://www.6miu.com/read-2622532.html

最新回复(0)