Pagini recente » Cod sursa (job #1969243) | Cod sursa (job #1975177) | Cod sursa (job #2653023) | Cod sursa (job #2767021) | Cod sursa (job #3002265)
#include <fstream>
using namespace std;
const int MOD = 9901;
ifstream f("sumdiv.in");
ofstream g("sumdiv.out");
int powlg(int x, int p)
{
int v = 1, a = x % MOD;
while(p)
{
if(p & 1)
v = v * a % MOD;
a = a * a % MOD;
p >>= 1;
}
return v;
}
inline int invMod(int x)
{
return powlg(x, MOD - 2);
}
int sFact(int d, int exp)
{
int x = d % MOD;
if(x == 0)
return 1;
if(x == 1)
return (exp + 1) % MOD;
return (powlg(d, exp + 1) + MOD - 1) % MOD * invMod(d - 1) % MOD;
}
int main()
{
int A, B, S = 1, p;
f >> A >> B;
for(int d = 2; d * d <= A; d++)
if(A % d == 0)
{
p = 0;
do
{
p++;
A /= d;
}
while(A % d == 0);
S = S * sFact(d, B * p) % MOD;
}
if(A > 1) S = S * sFact(A, B) % MOD;
g << S;
f.close();
g.close();
return 0;
}