[NOIP 2014]解方程 秦九韶算法+质数取模

代码

#include <cstdio>
#include <iostream>
#include <cstring>
using namespace std;

int n,m,len,prime[5]={10007,11261,14843,19997,21893},sum(1);
long long a[105][10],f[30000][10],ans[1001001];

long long Calculate(int x,int j){
    long long t=0;
    for(int i=n;i>=0;i--){
        t=(t*x+a[i][j])%prime[j];
    }
    return t;
}

void Read(int x) {
    static char s[10100];int flag=0;
    scanf("%s",s);
    len=strlen(s);
    for(int i=0;i<len;i++){
        if(s[i]=='-'){
            flag=1;
        }else{
            for(int j=0;j<5;j++){
                a[x][j]=(a[x][j]*10+s[i]-'0')%prime[j];
            }           
        }

    }
    if(flag==1){
        for(int j=0;j<5;j++){
            a[x][j]=prime[j]-a[x][j];
        }
    }
}

int main() {
    scanf("%d%d",&n,&m);
    for(int i=0;i<=n;i++) {
        Read(i);
    }
    for(int j=0;j<5;j++){
        for(int i=0;i<=prime[j];i++){
            f[i][j]=Calculate(i,j);
        }
    }   
    for(int i=1;i<=m;i++){
        int j=0;
        for(j=0;j<5;j++){
            if(f[i%prime[j]][j]){
                break;
            }
        }
        if(j==5){
            ans[sum++]=i;
        }
    }
    printf("%d\n",sum-1);
    for(int i=1;i<=sum-1;i++){
        printf("%lld\n",ans[i]);
    }
    return 0;
}