【codevs1038 一元三次方程求解】分治

这道题WA了将近两天！！！！！我最近一写题就WA，一WA就一天！！！！简直有毒。

题目里说了根的范围是$(-100, 100)$，而且根与根之差的绝对值>=1，那么思想就是把$(-100, 100)$分成$[-100, -99]$、$[-99, -98]$……这样一个一个的区间，看看区间端点a、b，f(a) * f(b)是否<0，是的话就对这个区间二分，找出答案。

这道题不能直接二分，因为它不是一个单调区间。超超超超超级棒的ilern跟我说如果是单调区间就可以直接二分，如果不是的话就把它划分成几个单调区间来二分。

知道了思想之后还是WA，原因是double用的不熟练。
double类型，输入用%lf， 输出是%.2lf, 2是要求保留2位小数。
float用%f输入，%.2f输出。
printf（“%03d”， 25）的结果是025，%3d的结果是25。就可以控制输出的格式。

需要注意的是用double需要注意精度，比较和常数的大小啊啥的都要允许有误差。可以定义 const double EPS = 1e-7；

#include <cstdio>
#include <cmath>

const double EPS = 1e-7;

double a, b, c, d, x1, x2, x3, l[4], r[4];

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

double find(double l, double r)
{
    //printf("find(%lf, %lf)\n", l, r);
    if (fabs(l - r) <= EPS || fabs(f((l + r) / 2)) <= EPS)
    {
        return (l + r) / 2;
    }
    else
    {
        if (f(l) * f((l + r) / 2) <= EPS)
        {
            return find(l, (l + r) / 2);
        }
        else
        {
            return find((l + r) / 2, r);
        }
    }
}

bool have(double a, double b)
{
    //double fa = f(a), fb = f(b);
    //printf("f(a = %.3lf) = %.3lf, f(b = %.3lf) = %.3lf\n", a, fa, b, fb);
    if (f(a) * f(b) <= EPS)
    {
        return true;
    }
    else
    {
        return false;
    }
}

int main()
{
    scanf("%lf %lf %lf %lf", &a, &b, &c, &d);
    int cnt = 0;

    for (int i = -100; i < 100; i++)
    {
        double j = i + 1 - EPS;
        if (have(i, j) == true)
        {
            l[cnt] = i;
            r[cnt] = j;
            cnt++;
        }
    }

    x1 = find(l[0], r[0]);
    x2 = find(l[1], r[1]);
    x3 = find(l[2], r[2]);

    printf("%.2lf %.2lf %.2lf\n", x1, x2, x3);
    return 0;
}

第一种。。。。。。我自己的bug找不出来，然后很表脸的去找了【dalao】debug，然后输出中间变量的调试方法get√【最近用gdb就只会设断点，跑跑跑，都不会输出中间变量了= =】
有的时候要加fabs【实数去绝对值，整数是abs】【例如f(d) == 0就要改成fabs(f(d)) <= EPS】，但是有的时候不用加。

#include <cstdio>
#include <cmath>

const double EPS = 1e-7;
double a, b, c, d, x1, x2, x3;

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

bool have(double a, double b)
{
    // double fa = f(a), fb = f(b);
    // printf("f(a = %.3lf) = %.3lf, f(b = %.3lf) = %.3lf\n", a, fa, b, fb);

    if (f(a) * f(b) <= EPS)
    {
        return true;
    }
    else
    {
        return false;
    }
}

void find(double l, double r)
{
    //printf("find(%lf, %lf)\n", l, r);
    double d = (l + r) / 2;
    if ((r - l) <= EPS || fabs(f(d)) <= EPS)
    {
        printf("%.2lf ", d);
        return;
    }
    else
    {
        if (f(l) * f(d) <= EPS)
        {
            find(l, d);
        }
        else
        {
            find(d, r);
        }
    }
}

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

    for (int i = -100; i < 100; i++)
    {
        if (have(i, i + 1 - EPS) == true)
        {
            find(i, i + 1 - EPS);
            //printf("%d\n", i);
        }
        else continue;
    }
    return 0;
}

第二种。。。。。。【然而并没有太大的差别(:з」∠)】
这个WA了很久的原因是主函数里i在循环的时候如果解是整数，就当i= i和i+ 1的时候都符合have，它都会输出。然后把边界定成（i， i + 1 - EPS）就行了。

【我还是太差了，这道并不算很难的题写了好久，有的bug自己还看不出来qnq】