noip2001 一元三次方程 （枚举+二分）

P1116一元三次方程求解

Accepted

标签： 搜索 枚举 NOIP提高组2001

描述

有形如：ax^3+bx^2+cx+d=0 这样的一个一元三次方程。给出该方程中各项的系数(a，b，c，d 均为实数)，并约定该方程存在三个不同实根(根的范围在-100至100之间)，且根与根之差的绝对值>=1。要求由小到大依次在同一行输出这三个实根(根与根之间留有空格)，并精确到小数点后2位。

格式

输入格式

输入该方程中各项的系数(a，b，c，d 均为实数)，

输出格式

由小到大依次在同一行输出这三个实根(根与根之间留有空格)，并精确到小数点后2位。

样例1

样例输入1[复制]

1 -5 -4 20

样例输出1[复制]

-2.00 2.00 5.00

限制

每个测试点1s

提示

提示：记函数f(x)，若存在2个实数x1和x2，且x1<x2，使f(x1)*(x2)<0，则在(x1, x2)之间一定存在实数x0使得f(x0)=0。

来源

NOIP2001第一题

解析：枚举+二分。主要是对枚举区间的端点要特别处理。

         ps：吐槽一下，我交vijos与tyvj时都出现错误，交了n都不不过，幸好tyvj提交时还能看到部分标准答案与自己的错误输出。我就发现不知为什么我用printf（）输出的都是0，直接输出答案都不行。后来改成cout就ok了。。。。。。。

              我*****

            附：已解决，double的输出采用%f即可（用lf用了4年，今天才遇到这种问题，是我老了？）

代码：

#include<cstdio>
#include<cmath>
using namespace std;

const double precision=0.000001;
double a,b,c,d;

double f(double x)
{
  return a*x*x*x+b*x*x+c*x+d;
}

double get(double l,double r)
{
  double mid;
  while(r-l>=0.001)
    {
      mid=(l+r)/2;
      if(f(l)*f(mid)<=precision)r=mid;
      else l=mid;
	}
  return l;
}

int main()
{
  scanf("%lf%lf%lf%lf",&a,&b,&c,&d);
  int k=0; double ans[3],x=-100.0;
  while(x<100.5 && k<3)
    {
      if(abs(f(x))<=precision)
	    {
		  ans[k++]=x,x+=1;
		  continue;
		}
	  if(abs(f(x+1))<=precision)
	    {
		  ans[k++]=x+1,x+=2;
		  continue;
		}
	  if(f(x)*f(x+1)<0)ans[k++]=get(x,x+1);
	  x+=1;
	}
   printf("%.2f %.2f %.2f\n",ans[0],ans[1],ans[2]);
  return 0;
}