【noip】解方程 秦九韶算法

解方程

描述

已知多项式方程：
$a_0+a_1x^1+a_2x^2……+a_{n-1}x^{n-1}+a_nx^n=0$
求这个方程在[1, m]内的整数解（n 和 m 均为正整数）。

输入格式

输入共 n+2 行。
第一行包含 2 个整数 n、m，每两个整数之间用一个空格隔开。
接下来的 n+1 行每行包含一个整数，依次为$a_0,a_1,a_2,...,a_n$。

输出格式

第一行输出方程在[1, m]内的整数解的个数。
接下来每行一个整数，按照从小到大的顺序依次输出方程在[1, m]内的一个整数解。

样例输入1

2 10
1
-2
1

样例输出1

1
1

样例输入2

2 10
2
-3
1

样例输出2

2
1
2

样例输入3

2 10
1
3
2

样例输出3

0

限制

对于 30%的数据，0 < n ≤ 2， |$a_i$| ≤ 100，anan ≠ 0， m ≤ 100；
对于 50%的数据，0 < n ≤ 100，|$a_i$| ≤ $10^{100}$ ，$a_n$ ≠ 0，m ≤ 100；
对于 70%的数据，0 < n ≤ 100，|$a_i$| ≤ $10^{10000}$，$a_n$ ≠ 0，m ≤ 10000；
对于 100%的数据，0 < n ≤ 100，|$a_i$| ≤ $10^{10000}$，$a_n$ ≠ 0，m ≤ 1000000。

来源

NOIP2014 提高组 Day2

这个题一看很吓人，$10^{10000}$？！多大一个数啊，高精度？肯定超时.
怎么做？
先是一定要用秦九韶算法，不然写起来不仅慢还很辛苦。

秦九韶算法：

$f(x)=a_0+a_1x^1+a_2x^2……+a_{n-1}x^{n-1}+a_nx^n=a_0+x(a_1+x(a_2+x(a_3+……x(a_{n-1}+xa_n))…)$
这样处理以后，枚举x，从内向外算就可以在O（n）验证方程，而且写起来也轻松
但是数字很大，高精度不可能
写不来set的我，天天暴力哈希，这让我直接就想到可以都同余方程，两边各模一个质数，若原方程计算结果为0，模了以后也为0，然后多取几个质数可以避免冲突。过了70%，T了30%。
也就是说还要优化。
都想到了同余方程那自然就可以想到后面的这个了：$f(x)\equiv f(x+p)(mod p)$
这样我们就只用算1~p而不是1~m了，p取3~5个2万左右的质数就可以避免冲突了。
代码见下，很简单，就没有写注释了。

#include<iostream>
#include<cstdio>
#include<string>
#define ll long long

using namespace std;
ll a[4][105],x,m,n,k,b[4][1000005],p[4]={0,10007,23333,11261};
int main()
{
    cin>>n>>m;
    for(int i=0;i<=n;i++)
    {
        string s;
        cin>>s;
        int len=s.size(),flag=0;
        if(s[0]=='-')flag=1;
        for(int j=flag;j<len;j++)
        {
            a[1][i]=(a[1][i]*10+s[j]-'0')%p[1];
            a[2][i]=(a[2][i]*10+s[j]-'0')%p[2];
            a[3][i]=(a[3][i]*10+s[j]-'0')%p[3]; 
        }
        if(flag)
        {
            a[1][i]=-a[1][i];
            a[2][i]=-a[2][i];
            a[3][i]=-a[3][i];
        }
    }
    for(int k=1;k<=3;k++)
    {
        for(int i=1;i<p[k];i++)
        {
            int t=a[k][n];
            for(int j=n-1;j>=0;j--)
                t=(t*i+a[k][j])%p[k];
            if(t==0)b[k][i]=1;
        }
    }
    for(int i=1;i<=m;i++)
    {
        b[0][i]=1;
        for(int j=1;j<=3;j++)
        if(!b[j][i%p[j]])b[0][i]=0;
    }
    int num=0;
    for(int i=1;i<=m;i++)
        if(b[0][i])num++;
    cout<<num<<'\n';
    for(int i=1;i<=m;i++)
        if(b[0][i])printf("%d\n",i);
    return 0;
}

大概就是这个样子，如果有什么问题，或错误，请在评论区提出，谢谢。