一元三次方程求解
来源:互联网 发布:nativeshare.js下载 编辑:程序博客网 时间:2024/04/27 19:40
Problem Description
有形如:ax^3+bx^2+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)之间一定有一个根。
要求由小到大依次在同一行输出这三个实根(根与根之间留有空格),并精确到小数点后2位。
提示:记方程f(x)=0,若存在2个数x1和x2,且x1<x2,f(x1)*f(x2)<0,则在(x1,x2)之间一定有一个根。
Input
输入有多行测试数据,每行为四个系数a,b,c,d,输入以0 0 0 0结束。
Output
对于每组测试数据,输出一个三个实根(根与根之间留有空格),并精确到小数点后2位。
Sample Input
1 -5 -4 200 0 0 0
Sample Output
-2.00 2.00 5.00
/* 利用公式把x1,x2求出。在(-100,x1)、(x1,x2)、(x2,100)范围内找出三个不同实根 利用二分法求出三个不同实根*/#include <stdio.h>#include <math.h>double a,b,c,d;double fun(double x){ return a*x*x*x+b*x*x+c*x+d;}double find(double min,double max) // 二分法{ double mid; while(1) { mid=(min+max)/2; if( fabs(min-max)<1e-8 ) return mid; else{ if(fun(mid)==0) return mid; else if( fun(min) * fun(mid)<0 ) max=mid; else min=mid;} }}int main(){ double x1,x2; while(scanf("%lf%lf%lf%lf",&a,&b,&c,&d)!=EOF) { x1=(-2*b-sqrt(4*b*b-12*a*c))/(a*6); // 利用公式把x1,x2求出 x2=(-2*b+sqrt(4*b*b-12*a*c))/(a*6); printf("%.2lf %.2lf %.2lf\n",find(-100,x1),find(x1,x2),find(x2,100)); } return 0;}
0 0
- 一元三次方程求解
- 一元三次方程求解
- 一元三次方程求解
- 一元三次方程求解
- 一元三次方程求解
- 一元三次方程求解
- 一元三次方程求解
- 一元三次方程求解
- 一元三次方程求解
- 一元三次方程求解
- 一元三次方程求解
- 一元三次方程求解
- 一元三次方程求解
- 一元三次方程的求解
- 求解一元三次方程--c++
- C++求解一元三次方程
- NOIP2011 一元三次方程求解
- NOIP2001 一元三次方程求解
- sendRedirect&forword的区别2
- HDU1162 Eddy's picture 【最小生成树Prim】
- X86和ARM:函数调用参数传递和返回值
- webservice应用举例
- binder的机制和原理
- 一元三次方程求解
- bluedroid 协议 学习
- poj 1258 Agri-Net (最小生成树 prim)
- 多线程和多进程的区别(小结)
- Android ActionBar完全解析,使用官方推荐的最佳导航栏(下)
- hdu 1372 || poj 2243 Knight Moves(搜索:BFS+优先队列水题)
- jquery1.11源码学习——揭秘jquery的extend方法
- 关于C#字符串的截取
- hdoj.1286 找新朋友 20140724