Cod sursa(job #2341371)

Utilizator pistvanPeter Istvan pistvan Data 11 februarie 2019 19:33:15
Problema Suma si numarul divizorilor Scor 0
Compilator cpp-64 Status done
Runda Arhiva educationala Marime 1.64 kb
#include <iostream>
#include <fstream>
#include <cmath>
#define M 9973
using namespace std;

bool prim[1000001];
long primlist[79000];
ofstream g("ssnd.out");


void precalc()
{
    primlist[1]=2;
    primlist[0]=1;
    for (long i=2;i<1000001;i+=2)
        prim[i]=1;
    for (long i=3;i<1000001;i+=2)
    {
        if (prim[i]==0)
            primlist[++primlist[0]]=i;
        else
            continue;
        if (i>1000)
            continue;
        for (long j=i*i;j<1000001;j+=i)
            prim[j]=1;

    }
}

void dioph(long a, long b, long& alpha, long& beta)
{
    long r=a%b, alpha_p, beta_p;
    if (r==0)
    {
        alpha=0;
        beta=1;
        return;
    }
    else
    {
        dioph(b, r, alpha_p, beta_p);
        alpha=beta_p;
        beta=alpha_p-beta_p*(a/b);
        return;
    }
}

void solve(long long N)
{
    long p=1, s=1, lim=sqrt(N);
    for (long i=1;i<=primlist[0] && i<=lim;i++)
    {
        if (N%primlist[i])
            continue;
        long d=0, pr=primlist[i], prexp=1, alpha, beta;
        dioph(pr-1, M, alpha, beta);
        while (alpha<0)
            alpha+=M;
        while (N%pr==0)
        {
            d++;
            N/=pr;
            prexp*=pr;
        }
        p=(p*(d+1))%M;
        s=(s*(prexp*pr-1))%M;
        s=(s*alpha)%M;
    }
    if (N>1)
    {
        p=(2*p)%M;
        s=(s*(N+1))%M;
    }
    g<<p<<' '<<s<<'\n';
}

void read()
{
    long T;
    long long N;
    ifstream f("ssnd.in");
    f>>T;
    while (T--)
    {
        f>>N;
        solve(N);
    }
}

long main()
{
    precalc();
    read();
}