luogu2312【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 longinline 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[3]={10007,10917,30071},a[N][3];bool f[40000][3];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<3;++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<3;++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<3;++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;}