Cod sursa(job #1909724)

Utilizator TudoseSanzianaTudose Sanziana TudoseSanziana Data 7 martie 2017 13:44:36
Problema Suma divizorilor Scor 30
Compilator cpp Status done
Runda Arhiva de probleme Marime 0.71 kb
#include <fstream>
using namespace std;

ifstream in("sumdiv.in");
ofstream out("sumdiv.out");
const int MOD=9901;
int a,b;
int ans=1;

int pow(int a,int b)
{
    int a1=a,ans;
    for(ans=1; b; b>>=1)
    {
        if(b&1)
            ans=(ans*a1)%MOD;
        a1=(a1*a1)%MOD;
    }
    return ans;

}
int main()
{
    int p;
    in>>a>>b;

    for(int i=2; i*i<=a; i++)
    {
        p=0;
        while(a%i==0)
        {
            p++;
            a/=i;
        }
        p*=b;
        if(p) ans*=((pow(i,p+1)-1)/(i-1));
        ans%=MOD;
    }

    if(a>1)
    {
        p=b;
        ans*=((pow(a,p+1)-1)/(a-1));
        ans%=MOD;
    }
    out<<ans%MOD;
    return 0;
}