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

来源：互联网发布：黄静案知乎编辑：程序博客网时间：2024/05/24 05:31

题目点这里

网上流行的解法是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;}

0 0