CODEVS 3732 解方程 NOIP 2014

题目大意

已知多项式方程：
a0+a1x+a2x^2+…+anx^n=0
求这个方程在[1, m]内的整数解（n 和 m 均为正整数）

• 输入描述

  输入共 n+2 行。
  第一行包含 2 个整数 n、m，每两个整数之间用一个空格隔开
  接下来的 n+1 行每行包含一个整数，依次为a0,a1,a2,…,an

• 输出描述

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

• 样例输入

  2 10
  1
  -2
  1

• 样例输出

  1
  1

• 数据范围

  对于 30%的数据，0 < n ≤ 2， |ai| ≤ 100，an ≠ 0， m ≤ 100；
  对于 50%的数据，0 < n ≤ 100， |ai| ≤ 10^100 ，an ≠ 0，m ≤ 100；
  对于 70%的数据，0 < n ≤ 100， |ai| ≤ 10^10000 ，an ≠ 0，m ≤ 10000；
  对于 100%的数据，0 < n ≤ 100， |ai| ≤ 10^10000 ，an ≠ 0，m ≤ 1000000。

分析

很少做这样的数学题，看了别人的题解A的。
高精度的话，还是算了（ＴＡＴ）
用神奇的模运算吧！把整个方程取分别取几个模，再枚举答案就可以了。如果枚举的ｘ都满足这几个方程，那么ｘ就是一个根。
再一个很神奇的就是，我用的那几个MOD值，提交后居然WA了！然后再用被人题解里的MOD值，就Ａ了！
（一脸懵逼

#include <cstdio>
#include <cstdlib>
#include <cstring>
#include <algorithm>
using namespace std;
const int N = 110, M = 1000010;
int MOD[5] = {11261,19997,22877,21893,14843};   //这是我WA掉的MOD值　30011,40009,49999,59999,70001
int pre[5][N], a[5][N], res[5][80000];
int ans[M];
int n, m;
char st[10010];
int readin(char st[], int mod) {

    int sum = 0;
    int flag = 1;
    int l = strlen(st);
    if(st[0] == '-') flag = -1;
    else sum += st[0] - '0';
    for(int i = 1; i < l; i++) {
        sum = ( sum * 10 + (st[i] - '0') ) % mod;
    }
    return sum*flag;        
}

int cal(int x, int mod, int j) {
    int sum = 0;
    for(int i = 0; i <= n; i++) {
        sum = (sum + pre[j][i]*a[j][i]) % mod;
    }
    if(sum < 0) sum += mod;
    return sum;
}

bool judge(int x) {
    for(int j = 0; j < 5; j++) {
        if(res[j][x % MOD[j]] != 0) return 0;
    }
    return 1;
}

int main() {
    scanf("%d%d", &n, &m);
    for(int i = 0; i <= n; i++) {
        scanf("%s", st);
        for(int j = 0; j < 5; j++) {
            a[j][i] = readin(st, MOD[j]);
        }
    }
    for(int j = 0; j < 5; j++) {
        pre[j][0] = 1;
        for(int x = 1; x < MOD[j]; x++) {
            for(int k = 1; k <= n; k++) {
                pre[j][k] = ( pre[j][k-1] * x ) % MOD[j];
            }
            res[j][x] = cal(x, MOD[j], j);
        }
    }
    for(int x = 1; x <= m; x++) {
        if(judge(x)) ans[++ans[0]] = x;
    }
    printf("%d\n", ans[0]);
    for(int i = 1; i <= ans[0]; i++) {
        printf("%d\n", ans[i]);
    }
    return 0;
}