【NOIP】提高组2014 解方程

【题意】已知n次方程（n<=100）及其所有系数（|ai|<=10^10000），求[1,m]中整数解的个数（m<=10^6）。

【算法】数论

【题解】如果f(x)=0，则有f(x)%p=0。

所以取若干个素数p，将所有数字读入取模并快速计算出所有f(x)%p，若均为0则认为f(x)=0。

优化：利用f(x)%p=f(x%p)，可以将枚举范围缩小。

#include<cstdio>
#include<cstring>
#include<algorithm>
using namespace std;
const int p[]={,,,,,,,,,,};
int n,m,a[][],f[];
bool b[];
char s[];
int power(int x,int k,int mod){
    int ans=;
    while(k){
        if(k&)ans=ans*x%mod;
        x=x*x%mod;
        k>>=;
    }
    return ans;
}
int main(){
    scanf("%d%d",&n,&m);
    for(int i=;i<=n;i++){
        scanf("%s",s+);
        int len=strlen(s+);
        for(int j=(s[]=='-'?:);j<=len;j++){
            for(int k=;k<=;k++)a[k][i]=(a[k][i]*+(s[j]-''))%p[k];
        }
        if(s[]=='-')for(int k=;k<=;k++)a[k][i]=-a[k][i];
    }
    for(int i=;i<=m;i++)b[i]=;
    for(int k=;k<=;k++){
        for(int i=;i<p[k];i++){//nai xin man man kan
            f[i]=;
            for(int j=;j<=n;j++){
                f[i]=(f[i]+a[k][j]*power(i,j,p[k]))%p[k];//
            }
            if(f[i]!=)b[i]=;
        }
        for(int i=p[k];i<=m;i++){
            if(f[i%p[k]]!=)b[i]=;
        }
    }
    int ans=;
    for(int i=;i<=m;i++)if(b[i])ans++;
    printf("%d\n",ans);
    for(int i=;i<=m;i++)if(b[i])printf("%d\n",i);
    return ;
}