一元三次方程组求解 luogu P1024

题目传送门

首先，要明确题目信息，f(x1) * f(x2) < 0, 则一定存在实数根在区间(x1, x2)。且所有的根都在[-100, 100)之间。根与根的绝对值之差 >= 1

那么，我们是否可以找到所有的x1 和 x2 呢？   当然可以。

根的取值范围这么小 （每一个区间都枚举一次的话也只有200个区间），那么我们当然就愉快的枚举啊！

在 [-100, 100)的区间里面枚举符合条件的 x1, x2。然后再对该区间进行二分查找（二分赛高！）

#include <bits/stdc++.h>
using namespace std;
#define du double 

du a, b, c, d;
int now;
du ans;

inline du cac(du x){ //x 对应的 f 的值
    return a * x * x * x + b * x * x + c * x + d;
}

inline void search(du l, du r){

    if(now == )return;
    if(r - l <= 0.001){//差这么小，已经可以说是相等了。直接输出
        printf("%.2lf ", l);
        now++;
        return ;
    }
    du mid = l + (r - l) / ;
    du ansl, ansr;
    ansl = cac(l) * cac(mid);
    ansr = cac(mid) * cac(r);
    if(cac(r) == ){     //是否存在端点为根的情况？
        printf("%.2lf ", r);
        now++;
    }
    if(cac(mid) == ){
        printf("%.2lf ", mid);
        now++;
    }
    if(ansl < )
        search(l, mid);//继续二分搜索， 注意不一定只有一个根
    else if(ansr < )search(mid, r);
    return;
}

int main(){
    freopen("count.in", "r", stdin);
    freopen("count.out", "w", stdout);
    cin >> a >> b >> c >> d;
    for(du i = -100.0; i <= ; i++){
        if(cac(i) * cac(i + ) <= )//枚举符合条件的区间
            search(i, i + );
    }
    return ;
}