Pagini recente » Cod sursa (job #2062227) | Cod sursa (job #303259) | Cod sursa (job #2175045) | Cod sursa (job #1501497) | Cod sursa (job #434313)
Cod sursa(job #434313)
#include <cstdio>
const int MOD = 9901;
long long A, B;
long long pw(long long a, long long p)
{
long long rez = 1;
a = a % MOD;
for(; p ; p = p>>1)
{
if(p & 1)
rez = (rez * a) % MOD , p--;
a = (a * a) % MOD;
}
return rez;
}
void rezolva()
{
long long n = A;
long long S = 1;
for(long long i = 2 ; i * i <= A ; i++)
if(n % i == 0)
{
int exp = 0;
while(n % i == 0)
{
n /= i;
exp++;
}
S = ((S * pw(i, B * exp + 1) + MOD - 1) % MOD * pw(i - 1, MOD - 2)) % MOD;
}
if(n > 1)
S = ((S * (pw(n, B + 1) + MOD - 1) % MOD) * pw(n - 1, MOD - 2)) % MOD;
printf("%lld", S);
}
int main()
{
freopen("sumdiv.in","r",stdin);
freopen("sumdiv.out","w",stdout);
scanf("%lld%lld",&A,&B);
rezolva();
return 0;
}