一元三次方程 double输出 -0.00

求一个 a*x*x*x+b*x*x+c*x+d 的解

题目很简单，但是我输出了-0.00，然后就一直卡着，这个问题以后要注意。

让0.00 编程-0.00的方法有很多。

第一种就是直接特判

if(fabs(x-0)<eps) x=0;

第二种就是+eps，但是这样可能会有问题

第三种就是二分的时候让答案往正的输出（这个是针对我这个题目的）

我觉得第三种比较靠谱，所以就写了第三个

然后就是让这个eps的精度变高一点，这样的话，浮点误差就比较小。

#include <cstdio>

#include <iostream>

#include <cstring>

#include <algorithm>

#include <cmath>

#include <queue>

#include <map>

#include <stack>

#define inf 0x3f3f3f3f

using namespace std;

const int maxn=2e4+;

typedef long long ll;

const double eps=1e-;

double a,b,c,d;

double cal(double x){

    return a*x*x*x+b*x*x+c*x+d;

}

double ok1(double l,double r){

    while(r-l>eps){

        double mid=(l+r)/;

        if(cal(mid)<) l=mid;//!!!

        else r=mid;

    }

    return r;//print r   not l

}

double ok2(double l,double r){

    while(r-l>eps){

        double mid=(l+r)/;

        if(cal(mid)>) l=mid;//!!!

        else r=mid;

    }

    return l;//print l not r

}

int main(){

    scanf("%lf%lf%lf%lf",&a,&b,&c,&d);

    double a1=*a,b1=*b,c1=c;

    double x1=(-b1+sqrt(b1*b1-*a1*c1))/(*a1);

    double x2=(-b1-sqrt(b1*b1-*a1*c1))/(*a1);

    if(x1>x2) swap(x1,x2);

    double ans1=-,ans2=,ans3=;

    if(a>) ans1=ok1(-,x1),ans2=ok2(x1,x2),ans3=ok1(x2,);

    else ans1=ok2(-,x1),ans2=ok1(x1,x2),ans3=ok2(x2,);

    printf("%.2f %.2f %.2f\n",ans1,ans2,ans3);

    return ;

}

//1.000000 134.810000 4749.016919 49980.676731

//-0.715  -0.71

if(fabs(x-)<eps) x=;

巴特西

一元三次方程 double输出 -0.00

最新文章

热门文章