HLG 1215 一元三次方程求解
来源:互联网 发布:java getclassname 编辑:程序博客网 时间:2024/05/16 11:51
Description
有形如:a*x^3+b*x^2+c*x+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
有多组测试数据,对于每组测试数据,仅有一行,包括四个数,a, b, c, d,相邻两个数由空格隔开。
Output
对于每组测试数据,输出三个解,要求从小到大输出,相邻两个数用空格隔开。
Sample Input
1 -5 -4 20
Sample Output
-2.00 2.00 5.00
#include <iostream>
#include <iomanip>
#include <math.h>
using namespace std;
double re[3];
double a,b,c,d,u,v;
int i,t;
double f(double x)
{
return ((a*x+b)*x+c)*x+d;
}
int main()
{
while(cin>>a>>b>>c>>d)
{
t=-1;
for(i=-100;i<=100;i++)
{
u=double(i);
v=u+0.99999;
if(fabs(f(u))<0.00001||f(u)*f(v)<=0)
{
t++;
if(fabs(f(u))<0.00001)
re[t]=u;
else
{
while((u+0.001<v)&&fabs(f((u+v)/2))>=0.00001)
{
if(f(u)*f((u+v)/2)<0)
v=(u+v)/2;
else
u=(u+v)/2;
}
re[t]=(u+v)/2;
}
}
}
cout<<setiosflags(ios::fixed)<<setprecision(2);
cout<<re[0]<<" "<<re[1]<<" "<<re[2]<<endl;
}
return 0;
}
- HLG 1215 一元三次方程求解
- 一元三次方程求解
- 一元三次方程求解
- 一元三次方程求解
- 一元三次方程求解
- 一元三次方程求解
- 一元三次方程求解
- 一元三次方程求解
- 一元三次方程求解
- 一元三次方程求解
- 一元三次方程求解
- 一元三次方程求解
- 一元三次方程求解
- 一元三次方程求解
- 一元三次方程的求解
- 求解一元三次方程--c++
- C++求解一元三次方程
- NOIP2011 一元三次方程求解
- 报道
- 新的起点
- 青青购返利网 www.51shzyt.com:网络购物背后孕育新的商业模式
- Resource '/servers' does not exist 问题的解决
- 我要改变
- HLG 1215 一元三次方程求解
- Windows Mobile6.5系统 SYS移植精简
- Windows Mobile DIY ROM 常用工具
- linux学习笔记(更新中。。。)
- Apktool v1.5.2 - APK反编译工具
- 【爱机控首发】HTC_Desire_606w_1.…
- iOS7不越狱干掉新浪微博客户端广告…
- Ubuntu 12.04创建快捷方式
- 【Android Kitchen v0.227】ROM定制厨房越狱修改汉化版