【数论】[NOIP2014]bzoj3751 codevs3732 解方程

最新推荐文章于 2020-10-07 16:14:14 发布

MoeO3

最新推荐文章于 2020-10-07 16:14:14 发布

阅读量892

点赞数

分类专栏：数学文章标签： noip2014 数论

本文链接：https://blog.csdn.net/qq_21841245/article/details/48595047

版权

数学专栏收录该内容

8 篇文章 0 订阅

订阅专栏

题目点这里

网上流行的解法是mod五个小素数，因为十个数据点嘛除非人品太好了不然一般都不会炸吧……

然后之前vfk讲了万无一失的做法就是先用一个小素数筛，再用一个大素数筛，然后高精度带回去验证……

本来打算把完整解法写完的不过没力气写高精度了……想当初我是现场写高精度的啊……

选了一个颜值比较高的素数和UOJ素数然后就水过了…………那就这样吧………………作为一只高三党狗还是好好刷常规吧TAT

为了防止炸int改了半天……秦九韶算法还去查了公式不然真的不知道怎么玩了……

#include <cstdio>
#include <iostream>

using namespace std;

int n, m;
int p[] = {38833, 998244353};
int mod[2][105], ans[1000005], first[1000005], second[1000005];

char s[10005];
#define getn(i) (i - '0')

int main()
{
    cin >> n >> m; bool flag;
    for (int i = 0; i <= n; ++ i) {
        scanf("%s", s + 1);
        int &a = mod[0][i], &b = mod[1][i];
        if (s[1] == '-') flag = 1;
        else { flag = 0; a = b = getn(s[1]);}
        for (int j = 2; s[j]; ++ j) {
            a = a * 10 + getn(s[j]); if (a > p[0]) a %= p[0];
            long long t = (long long)b * 10 + getn(s[j]); b = (t > p[1]) ? t % p[1] : t;
        }
        if (flag) { a = -a; b = -b; }
    }
    // 第一次筛选
    for (int x = 1; x < p[0]; ++ x) {
        ans[x] = mod[0][n];
        for (int i = n - 1; ~i; -- i) ans[x] = (ans[x] * x + mod[0][i]) % p[0];
    }
    for (int i = 1; i <= m; ++ i) if (!ans[i % p[0]]) first[++first[0]] = i;
    // 第二次筛选
    for (int i = 1; i <= first[0]; ++ i) {
        int x = first[i]; ans[i] = mod[1][n];
        for (int j = n - 1; ~j; -- j) ans[i] = ((long long)ans[i] * x + mod[1][j]) % p[1];
    }
    for (int i = 1; i <= first[0]; ++ i) if (!ans[i]) second[++second[0]] = first[i];
    //此处应有高精度验证……
    cout << second[0] << endl;
    for (int i = 1; i <= second[0]; ++ i) printf("%d\n", second[i]);
    return 0;
}