bzoj3751【2014提高】解方程（数论+筛法+枚举）

这题强啊。。。蒟蒻弱到去写带负数的高精度了，结果玩挂。
因为数无比的大，所以我们发现给两边都模个数好了hh。为了更准确些，让他们摸个大质数吧！一个正确性还不太够，多模几个吧！于是诞生了正解。读入的时候就顺道都模了，我这里是用了三个质数。然后枚举x，带入，用秦九韶公式计算，判断是否得0.如果模3个质数的情况下均为0，基本可以说他就是得0了（错误的概率极小）。然后这样枚举是O(nm)的，可能会卡常。考虑一个小优化，x,x+p,x+p*k，（p为模的质数）带入后的答案应该是相同的！因此我们只需计算出0…p-1带入后的答案是否满足即可，其他的数就像素数筛那样直接推出即可。复杂度(n*p)。素数选的好，不超ll，跑的贼快，算的贼准。素数选的不好，wawawa，TTT.

#include <bits/stdc++.h>
using namespace std;
#define N 110
#define inf 0x3f3f3f3f
#define ll long long
inline int read(){
    int x=0,f=1;char ch=getchar();
    while(ch<'0'||ch>'9'){if(ch=='-')f=-1;ch=getchar();}
    while(ch>='0'&&ch<='9') x=x*10+ch-'0',ch=getchar();
    return x*f;
}
int n,m,prime[5]={30011,11261,14843,19997,21893},a[N][5];
bool f[40000][5];
char s[10010];
vector<int>ans;
bool calc(int x,int j){
    int tmp=a[n][j];
    for(int i=n-1;i>=0;--i)
        tmp=(tmp*x+a[i][j])%prime[j];
    return !tmp;
}
int main(){
//  freopen("a.in","r",stdin);
    n=read();m=read();
    for(int i=0;i<=n;++i){
        scanf("%s",s+1);int len=strlen(s+1);
        for(int j=0;j<5;++j){
            if(s[1]=='-'){
                int x=0;
                for(int ii=2;ii<=len;++ii)
                    x=x*10+s[ii]-'0',x%=prime[j];
                a[i][j]=-x;
            }else{
                int x=0;
                for(int ii=1;ii<=len;++ii)
                    x=x*10+s[ii]-'0',x%=prime[j];
                a[i][j]=x;
            }
        }       
    }
    for(int j=0;j<5;++j)
        for(int i=0;i<prime[j];++i)
            f[i][j]=calc(i,j);
    for(int i=1;i<=m;++i){
        bool flag=1;
        for(int j=0;j<5;++j)
            if(!f[i%prime[j]][j]) flag=0;
        if(flag) ans.push_back(i);
    }
    printf("%d\n",ans.size());
    for(int i=0;i<ans.size();++i) printf("%d\n",ans[i]);
    return 0;
}