Cod sursa(job #2082785)

Utilizator xRoALexBirtoiu Alexandru xRoALex Data 6 decembrie 2017 19:31:05
Problema Suma si numarul divizorilor Scor 10
Compilator cpp Status done
Runda Arhiva educationala Marime 1.41 kb
#include <cstdio>
#include <vector>
#define MOD 9973
using namespace std;

FILE *in = fopen("ssnd.in","r");
FILE *out = fopen("ssnd.out","w");

int prim[1000005];
int t;
vector <int> num;

void ciur()
{
    prim[1]=1;
    num.push_back(2);
    for(int i=3; i*i <= 1000000; i+=2)
        if(prim[i]==0)
            {
                num.push_back(i);
                for(int j=i+i; j <= 1000000; j+=i)
                prim[j]=1;
                }

}

int putere(long long x,long long e)
{
    long long r=1;
    for(int d=0; (1<<d) <= e; d++)
    {
        if((1<<d)&e)
        {
            r=(r*x)%MOD;
        }
        x=(x*x)%MOD;
    }
    return r%MOD;
}

int invmod(long long x)
{
    return putere(x,MOD-2)%MOD;
}

void desc(long long x)
{
    int e,k=0,p=1,suma=1;
    while(x!=1 && num[k]*num[k] <= x)
    {
        e=0;
        while(x%num[k]==0)
        {
            e++;
            x/=num[k];
        }
        p*=e+1;
        if(e>0)
        suma = (suma * ((putere(num[k],e + 1)-1) * invmod(num[k]-1))) % MOD;
        k++;
    }
    if(x!=1)
    {
        suma = (suma * ((putere(x,2)-1) * invmod(x-1))) % MOD;
        p*=2;
    }
    fprintf(out,"%d %d\n",p,suma);
}

int main()
{
    fscanf(in,"%d",&t);
    ciur();
    while(t)
    {
        long long x;
        fscanf(in,"%lld",&x);
        desc(x);
        t--;
    }
    return 0;
}