Pagini recente » Cod sursa (job #2805615) | Cod sursa (job #1630915) | Cod sursa (job #2382320) | Cod sursa (job #945398) | Cod sursa (job #2468806)
#include <fstream>
#define MOD 9901
using namespace std;
ifstream fin("sumdiv.in");
ofstream fout("sumdiv.out");
int a,b;
int lgput(int x, int y)
{
if (y == 1)
return x;
int val = lgput(x, y/2)%MOD;
if (y%2 == 0)
return (val*val)%MOD;
else
return (((val*val)%MOD)*x)%MOD;
}
int main()
{
fin >> a >> b; int sol = 1;
for (int d=2; d<=a/d; d++)
if (a%d == 0)
{
int exponent = 0;
while (a%d == 0)
{
a /= d;
exponent++;
}
sol *= lgput(d, exponent*b+1)-1;
sol *= lgput(d-1, MOD-2); sol %= MOD;
}
if (a != 1)
{
sol *= lgput(a, b+1)-1; sol %= MOD;
sol *= lgput(a-1, MOD-2); sol %= MOD;
}
fout << sol;
return 0;
}