BZOJ 3751: [NOIP2014]解方程 数学

3751: [NOIP2014]解方程

题目连接：

http://www.lydsy.com/JudgeOnline/problem.php?id=3751

Description

已知多项式方程：

a0+a1x+a2x^2+...+an*xn=0
求这个方程在[1,m]内的整数解（n和m均为正整数）。

Input

第一行包含2个整数n、m，每两个整数之间用一个空格隔开。
接下来的n+1行每行包含一个整数，依次为a0,a1,a2,...,an。

Output

第一行输出方程在[1,m]内的整数解的个数。

接下来每行一个整数，按照从小到大的顺序依次输出方程在[1,m]内的一个整数解。

Sample Input

2 10

2

-3

1

Sample Output

2

1

2

Hint

对于100%的数据，0<n≤100，|ai|≤1010000，an≠0，m≤1000000。

题意

题解：

找6个10000左右的素数，然后把这个素数范围内的数拿去测试就好了

因为 (a % p) = (a+p % p)，所以a不合适，那么a+p肯定也是不合适的。

于是用这个方法去xjb筛一波就好了。

代码

include<bits/stdc++.h>

using namespace std;

const int maxn = 105;
const int N = 1e6+6;
int p[10]={19997,22877,14843,10007,11261,21893};
char s[105][10005];
int flag[N];
int a[maxn],n,m,len[maxn];
int main(){
scanf("%d%d",&n,&m);
for(int i=0;i<=n;i++)
scanf("%s",s[i]);
for(int i=0;i<=n;i++)
len[i]=strlen(s[i]);
for(int o=0;o<6;o++){
for(int i=0;i<=n;i++){
int flag = (s[i][0]=='-')?1:0;
a[i]=0;
for(int j=flag;j<len[i];j++){
a[i]=(a[i]10+(s[i][j]-'0'))%p[o];
}
if(flag)a[i]=-a[i];
}
for(int i=1;i<=p[o];i++){
int tmp=0;
for(int j=n;j>=0;j--)
tmp=(tmpi+a[j])%p[o];
if(tmp)for(int j=0;i+jp[o]<=m;j++)
flag[i+jp[o]]=1;
}
}
int ans = 0;
for(int i=1;i<=m;i++)if(!flag[i])ans++;
cout<<ans<<endl;
for(int i=1;i<=m;i++)
if(!flag[i])
cout<<i<<endl;
}