Cod sursa(job #2130162)

Utilizator MattCMatei Coroiu MattC Data 13 februarie 2018 14:56:57
Problema Suma si numarul divizorilor Scor 40
Compilator cpp Status done
Runda Arhiva educationala Marime 1.31 kb
#include <fstream>
#include <bitset>

using namespace std;
const int Nmax = 1e6 + 5;
const int MOD = 9973;

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

bitset<Nmax>ciur;
int P[Nmax],K,x;

inline int pow(int x, int p)
{
    int rez = 1;
    x %= MOD;

    for(; p; p >>= 1)
    {
        if(p & 1)
        {
            rez *= x;
            rez %= MOD;
        }

        x *= x;
        x %= MOD;
    }

    return rez;
}

void solve()
{
    int i,j,nd=1,sd=1,p;
    fin >> x;
    for (i=1; i<=K && 1LL*P[i]*P[i]<=x; ++i)
    {
        if (x%P[i]) continue;
        p=0;
        while (x%P[i]==0)
        {
            x /= P[i];
            ++p;
        }
        nd *= (p+1);
        int p1 = (pow(P[i], p+1) - 1) % MOD;
        int p2 = pow(P[i]-1, MOD-2) % MOD;
        sd = (((sd * p1) % MOD) * p2) % MOD;
    }
    if (x>1)
    {
        nd = nd << 1;
        sd = (1LL*sd*(x+1) % MOD);
    }
    fout << nd << " " << sd << "\n";
}

int main()
{
    int i,j,n;
    ciur[0]=ciur[1]=1;
    for (i=2; i<Nmax; ++i)
    {
        if (!ciur[i])
        {
            P[++K]=i;
            for (j= (i << 1); j<Nmax; j+=i)
            {
                ciur[j]=1;
            }
        }
    }
    fin >> n;
    while (n--)
        solve();
    return 0;
}