noip2014解方程

3．解方程
(equation.cpp/c/pas)
【问题描述】
已知多项式方程：
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]内的整数解的个数。

#include<cstdio>
#include<cstring>
#include<algorithm>
using namespace std;
const int maxn=10010;
char st[maxn]; 
int mod[5]={9973,9931,9941,9949,9967};
int pos[5][maxn];
int n,m;
int ans[maxn];
int res[5][maxn];
int calc(int t,int x){
    int sum=0;
    for(int i=n;i>=0;i--){
        sum=((sum*x)%mod[t]+pos[t][i])%mod[t];
    }
    return sum;
}
bool jud(int x){
    for(int i=0;i<5;i++){
        if(res[i][x%mod[i]]!=0) return false;
    }
    return true;
}
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++){
        scanf("%s\n",st);
        int l=strlen(st);
        int flag=0;
        if(st[0]!='-') for(int k=0;k<5;k++) pos[k][i]=st[0]-'0';
        else flag=1;
        for(int j=1;j<l;j++){
            for(int k=0;k<5;k++){
                pos[k][i]=(pos[k][i]*10+st[j]-'0')%mod[k];
            }
        }
        if(flag==1)
        for(int k=0;k<5;k++) pos[k][i]*=-1;
       /* for(int i=0;i<=n;i++) printf("%d ",pos[0][i]);
        printf("\n");*/
    }
    for(int x=1;x<min(mod[0],m+1);x++){
        for(int t=0;t<5;t++){
            res[t][x]=calc(t,x);
        }
    }
    //printf("%d",calc(0,1));
    int pus=0;
    for(int i=1;i<=m;i++){
        if(jud(i)) ans[++pus]=i; 
    }
    printf("%d\n",pus);
    for(int i=1;i<=pus;i++){
        printf("%d\n",ans[i]);
    }
    return 0;
}