AcWing - 515 - 解方程 = 简单数学

标签：write   top   using   namespace   else   include   而不是   https   整数   

https://www.acwing.com/problem/content/517/

给一个多项式，求他的整数根。

首先稳妥的办法应该是整一大堆质数，然后用中国剩余定理合并（当然不是真的合并）。

奇怪一点的就用几个就可以了，卡掉的概率极低，加上是OI其实没问题的。

但是这题卡常，要用秦九韶公式卡掉一半的乘法（和取模），还卡读入，真的垃圾题。

注意读入的时候要%p而不是-p，找到一个写法的bug不错。

#include <bits/stdc++.h>
using namespace std;
typedef long long ll;

const int MAXP = 1;
const int p = 19260817;

int ap[105];

inline int read() {
    int b = 0;
    char c = getchar();
    while(true) {
        if(c == '-' || (c >= '0') && (c <= '9'))
            break;
        c = getchar();
    }
    if(c == '-') {
        b = 1;
        c = getchar();
    }
    int cur = 0;
    while(c >= '0' && c <= '9') {
        cur = (cur << 3) + (cur << 1) + (c - '0');
        if(cur >= p)
            cur %= p;
        c = getchar();
    }
    if(b)
        cur = (p - cur) % p;
    return cur;
}

inline void write(int x) {
    if(x <= 9) {
        putchar('0' + x);
    } else {
        write(x / 10);
        putchar('0' + x % 10);
    }
}

int ans[1000005];
int atop = 0;

int main() {
#ifdef Yinku
    freopen("Yinku.in", "r", stdin);
#endif // Yinku
    int n = read(), m = read();
    for(int i = 0; i <= n; ++i)
        ap[i] = read();

    for(int x = 1; x <= m; ++x) {
        ll cur = 0;
        for(int i = n; i >= 1; --i) {
            cur = (cur + ap[i]) * x;
            if(cur >= p)
                cur %= p;
        }
        cur = (cur + ap[0]) % p;
        if(cur == 0)
            ans[++atop] = x;
    }

    write(atop);
    putchar('\n');
    for(int i = 1; i <= atop; ++i) {
        write(ans[i]);
        putchar('\n');
    }
    return 0;
}

AcWing - 515 - 解方程 = 简单数学

标签：write   top   using   namespace   else   include   而不是   https   整数   

原文地址：https://www.cnblogs.com/Inko/p/11664463.html