Pagini recente » Cod sursa (job #1894913) | Cod sursa (job #1069634) | Cod sursa (job #784771) | Cod sursa (job #834516) | Cod sursa (job #756933)
Cod sursa(job #756933)
#include <fstream>
using namespace std;
int put (int a, int b, int mod) {
if (a == 1) return 1 % mod;
if (b == 0) return 1 % mod;
if (b == 1) return a % mod;
int temp = put (a, b >> 1, mod);
temp = (temp * temp) % mod;
if (b & 1) temp = (temp * a) % mod;
return temp;
}
int invers (int val, int prim) {
return put (val, prim - 2, prim);
}
int main () {
int a, b, t, p, s = 1, mod = 9901;
ifstream f ("sumdiv.in");
ofstream g ("sumdiv.out");
f >> a >> b;
for (p = 2; p * p <= a; ++ p) {
t = 1;
while (a % p == 0) {
a /= p;
t *= p;
}
if (t > 1) {
t = t % mod;
t = (t * put (t, b, mod)) % mod;
-- t;
if (t < 0) t += mod;
s = (s * t) % mod;
s = (s * invers (p - 1, mod)) % mod;
}
}
if (a > 1) {
t = a;
t = t % mod;
t = (t * put (t, b, mod)) % mod;
-- t;
if (t < 0) t += mod;
s = (s * t) % mod;
s = (s * invers (a - 1, mod)) % mod;
}
g << s << endl;
return 0;
}