NOI:7891 一元三次方程求解

题目链接


题解:二分法求根,浮点数精度

#include <stdio.h>
double a,b,c,d;
inline double sol(double x){
    return a*x*x*x+b*x*x+c*x+d;
}
double query(double l,double r,bool f){//二分法求根
    double mid=0.00;
    while(l+0.0001<r){//精度
        mid=(l+r)/2;
        bool b=sol(mid)>0;
        if(f?(!b):b)r=mid;//二分法
        else l=mid;
    }
    return mid;
}
int main(){
    
    scanf("%lf%lf%lf%lf",&a,&b,&c,&d);
    if(a>0){
        double i;
        for(i=-100;i<=100&&sol(i)<0;i++);
        i+=0.00001;//精度问题 得加这一句 不然判断时可能sol(i)<0==1或==0
        printf("%.2lf ",query(i-1,i,false));//up
        for(;i<=100&&sol(i)>=0;i++);
        i+=0.00001;
        printf("%.2lf ",query(i-1,i,true));//down
        for(;i<=100&&sol(i)<0;i++);
        i+=0.00001;
        printf("%.2lf\n",query(i-1,i,false));//up
    }else {
        double i;
        for(i=-100;i<=100&&sol(i)>=0;i++);
        i+=0.00001;
        printf("%.2lf ",query(i-1,i,true));//down
        for(;i<=100&&sol(i)<0;i++);
        i+=0.00001;
        printf("%.2lf ",query(i-1,i,false));//up
        for(;i<=100&&sol(i)>=0;i++);
        i+=0.00001;
        printf("%.2lf\n",query(i-1,i,true));//down
    }
    
}

猜你喜欢

转载自blog.csdn.net/wuzhenzi5193/article/details/80587160