多项式求解【模板】

#include<iostream>
#include<cstdio>
#include<cstring>
#include<string.h>
#include<algorithm>
#include<vector>
#include<cmath>
#include<stdlib.h>
#include<time.h>
#include<stack>
#include<set>
#include<map>
#include<queue>
#include<sstream>
using namespace std;
#define rep0(i,l,r) for(int i = (l);i < (r);i++)
#define rep1(i,l,r) for(int i = (l);i <= (r);i++)
#define rep_0(i,r,l) for(int i = (r);i > (l);i--)
#define rep_1(i,r,l) for(int i = (r);i >= (l);i--)
#define MS0(a) memset(a,0,sizeof(a))
#define MS_1(a) memset(a,-1,sizeof(a))
#define MSinf(a) memset(a,0x3f,sizeof(a))
#define MSfalse(a) memset(a,false,sizeof(a))
#define pin1(a) printf("%d",(a))
#define pin2(a,b) printf("%d%d",(a),(b))
#define pin3(a,b,c) printf("%d%d%d",(a),(b),(c))
#define pll1(a) printf("%lld",(a))
#define pll2(a,b) printf("%lld%lld",(a),(b))
#define pll3(a,b,c) printf("%lld%lld%lld",(a),(b),(c))
#define pdo1(a) printf("%f",(a))
#define pdo2(a,b) printf("%f%f",(a),(b))
#define pdo3(a,b,c) printf("%f%f%f",(a),(b),(c))
#define huiche puts("")
#define inf 0x3f3f3f3f
#define lson i<<1,l,mid
#define rson ((i<<1)|1),mid+1,r
#define uint unsigned int
typedef pair<int,int> PII;
#define A first
#define B second
#define pb push_back
#define MK make_pair
#define ll long long
#define eps 1e-8
inline void read1(int &num) {
    char in;
    bool IsN=false;
    in=getchar();
    while(in!='-'&&(in<'0'||in>'9')) in=getchar();
    if(in=='-') {
        IsN=true;
        num=0;
    } else num=in-'0';
    while(in=getchar(),in>='0'&&in<='9') {
        num*=10,num+=in-'0';
    }
    if(IsN) num=-num;
}
inline void read2(int &a,int &b) {
    read1(a);
    read1(b);
}
inline void read3(int &a,int &b,int &c)
{
    read1(a);
    read1(b);
    read1(c);
}
inline void read1(ll &num) {
    char in;
    bool IsN=false;
    in=getchar();
    while(in!='-'&&(in<'0'||in>'9')) in=getchar();
    if(in=='-') {
        IsN=true;
        num=0;
    } else num=in-'0';
    while(in=getchar(),in>='0'&&in<='9') {
        num*=10,num+=in-'0';
    }
    if(IsN) num=-num;
}
inline void read2(ll &a,ll &b) {
    read1(a);
    read1(b);
}
inline void read3(ll &a,ll &b,ll &c)
{
    read1(a);
    read1(b);
    read1(c);
}
inline void read1(double &num) {
    char in;
    double Dec=0.1;
    bool IsN=false,IsD=false;
    in=getchar();
    while(in!='-'&&in!='.'&&(in<'0'||in>'9'))
        in=getchar();
    if(in=='-') {
        IsN=true;
        num=0;
    } else if(in=='.') {
        IsD=true;
        num=0;
    } else num=in-'0';
    if(!IsD) {
        while(in=getchar(),in>='0'&&in<='9') {
            num*=10;
            num+=in-'0';
        }
    }
    if(in!='.') {
        if(IsN) num=-num;
        return ;
    } else {
        while(in=getchar(),in>='0'&&in<='9') {
            num+=Dec*(in-'0');
            Dec*=0.1;
        }
    }
    if(IsN) num=-num;
}
inline void read2(double &a,double &b)
{
    read1(a);
    read1(b);
}
inline void read3(double &a,double &b,double &c)
{
    read1(a);
    read1(b);
    read1(c);
}















int xi[1000];
int new_xi[1000];
int maxm;
double ep;
void init(void)
{
    for(int i = 0;i<maxm-1;i++)
    {
        new_xi[i] = xi[i]*maxm-i;
    }
}
double ff(double x)
{//求导
    double ret = new_xi[0];
    for(int i = 1;i<maxm-1;i++)
        ret = ret*x+xi[i];
    return ret;
}
double f(double x)
{
    //x*x+x*2 +1
    double ret = xi[0];
    for(int i = 1;i<maxm;i++)
        ret = ret*x+xi[i];
    return ret;
}
void ans(void)
{
    for(int i = 0;i<1000;i++)
    {
        bool flag = false;
        double ans = rand()/32767.0*(2e9)-1e9;
        double ans1 = ans - f(ans)/ff(ans);
        while(fabs(ans1-ans)>ep)
        {
            ans = ans1;
            ans1 = ans - f(ans)/ff(ans);
        }
        printf("%.4f\n",ans);
    }
}
int main()
{
    freopen("in.txt","r",stdin);
    srand(time(0));
    maxm = 0;
    puts("请分别输入精度，系数：");
    read1(ep);
    while(scanf("%d",&xi[maxm])!=EOF) maxm++;
    init();
    ans();
    return 0;
}