码迷,mamicode.com
首页 > 其他好文 > 详细

noip2001T1二分

时间:2015-10-27 23:25:15      阅读:247      评论:0      收藏:0      [点我收藏+]

标签:

不确定用二分写行不行。担心会不会卡精度,于是就设置一下精度下限就过了。

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;
}

----------------------------------------------------------------------------

noip2001T1二分

标签:

原文地址:http://www.cnblogs.com/fighting-to-the-end/p/4915674.html

(0)
(0)
   
举报
评论 一句话评论(0
登录后才能评论!
© 2014 mamicode.com 版权所有  联系我们:gaon5@hotmail.com
迷上了代码!