Cod sursa(job #2468806)

Utilizator LeVladzCiuperceanu Vlad LeVladz Data 5 octombrie 2019 23:05:44
Problema Suma divizorilor Scor 40
Compilator cpp-64 Status done
Runda Arhiva de probleme Marime 0.85 kb
#include <fstream>
#define MOD 9901

using namespace std;

ifstream fin("sumdiv.in");
ofstream fout("sumdiv.out");

int a,b;

int lgput(int x, int y)
{
    if (y == 1)
        return x;
    int val = lgput(x, y/2)%MOD;
    if (y%2 == 0)
        return (val*val)%MOD;
    else
        return (((val*val)%MOD)*x)%MOD;
}

int main()
{
    fin >> a >> b; int sol = 1;
    for (int d=2; d<=a/d; d++)
        if (a%d == 0)
        {
            int exponent = 0;
            while (a%d == 0)
            {
                a /= d;
                exponent++;
            }
            sol *= lgput(d, exponent*b+1)-1;
            sol *= lgput(d-1, MOD-2); sol %= MOD;
        }
    if (a != 1)
    {
        sol *= lgput(a, b+1)-1; sol %= MOD;
        sol *= lgput(a-1, MOD-2); sol %= MOD;
    }
    fout << sol;
    return 0;
}