Cod sursa(job #3199725)

Utilizator Mihai_OctMihai Octavian Mihai_Oct Data 2 februarie 2024 16:14:26
Problema Suma divizorilor Scor 10
Compilator cpp-64 Status done
Runda Arhiva de probleme Marime 0.92 kb
#include <bits/stdc++.h>

using namespace std;

ifstream fin("sumdiv.in");
ofstream fout("sumdiv.out");
const int mod = 9901;
int n, b, d, e, p;
long long sum = 1;

static inline int Put(int a, int n) {
    long long p = 1;
    while(n) {
        if(n & 1) p = (p * a) % mod;
        a = (a * a) % mod;
        n >>= 1;
    }
    return p;
}

static inline int Inv(int a) {
    return Put(a, mod - 2);
}

int main() {
    fin >> n >> b;
    d = 2;
    while(d * d <= n) {
        int e = 0;
        while(n % d == 0) {
            e++;
            n /= d;
        }

        if(e) {
            int p = (Put(d, e * b + 1) - 1 + mod) % mod;
            sum = (sum * ((p - 1) / Inv(d - 1))) % mod;
        }

        d++;
    }
    if(n > 1) {
        int nn = (Put(n, 1 * b + 1) - 1 + mod) % mod;
        sum = (sum * ((nn - 1) / Inv(n - 1))) % mod;
    }

    fout << sum;

    return 0;
}