P1024 [NOIP2001 提高组] 一元三次方程求解

题目描述

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

提示：记方程 �(�)=0f(x)=0，若存在 22 个数 �1x1​ 和 �2x2​，且 �1<�2x1​<x2​，�(�1)×�(�2)<0f(x1​)×f(x2​)<0，则在 (�1,�2)(x1​,x2​) 之间一定有一个根。

输入格式

一行，44 个实数 �,�,�,�a,b,c,d。

输出格式

一行，33 个实根，从小到大输出，并精确到小数点后 22 位。

输入输出样例

输入 #1复制

1 -5 -4 20

输出 #1复制

-2.00 2.00 5.00

说明/提示

【题目来源】

NOIP 2001 提高组第一题

【二分】
因为区间很大，所以可以二分。
三个答案都在[-100,100]范围内，两个根的差的绝对值>=1,保证了每一个大小为1的区间里至多有1个解，也就是说当区间的两个端点的函数值异号时区间内一定有一个解，同号时一定没有解。那么我们可以枚举互相不重叠的每一个长度为1的区间，在区间内进行二分查找。
【参见代码】

#include<cstdio>
double a,b,c,d;
double fc(double x)
{
    return a*x*x*x+b*x*x+c*x+d;
}
int main()
{
    double l,r,m,x1,x2;
    int s=0,i;
    scanf("%lf%lf%lf%lf",&a,&b,&c,&d);  //输入
    for (i=-100;i<100;i++)
    {
        l=i; 
        r=i+1;
        x1=fc(l); 
        x2=fc(r);
        if(!x1) 
        {
            printf("%.2lf ",l); 
            s++;
        }      //判断左端点，是零点直接输出。
                        
                        //不能判断右端点，会重复。
        if(x1*x2<0)                             //区间内有根。
        {
            while(r-l>=0.001)                     //二分控制精度。
            {
                m=(l+r)/2;  //middle
                if(fc(m)*fc(r)<=0) 
                   l=m; 
                else 
                   r=m;   //计算中点处函数值缩小区间。
            }
            printf("%.2lf ",r);  
            //输出右端点。
            s++;
        }
        if (s==3) 
            break;             
            //找到三个就退出大概会省一点时间
    }
    return 0;
}