【NOIP2014】解方程

题目链接：https://www.luogu.org/problemnew/show/P2312

//这道题一看数据范围就应该向取模的方向想。

30分做法，直接手推，带入即可。

50分听说可以高精，没搞。

100分一个很明显的思路就是[1,m]内枚举x，然后代入方程，看一下是否为零。但这样有两个问题：

(1)怎么搞这个系数。前面也说过，可以取模，那么我们的根据是什么呢？f(x)=0=>f(x)%p=0，就是这个咯。详细的做法需要用到读入优化，一边读一边模即可。

(2)时间复杂度太高：这里用的是秦九韶算法，无论是手动计算还是计算机计算都能够大幅度减少复杂度的算法，主要用来解一元n次多项式的问题，简略说下就是：一般地，一元n次多项式的求值需要经过(n+1)*n/2次乘法和n次加法，而秦九韶算法只需要n次乘法和n次加法。详细见：https://baike.baidu.com/item/秦九韶算法/449196?fr=aladdin

最后在说下模数的选择，大模数以及质数即可，可以多试几次。

code：

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

typedef long long ll; 
const int N=1e2+10, M=1e6+10, mod=1e9+7;
int n, m, cnt, ans[N];
ll a[N];

ll read() {
    ll x=0, w=1;
    char ch=0;
    while (ch<'0' || ch>'9') {
        if (ch == '-') w=-1;
        ch=getchar();
    }
    while (ch>='0' && ch<='9') {
        x=((x<<3)+(x<<1)+ch-'0')%mod;
        ch=getchar();
    }
    return x*w;
}

int main() {
    cin >> n >> m;
    for (int i=0; i<=n; i++)
        a[i]=read();
    int x=1;
    while (x<=m) {
        ll tmp=(a[n]*x+a[n-1])%mod;
        for (int i=n-2; i>=0; i--)
            tmp=((tmp*x)%mod+a[i])%mod;
        if (tmp==0) ans[++cnt]=x;
        x++;
    }
    cout << cnt << endl;
    for (int i=1; i<=cnt; i++)
        cout << ans[i] << endl;
    return 0;
}