Pagini recente » Cod sursa (job #438401) | Cod sursa (job #2251491) | Cod sursa (job #204830) | Istoria paginii runda/simulare_11_12/clasament | Cod sursa (job #539748)
Cod sursa(job #539748)
#include<stdio.h>
#define MOD 9901
int A,B,cnt,i,sol=1;
FILE *in,*out;
int putere(int a,int b)
{
int rez=1,n=a%MOD;
for(i=0;(1<<i)<=b;i++)
{
if((1<<i) & b)
rez=(rez*n)%MOD;
n=(n*n)%MOD;
}
return rez;
}
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))
{
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;
}