【Ybtoj 第29章例2】约数之和【同余问题】

最新推荐文章于 2022-02-10 14:06:14 发布

。十二。

最新推荐文章于 2022-02-10 14:06:14 发布

阅读量182

点赞数

分类专栏： Ybtoj 数学文章标签：算法哈希算法深度优先

本文链接：https://blog.csdn.net/kejin2019/article/details/119786741

版权

Ybtoj 同时被 2 个专栏收录

142 篇文章 3 订阅

订阅专栏

数学

119 篇文章 0 订阅

订阅专栏

这篇博客主要介绍了如何利用模指数运算（快速幂）和素数分解来解决大整数乘法问题。作者首先讲解了如何通过素数分解将大整数A表示为素数的乘积，然后使用模指数运算求解B的A次方模素数的逆元，从而高效地进行计算。代码示例中展示了如何实现这些算法，并在主函数中展示了整个计算过程。

摘要由CSDN通过智能技术生成

在这里插入图片描述

解题思路

在这里插入图片描述

代码

#include <iostream>
#include <cmath>
#include <cstring>
#include <cstdio>
#include <algorithm>
#define ll long long

using namespace std;

const int mod = 9901;
ll A, B, cnt, ans = 1ll, pri[20010], v[200000010], c[2010];

void prime(int x) {
    for (int i = 2; i * i <= x; i++) {
        if (x % i == 0)
            pri[++cnt] = i;
        while (x % i == 0) {
            c[cnt]++;
            x /= i;
        }
    }
    if (x > 1) {
        pri[++cnt] = x;
        c[cnt] = 1;
    }
}

ll ksm(ll a, ll x) {
    ll ans = 1, base = a;
    while (x) {
        if (x & 1)
            ans = ans * base % mod;
        base = base * base % mod;
        x = x >> 1;
    }
    return ans;
}

int main() {
    scanf("%lld%lld", &A, &B);
    prime(A);
    for (int i = 1; i <= cnt; i++) {
        if ((pri[i] - 1) % mod == 0) {
            ans = ans * (B * c[i] + 1) % mod;
            continue;
        }
        int w = ksm(pri[i], B * c[i] + 1);
        w = (w + mod - 1) % mod;
        int v = pri[i] - 1;
        v = ksm(v, mod - 2);
        ans = ans % mod * w * v % mod;
    }
    printf("%lld", ans);
}