Cod sursa(job #1951438)

Utilizator PaulStighiStiegelbauer Paul-Alexandru PaulStighi Data 3 aprilie 2017 16:48:38
Problema Suma si numarul divizorilor Scor 40
Compilator cpp Status done
Runda Arhiva educationala Marime 1.4 kb
#include<fstream>
#include<bitset>
using namespace std;
ifstream fin("ssnd.in");
ofstream fout("ssnd.out");

const int CMax = 1000005;
const int MOD = 9973;

int T,K,N,Nr,Sum;
int Prim[CMax];
bitset <CMax> Ciur;

void Sieve()
{
    for(int i = 2 ; i < CMax ; ++i)
        if(!Ciur[i])
        {
            Prim[++K] = i;

            for(int j = i + i ; j < CMax ; j += i)
                Ciur[j] = 1;
        }
}

inline int put(int x,int p)
{
    int sol = 1;

    x %= MOD;

    while(p)
    {
        if(p & 1)   sol = (sol * x) % MOD;
        x = (x * x) % MOD;
        p >>= 1;
    }

    return sol;
}

void Solve()
{
    Nr = Sum = 1;

    for(int i = 1 ; i <= K && (1LL * Prim[i] * Prim[i]) <= N ; ++i)
    {
        if(N % Prim[i]) continue;

        int pow = 1;

        while(!(N % Prim[i]))
        {
            N /= Prim[i];
            pow++;
        }

        Nr *= pow;

        int up,down;

        up = (put(Prim[i],pow) - 1) % MOD;
        down = put(Prim[i] - 1,MOD - 2) % MOD;

        Sum = (((Sum * up) % MOD) * down) % MOD;
    }

    if(N > 1)
    {
        Nr *= 2;
        Sum = (1LL * Sum * (N + 1)) % MOD;
    }
}

int main()
{
    Sieve();

    fin>>T;

    while(T--)
    {
        fin>>N;

        Solve();

        fout<<Nr<<" "<<Sum<<"\n";
    }

    fin.close();
    fout.close();
    return 0;
}