97. 约数之和

最新推荐文章于 2024-03-13 19:37:13 发布

追寻远方的人

最新推荐文章于 2024-03-13 19:37:13 发布

阅读量45

点赞数

分类专栏：算法提高复习文章标签：等比数列约数

本文链接：https://blog.csdn.net/segegse/article/details/130936880

版权

算法提高复习专栏收录该内容

167 篇文章 3 订阅

订阅专栏

假设现在有两个自然数 A 和 B，S 是 AB 的所有约数之和。

请你求出 Smod9901的值是多少。

输入格式

在一行中输入用空格隔开的两个整数 A 和 B。

输出格式

输出一个整数，代表 Smod9901 的值。

数据范围

0≤A,B≤5×107

输入样例：

2 3

输出样例：

注意: A 和 B 不会同时为 00。

$[外链图片转存失败,源站可能有防盗链机制,建议将图片保存下来直接上传(img-yKn4BiNs-1685366398827)(C:\Users\myp\AppData\Local\Temp\WeChat Files\1ff38ab283eba1105dbb64f10f9fdbe.png)]$

代码:

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

#define x first
#define y second
#define int long long
#define endl '\n'
#define IOS ios::sync_with_stdio(false); cin.tie(0);cout.tie(0);
typedef unsigned long long ull;
typedef pair<int, int> PII;

int lowbit(int x) { return x & (- x); }
int  gcd(int a, int b) { return b ? gcd(b, a % b) : a; }
const int mod = 9901;

int qmi(int a, int k)
{
    int res = 1;
    a = a % mod;

    while (k)
    {
        if (k & 1) res = res * a % mod;

        a = a * a % mod;
        k >>= 1;
    }

    return res;
}

int sum(int p, int k)
{
    if (k == 0) return 1;

    if (k % 2 == 0) return (1 + p % mod * sum(p, k - 1)) % mod;

    if (k % 2) return (1 + qmi(p, k / 2 + 1)) * sum(p, k / 2) % mod;
}


signed main()
{
    IOS;
    int a, b;
    cin >> a >> b;
    int res = 1;

    for (int i = 2; i <= a / i; i++)
    {
        if (a % i == 0)
        {
            int s = 0;

            while (a % i == 0)
            {
                a /= i, s++;
            }

            res = res * (sum(i, b * s)) % mod;
        }
    }

    if (a > 1) res = res * (sum(a, b)) % mod;

    if (a == 0) res = 0;

    cout << res;
    return 0;
}