Cod sursa(job #1203060)

Utilizator cristian.caldareaCaldarea Cristian Daniel cristian.caldarea Data 30 iunie 2014 16:04:44
Problema Suma si numarul divizorilor Scor 0
Compilator cpp Status done
Runda Arhiva educationala Marime 1.36 kb
#include <fstream>
#include <bitset>

using namespace std;

const int MAX_N = 1000005;
const int MOD = 9973;

ifstream f("ssnd.in");
ofstream g("ssnd.out");

long long n;
int t, k, P[MAX_N];
bitset <MAX_N> viz;

void ciur();
inline int pow(int x, int p);
void solve();
int main() {
    ciur();
    f >> t;
        solve();
}
void ciur()
{
    for(int i = 2; i < MAX_N; ++i)
    {
        if(viz[i] == 0)
        {
            P[++k] = i;

            for(int j = i+i; j < MAX_N; j += i)
            {
                viz[j] = 1;
            }
        }
    }
}
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()
{
    f >> n;

    int nd = 1, sd = 1;

    for(int i = 1; i <= k && 1LL * P[i] * P[i] <= n; ++i) {
        if(n % P[i]) continue;
        int p = 0;

        while(n% P[i] == 0) {
            n /= 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(n > 1) {
        nd *= 2;
        sd = (1LL*sd*(n+1) % MOD);
    }

    g << nd << " " << sd << "\n";
}