Pagini recente » Cod sursa (job #2899765) | Cod sursa (job #479980) | Cod sursa (job #1761497) | Istoria paginii utilizator/gugustiucfan | Cod sursa (job #539920)
Cod sursa(job #539920)
#include<stdio.h>
#define MOD 9901
int cnt,sol=1,i;
int a,b,n;
FILE *in,*out;
int putere(int a,int b)
{
int sol=1;
n=a%MOD;
for(int q=0;(1<<q)<=b;q++)
{
if((1<<q) & b)
sol*=n , sol%=MOD;
n=(n*n)%MOD;
}
return sol;
}
int sum(int a,int b)
{
if(b==1)
return a%MOD;
if(b & 1)
return (sum(a,b-1)+putere(a,b))%MOD;
return ((putere(a,b/2)+1)*sum(a,b/2))%MOD;
}
int main()
{
in=fopen("sumdiv.in","rt");
out=fopen("sumdiv.out","wt");
fscanf(in,"%d %d",&a,&b);
if(!b)
{
fprintf(out,"1\n");
return 0;
}
for(i=2;i*i<=a;i++)
{
if(!(a%i))
{
for(;a%i==0;)
{
a/=i;
cnt++;
}
sol=sol*(sum(i,cnt*b)+1)%MOD;
}
}
if(a!=1)
sol=sol*(sum(a,b)+1)%MOD;
fprintf(out,"%d",sol);
return 0;
}