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);
}