Cod sursa(job #2663922)

Utilizator metallidethantralayerIon Cojocaru metallidethantralayer Data 27 octombrie 2020 16:44:51
Problema Suma divizorilor Scor 30
Compilator cpp-64 Status done
Runda Arhiva de probleme Marime 0.84 kb
#include <bits/stdc++.h>
#define MOD 9901
using namespace std;
ifstream f("sumdiv.in");
ofstream g("sumdiv.out");
long long int n,m;
long long int PutereExp(long long int a,long long int b)
{
    long long int rez=1;
    while(b)
    {
        if(b&1)
            rez=(rez*a)%MOD;
        a=(a*a)%MOD;
        b>>=1;
    }
    return rez;
}
long long int CalculFactorPrim(long long int x)
{
    return (1LL*((PutereExp(x,m+1)-1+MOD)%MOD)*PutereExp(x-1,MOD-2))%MOD;
}
int main()
{
    f>>n>>m;
    long long int d=2,rez=1;
    while(n>1)
    {
        bool ok=true;
        if(!(n%d))
        {
            ok=false;
            while(!(n%d))
                n/=d;
        }
        if(!ok)
            rez*=CalculFactorPrim(d);
        d++;
        if(d*d>n&&n>1)
            d=n;

    }
    g<<rez;
    return 0;
}