Pagini recente » Cod sursa (job #845390) | Cod sursa (job #1075216) | Cod sursa (job #939002) | Cod sursa (job #2528532) | Cod sursa (job #2766507)
#include <fstream>
#include <iostream>
using namespace std;
const int MOD = 9901;
int a, b;
void read() {
ifstream f("sumdiv.in");
f >> a >> b;
f.close();
}
int S = 1;
int exp(int a, int b) {
int P = 1;
a %= MOD;
while (b) {
if (b & 1)
P = (1LL * P * a) % MOD;
a = (1LL * a * a) % MOD;
b /= 2;
}
return P;
}
void solve() {
int d = 2, p;
while (a != 1) {
p = 0;
while (a % d == 0) {
p++;
a /= d;
}
if (p != 0)
S = (1LL * S * (exp(d, p * b + 1) - 1 + MOD) % MOD * exp(d - 1, MOD - 2) % MOD) % MOD;
d++;
if (d * d > a && a != 1) {
S = (1LL * S * (exp(a, b + 1) - 1 + MOD) % MOD * exp(a - 1, MOD - 2) % MOD) % MOD;
break;
}
}
}
void output() {
ofstream g("sumdiv.out");
g << (a == 0 ? 0 : S);
g.close();
}
int main() {
read();
solve();
output();
return 0;
}