NOIP2014 解方程

题目描述

已知多项式方程：

a0+a1x+a2x^2+..+anx^n=0

求这个方程在[1, m ] 内的整数解（n 和m 均为正整数）

输入输出格式

输入格式：
输入文件名为equation .in。
输入共n + 2 行。
第一行包含2 个整数n 、m ，每两个整数之间用一个空格隔开。
接下来的n+1 行每行包含一个整数，依次为a0,a1,a2..an
输出格式：
输出文件名为equation .out 。
第一行输出方程在[1, m ] 内的整数解的个数。
接下来每行一个整数，按照从小到大的顺序依次输出方程在[1, m ] 内的一个整数解。

输入输出样例

输入样例#1：
2 10
1
-2
1
输出样例#1：
1
1
输入样例#2：
2 10
2
-3
1
输出样例#2：
2
1
2
输入样例#3：
2 10
1
3
2

输出样例#3：
0

代码

#include<cstdio>#include<cstring>#define M 1000001#include<algorithm>using namespace std;int p[6]={0,7457,7829,8623,8971,9781};int pre[6][105],a[6][105],ans[M],res[6][30005];char c[10000];int n,m;inline int read() {    int x=0,f=1;char c=getchar();    while(c>'9'||c<'0') {if(c=='-')f=-1;c=getchar();}     while(c>='0'&&c<='9') x=10*x+c-48,c=getchar();    return x*f;}inline int cal(int t) {    int sum=0;    for(int i=0;i<=n;i++) {        sum+=(a[t][i]*pre[t][i])%p[t];        sum%=p[t];    }    if(sum<0) sum+=p[t];    return sum;}inline bool judge(int x) {    for(int i=1;i<=5;i++)      if(res[i][x%p[i]]) return false;    return true;}int main() {    bool flag=false;    n=read();m=read();    for(int i=0;i<=n;i++) {        flag=false;        scanf("%s",c+1);        int l=strlen(c+1);        for(int t=1;t<=5;t++)          if(c[1]!='-') a[t][i]=c[1]-48;            else a[t][i]=0,flag=true;        for(int t=1;t<=5;t++) {            for(int k=2;k<=l;k++)                a[t][i]=(a[t][i]*10+c[k]-48)%p[t];            if(flag) a[t][i]=-a[t][i];        }    }    for(int j=1;j<=5;j++)       for(int k=1;k<p[j];k++) {        pre[j][0]=1;        for(int i=1;i<=n;i++)          pre[j][i]=(pre[j][i-1]*k)%p[j];        res[j][k]=cal(j);      }    for(int i=1;i<=m;i++)      if(judge(i)) ans[++ans[0]]=i;    printf("%d\n",ans[0]);    for(int i=1;i<=ans[0];i++)      printf("%d\n",ans[i]);    return 0;}