NOIP 2014 解方程

题目描述 Description

输入描述 Input Description

输入文件名为equation.in。

输入共n+2行。

第一行包含2个整数n、m，每两个整数之间用一个空格隔开。

接下来的n+1行每行包含一个整数，依次为a0,a1,a2,……,an。

输出描述 Output Description

输出文件名为equation.out。

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

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

样例输入 Sample Input

equation.in	equation.out
2 10 1 -2 1	1 1

equation.in	equation.out
2 10 2  -3 1	2 1 2

样例输出 Sample Output

equation.in	equation.out
2 10 1 3 2	0

数据范围及提示 Data Size & Hint

http://codevs.cn/problem/3732/

30% 我们只需要枚举1到m并且计算f(i)是否为0即可

50% 与30%思路相同，只是改为高精度（建议使用重载运算符）

70% 取5，6个较大的模数，将ai输入的同时进行取模，如果都为0就记为答案

100% 对于f(i)的值取模p时,f(i)=f(i%p)所以不需要枚举到m,只要到每个质数-1,即可，而较大的数的值只需找对应的5个值查看即可。

#include<cstdio>
#include<cstring>
int n,pr[5]={11261,19997,22877,21893,14843},qm[105][5],jl[1000005];
bool vis[5][30005];
char c[10005];
int f(int a,int b,int c)
{
    if(a==n) return qm[n][b];
    return (qm[a][b]+f(a+1,b,c)*c%pr[b])%pr[b];
}
int main()
{
    int m,ans=0;
    scanf("%d%d",&n,&m);
    for(int i=0;i<=n;++i)
    {
        scanf("%s",c);
        int l=strlen(c),f=c[0]=='-'?1:0;
        for(int j=0;j<5;++j)
        {
            for(int k=f;k<l;++k)
                qm[i][j]=(qm[i][j]*10+c[k]-'0')%pr[j];
            if(f) qm[i][j]=-qm[i][j];
        }
    }
    for(int i=0;i<5;++i)
        for(int j=0;j<pr[i];++j)
            if(f(0,i,j)) vis[i][j]=1;
    for(int i=1;i<=m;++i)
    {
        bool sf=1;
        for(int j=0;j<5;++j)
            if(vis[j][i%pr[j]]) sf=0;
        if(sf) jl[ans++]=i;
    }
    printf("%d
",ans);
    for(int i=0;i<ans;++i) printf("%d
",jl[i]);
    return 0;
}