noip2014 解方程 Hash+数学

f(x)=0 ==> f(x)%p=0。但是f(x)%p=0不能推出f(x)=0。
如果对于许多个不同的p，f(x)%p=0 ，那么可看做f(x)=0；
所以选许多个（3个以上）质数p，f(x)都=0，那么x可以认为就是答案了。
枚举x，O（N）计算f（x）的值，时间复杂度O（N*M）；
虽然说只要机子吊，十亿随便跑，但是CCF的算盘对这个复杂度还是很难接受的。
如果f（x）%p！=0，那么f（x+p）%p！=0，所以这是一个剪枝……

#include<iostream>
#include<cstdio>
#define LL long long
using namespace std;
LL pri[]={0,100003,63443,79087,99901,143651};
int T=5;
LL a[105][10];
bool use[1000007];
bool can[1000007][10];
void read(LL i)
{
    LL f=1;
    char c=getchar();
    while(c<'0'||c>'9')
    {
        if(c=='-') f=-1;
        c=getchar();
    }
    while(c>='0'&&c<='9')
    {
        for(int k=1;k<=T;k++)
        {
            a[i][k]*=10;
            a[i][k]+=c-'0';
            a[i][k]%=pri[k];
        }
        c=getchar();
    }
    for(int k=1;k<=T;k++)
        a[i][k]*=f;
}
int n,m;
bool f(LL x)
{
    bool b=0;
    for(int k=1;k<=T;k++)
    {
        LL ans=0;
        for(int i=n;i>=1;i--)
        {
            ans+=a[i][k];
            ans*=x;ans%=pri[k];
        }
        ans+=a[0][k];
        ans%=pri[k];
        if(ans!=0) 
        {
            can[x][k]=1;
            b=true;
        }
    }
    if(b) return false;
    return true;
}
int main()
{
    scanf("%d%d",&n,&m);
    for(int i=0;i<=n;i++)
        read(i);
    int cnt=0;
    for(LL i=1;i<=m;i++)
    {
        if(!use[i])
        {
            bool h=0;
            for(int k=1;k<=T;k++)
            {
                LL p=i%pri[k];
                if(can[p][k]==1)
                {use[i]=1;h=1;break;}
            }
            if(h) continue;
            if(f(i)==false)
                use[i]=1;
            else
                cnt++;
        }
    }
    printf("%d\n",cnt);
    for(int i=1;i<=m;i++)
        if(!use[i]) printf("%d\n",i);
}