Pagini recente » Cod sursa (job #956096) | Cod sursa (job #3338101) | Cod sursa (job #2334414) | Cod sursa (job #1544543) | Cod sursa (job #3322745)
#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;
while ( e > 0 ) {
if ( e % 2 == 1 )
rez = ((long long)rez * p) % MOD;
p = ((long long)p * p) % 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, e * b + 1 ) - 1)) % MOD;
imp = (imp * (i - 1)) % MOD;
}
}
if ( a > 1 ) {
sdiv = (sdiv * (pow_mod( a, 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;
}