Cod sursa(job #1937775)

Utilizator GoogalAbabei Daniel Googal Data 24 martie 2017 11:31:09
Problema Suma si numarul divizorilor Scor 100
Compilator cpp Status done
Runda Arhiva educationala Marime 1.71 kb
#include <fstream>
#define MOD 9973
#define NMAX 1000005

using namespace std;

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

long long t;
long long n,k,prim[NMAX],r1,r2;
bool ciur[NMAX];

void ciur_er()
{
    long long i,j;
    for(i=2; i<NMAX; i++)
    {
        if(!ciur[i])
        {
            prim[++k]=i;
            for(j=i+i; j<=NMAX; j+=i)
                ciur[j]=true;
        }
    }
}

long long ridic_log(long long n,long long p)
{
    long long a=n,i,sol=1;
    for(i=0;(1<<i)<=p;i++)
    {
        if(((1<<i)&p)>0)
            sol=(sol*a)%MOD;
        a=(a*a)%MOD;
    }

    return sol;
}

void inv_mod(long long &x, long long &y, long long a, long long b)
{
    if(!b)
    {
        x=1;
        y=0;
    }

    else
    {
        inv_mod(x,y,b,a%b);
        swap(x,y);
        y-=x*(a/b);
    }
}

long long turn_it(long long a)
{
    long long inv=0,ins;
    inv_mod(inv,ins,a,MOD);
    if(inv<=0)
        inv=MOD+inv%MOD;
    return inv;
}

int main()
{
    long long i,j,p,p1,p2;
    ciur_er();
    fin>>t;
    for(i=1; i<=t; i++)
    {
        fin>>n;
        r1=r2=1;
        for(j=1; j<=k && prim[j]*prim[j]<=n; j++)
        {
            if(n%prim[j])
                continue;
            p=0;
            while(!(n%prim[j]))
            {
                p++;
                n/=prim[j];
            }

            r1*=(p+1);

            p1=(ridic_log(prim[j],p+1)-1)%MOD;
            p2=turn_it(prim[j]-1)%MOD;

            r2=(((r2*p1)%MOD )* p2)%MOD;
        }

        if(n>1)
        {
            r1*=2;
            r2=(r2*(n+1))%MOD;
        }

        fout<<r1<<' '<<r2<<'\n';
    }
    return 0;
}