Cod sursa(job #1827353)

Utilizator AndreiMaximIonutMaxim Andrei AndreiMaximIonut Data 11 decembrie 2016 19:20:02
Problema Suma si numarul divizorilor Scor 10
Compilator cpp Status done
Runda Arhiva educationala Marime 1.11 kb
#include <fstream>
#include <bitset>
using namespace std;
ifstream fin("ssnd.in");
ofstream fout("ssnd.out");
#define Nmax 1000005
int t,k,P[Nmax];
bitset <Nmax> v;
long long n;
const int MOD = 9973;
void ciur()
{
    long i,j;
    v[0]=v[1]=1;
    for(i=2; i*i<Nmax; i++)
        if(v[i]==0)
        {
            P[++k]=i;
            for(j=i; j*i<Nmax; j++)
                v[i*j]=1;
        }
}

inline int pow(int x, int p)
{
    int d=1,i;
    for(i=1; i<=p; i++)
        d=d*x,d=d%MOD;
    return d;
}
void solve()
{
    fin>>n;
    int i,nd=1,sd=1;
    for(i=1; i<=k && P[i]*P[i]<=n; i++)
        if(n%P[i]==0)
        {
            int e=0;
            while(n%P[i]==0)
                e++,n/=P[i];
            int a=pow(P[i],e+1)%MOD;
            nd=nd*(e+1);
            sd=sd*(((a-1)/(P[i]-1))%MOD)%MOD;
            sd=sd%MOD;

        }
    if(n>1)
    {
        nd*=2;
        sd=sd*((((n*n-1)%MOD)/(n-1))%MOD)%MOD;
    }
    fout<<nd<<' '<<sd<<endl;

}
int main()
{
    fin>>t;
    ciur();
    while(t)
    {
        solve();
        t--;
    }
    return 0;
}