解方程

#include<cstdio>
#include<queue>
#include<cstring>
using namespace std;
#define ll long long
const int M=11000;
int n,m;ll f[100000][5];
const int p[]={10007,10917,30071}; ll a[5][110],b[1110],s[1100000],x[]={1,1,1},tot;
bool calc(int value, int j)
{
long long tmp=0;
for(int i=n;i>=0;--i)
tmp=(tmp*value+a[j][i])%p[j];
return tmp!=0;
}

int main(){
//freopen("equation.in","r",stdin);
//freopen("equation.out","w",stdout);
scanf("%d%d\n",&n,&m);
for(int i=0;i<=n;i++){
char w[M];
gets(w);
for(int j=0;j<strlen(w);j++)
if(w[j]!='-'){
for(int k=0;k<3;k++)
a[k][i]=(a[k][i]*10+w[j]-'0')%p[k];
}
for(int k=0;k<3;k++)
if(w[0]=='-')a[k][i]=p[k]-a[k][i];
}
for(int j=0;j<3;++j)
for(int i=0;i<p[j];++i)
f[i][j]=calc(i,j);
for(int i=1;i<=m;i++){
ll ans=0;
int jhfd=0;
for(int k=0;k<3;k++)
if(f[i%p[k]][k]){
jhfd=1;
}
if(!jhfd)
s[++tot]=i;
}

printf("%d\n",tot);
for(int i=1;i<=tot;i++) printf("%lld\n",s[i]);
}