Pagini recente » Cod sursa (job #3326178) | Cod sursa (job #3327396) | Cod sursa (job #2309683) | Cod sursa (job #2653063) | Cod sursa (job #3322741)
#include <bits/stdc++.h>
using namespace std;
const int MOD = 9901;
int invMod( int b ) {
int inv = 1;
while ( inv < MOD && (b * inv) % MOD != 1 )
inv++;
return inv;
}
int pow_mod( int p, int e ) {
int rez = 1, a = p;
while ( e > 0 ) {
if ( e % 2 == 1 )
rez = ((long long)rez * a) % MOD;
a = ((long long)a * a) % MOD;
e /= 2;
}
return rez;
}
int main() {
FILE *fin, *fout;
int a, b, sdiv, i, e, imp;
fin = fopen( "sumdiv.in", "r" );
fscanf( fin, "%d%d", &a, &b );
fclose( fin );
sdiv = imp = 1;
for ( i = 2; i * i <= a; i++ ) {
if ( a % i == 0 ) {
e = 0;
while ( a % i == 0 ) {
e++;
a /= i;
}
sdiv = (sdiv * (pow_mod( i % MOD, e * b + 1 ) - 1)) % MOD;
imp = (imp * (i - 1)) % MOD;
}
}
if ( a > 1 ) {
sdiv = (sdiv * (pow_mod( a % MOD, b + 1 ) - 1)) % MOD;
imp = (imp * (a - 1)) % MOD;
}
fout = fopen( "sumdiv.out", "w" );
fprintf( fout, "%d\n", (sdiv * invMod( imp )) % MOD );
fclose( fout );
return 0;
}