NOIP--一元三次方程（二分）

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

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

输入输出格式

输入格式：

一行，4个实数A，B，C，D。

输出格式：

一行，三个实根，并精确到小数点后2位。

输入输出样例

输入样例#1：  复制

1 -5 -4 20

输出样例#1：  复制

-2.00 2.00 5.00

题解：因为解之间间隔大于1，所以可以从-200枚举整数到200，若f(i)*f(i+1)<0,则该区间内必有一个解，在此区间中二分查找符合值便可。复杂度O(400*logn)。

#include<stdio.h>
#include <algorithm>
#include<iostream>
#include<string.h>
#include<vector>
#include<stdlib.h>
#include<math.h>
#include<queue>
#include<deque>
#include<ctype.h>
#include<map>
#include<set>
#include<stack>
#include<string>
#include<algorithm>
#define INF 0x3f3f3f3f
#define FAST_IO ios::sync_with_stdio(false)
const double PI = acos(-1.0);
const double eps = 1e-6;
const int MAX=1e5+10;
const int mod=1e9+7;
typedef long long ll;
using namespace std;
#define gcd(a,b) __gcd(a,b)
inline ll lcm(ll a,ll b){return a/gcd(a,b)*b;}
inline ll qpow(ll a,ll b){ll r=1,t=a; while(b){if(b&1)r=(r*t)%mod;b>>=1;t=(t*t)%mod;}return r;}
inline ll inv1(ll b){return qpow(b,mod-2);}
inline ll exgcd(ll a,ll b,ll &x,ll &y){if(!b){x=1;y=0;return a;}ll r=exgcd(b,a%b,y,x);y-=(a/b)*x;return r;}
inline ll read(){ll x=0,f=1;char c=getchar();for(;!isdigit(c);c=getchar()) if(c=='-') f=-1;for(;isdigit(c);c=getchar()) x=x*10+c-'0';return x*f;}

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

int main()
{

    cin>>a>>b>>c>>d;

    for(int i=-100;i<100;i++)
    {
        double l=i,r=i+1,m,x1=f(i),x2=f(i+1);
        if(!x1) printf("%.2f ",l);

        if(x1*x2>=0) continue;

        while(l<=r-0.001)
        {
            m=(l+r)/2;
            if(f(m)*f(r)<=0)
                l=m;
            else
                r=m;
        }

        printf("%.2f ",r);
    }

    putchar('\n');
    return 0;
}