Pagini recente » Cod sursa (job #235224) | Cod sursa (job #547249) | Cod sursa (job #2565219) | Cod sursa (job #829819) | Cod sursa (job #1989067)
#include <iostream>
#include <cmath>
#include <fstream>
using namespace std;
ifstream f("sumdiv.in");
ofstream g("sumdiv.out");
int vp[78500];
long long x,nr,s,a,b,aa,bb;
long long Exp_Log_Mod(long long a,long long b,long long pr)
{
long long rez=1;
while(b)
{
if(b%2)
{
rez=(rez%pr)*(a%pr)%pr;
b--;
}
else
{
a=(a%pr)*(a%pr)%pr;
b/=2;
}
}
return rez;
}
long long Fractie(long long d,long long p,long long pr)
{
long long nmrt,numi;
nmrt=Exp_Log_Mod(d,p,pr)-1;
numi=Exp_Log_Mod(d-1,pr-2,pr);
return (nmrt*numi)%pr;
}
void Nr_Sum_Div(int x,int bb)
{
int r,p,i=1;
s=1;
r=sqrt(x);
while(vp[i]<=r)
{
p=0;
while(x%vp[i]==0)
{
p++;
x/=vp[i];
}
if(p)
{
s=s*Fractie(vp[i],p*bb+1,9901)%9901;
r=sqrt(x);
}
i++;
}
if(x!=1)
s*=Fractie(x,bb+1,9901)%9901;;
}
void Ciur()
{
int n=0;
int i,j;
bool c[1000000];
for(i=2;i<=1000000;i++)
if(!c[i])
{
vp[++n]=i;
for(j=2*i;j<=1000000;j+=i)
c[j]=1;
}
}
int main()
{
int aa,bb;
Ciur();
f>>aa>>bb;
Nr_Sum_Div(aa,bb);
g<<s;
f.close();
g.close();
return 0;
}