首页 > 代码库 > XidianOJ 1093 一元三次方程
XidianOJ 1093 一元三次方程
题目描述
有形如:ax3+bx2+cx+d=0 这样的一个一元三次方程。给出该方程中各项的系数(a,b,c,d 均为实数),并约定该方程存在三个不同实根(根的范围在-100至100之间),且根与根之差的绝对值>=1。要求由小到大依次在同一行输出这三个实根(根与根之间留有空格),并精确到小数点后2位。
输入
多组数据,一行,四个数
输出
三个根
样例输入
1 -5 -4 20
样例输出
-2.00 2.00 5.00
--正文
告诉你,根在-100 ~ 100间,2分找
#include<iostream> #include<cstdio> #include<cstring> #include<algorithm> #include<cmath> using namespace std; const double eps=0.0001; double a,b,c,d; double ans[5]; int sum; double calc(double x) { return a*x*x*x+b*x*x+c*x+d; } void find(double l,double r) { if(r-l <= eps) { ans[++sum] = (l+r)/2; return; } double mid = (l+r)/2,tmp = calc(mid),tmpl = calc(l),tmpr = calc(r); if(abs(tmp) < eps) { ans[++sum] = mid; return; } if(tmp*tmpl < 0) { find(l,mid); return; } find(mid,r); } int main() { while(scanf("%lf %lf %lf %lf",&a,&b,&c,&d) != EOF) { sum = 0; for(int i=-101;i<=101;i++) { double tmp1 = calc(i),tmp2 = calc(i+1); if(abs(tmp1) < eps) ans[++sum] = i; else if(abs(tmp2) < eps) { ans[++sum] = i+1; i++; } else if(tmp1*tmp2 < 0) find(i,i+1); if(sum == 3) break; } sort(ans+1,ans+4); printf("%.2lf %.2lf %.2lf\n",ans[1],ans[2],ans[3]); } return 0; }
XidianOJ 1093 一元三次方程
声明:以上内容来自用户投稿及互联网公开渠道收集整理发布,本网站不拥有所有权,未作人工编辑处理,也不承担相关法律责任,若内容有误或涉及侵权可进行投诉: 投诉/举报 工作人员会在5个工作日内联系你,一经查实,本站将立刻删除涉嫌侵权内容。