Cod sursa(job #1060861)

Utilizator auRSTARHreapca Aurelian auRSTAR Data 18 decembrie 2013 20:43:32
Problema Suma divizorilor Scor 0
Compilator cpp Status done
Runda Arhiva de probleme Marime 0.84 kb
#include<cstdio>
#define LL long long
#define MOD 9901
LL A,B,result,i,exp;
LL pow(LL x,LL y)
{
    int res=1;
    for(;y;y>>=1)
    {
        if(y&1)
            res=(res*x)%MOD;
        x=(x*x)%MOD;
    }
    return res;
}
int main()
{
    freopen("sumdiv.in","r",stdin);
    //freopen("sumdiv.out","w",stdout);
    scanf("%lld%lld",&A,&B);
    result=1;
    for(i=2;i*i<=A;i++)
    {
        for(exp=0;!(A%i);exp++,A/=i);
        if(exp)
        {
            LL g=(pow(i,exp*B+1)-1)%MOD;
            if(g<0)g+=MOD;
            result=(result*g)%MOD;
            result=(result*pow(i-1,MOD-2))%MOD;
        }
    }
    if(A>1)
    {
        LL g=(pow(A,B+1)-1)%MOD;
        if(g<0)g+=MOD;
        result=(result*g)%MOD;
        result=(result*pow(A-1,MOD-2))%MOD;
    }
    printf("%lld\n",result);
    return 0;
}