Cod sursa(job #539920)

Utilizator SadmannCornigeanu Calin Sadmann Data 23 februarie 2011 14:56:40
Problema Dezastru Scor 0
Compilator cpp Status done
Runda Arhiva de probleme Marime 0.96 kb
#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;
}