我觉得我考的话最多敲个高精了……
n=0 => n%p=0,显然是对的。
n%p=0 => n=0,显然是错的。
但p多选几个,就极有可能是对的……
所以说我们可以对输入给定的数模好几个质数,然后枚举1~m看看有几个能使方程成立,这样能拿70 。
满分做法:n%p=(n+kp)%p,所以每次只枚举[0,p-1],然后把倍数筛掉即可。
p不能选太大,否则会T。p不能选太小,否则会WA。p不能选太少,否则筛不干净……
总之我这样的脸黑的,做这种题简直丧病……
还有,这个题实测加不加秦九韶好像差不多
代码:
#include<cstdio>
#include<iostream>
#include<cstring>
#include<algorithm>
using namespace std;
const int SIZE=2000010;
typedef long long LL;
LL p[]={11279,11953,12133,12197,13831};
LL a[SIZE][10];
void add(LL &x,LL mod,char c)
{
x=(((x<<3)%mod+(x<<1)%mod)%mod+c-'0')%mod;
}
void change(LL &x,LL mod)
{
x=((-x)%mod+mod)%mod;
}
void scan(LL i)
{
char c=getchar();
bool flag=0;
while(c<'0'||c>'9') {if(c=='-') flag=1;c=getchar();}
while(c>='0'&&c<='9')
{
for(int j=0;j<5;j++)
add(a[i][j],p[j],c);
c=getchar();
}
if(flag)
{
for(int j=0;j<5;j++)
change(a[i][j],p[j]);
}
}
int n,m;
LL f(LL x,LL id)
{
LL mod=p[id];
LL ans=a[n][id];
for(int i=n;i>=1;i--)
{
ans=( ( ( (ans%mod)* (x%mod) )%mod +a[i-1][id]%mod) %mod)%mod;
}
return ans;
/* LL mod=p[id],ans=0,xx=1;
for(int i=0;i<=n;i++)
{
ans=(ans+((a[i][id]%mod)*(xx%mod))%mod)%mod;
xx=((xx%mod)*(x%mod))%mod;
}
return ans;*/
}
bool use[SIZE];
int ans[SIZE];
int main()
{
scanf("%d%d",&n,&m);
for(int i=0;i<=n;i++)
scan(i);
for(int j=0;j<5;j++)
{
for(int i=0;i<p[j];i++)
{
if(f(i,j)!=0)
{
for(int k=i;k<=m;k+=p[j])
use[k]=1;
}
}
}
for(int i=1;i<=m;i++)
{
if(!use[i])
ans[++ans[0]] = i;
}
for(int i=0;i<=ans[0];i++)
{
printf("%d\n",ans[i]);
}
return 0;
}