Cod sursa(job #1095216)

Utilizator Dddarius95Darius-Florentin Neatu Dddarius95 Data 30 ianuarie 2014 16:12:06
Problema Suma si numarul divizorilor Scor 100
Compilator cpp Status done
Runda Arhiva educationala Marime 1.22 kb
#include <fstream>
#include <bitset>
#define Nmax 1000000
#define X 9973
using namespace std;
ifstream f("ssnd.in");
ofstream g("ssnd.out");

int T,P[100000];
long long N,D,S;
bitset < Nmax+99 > v;

void Ciur()
{
    P[++P[0]]=2;
    for(int i=3;i<=Nmax;i+=2)
        if(!v[i])
        {
            P[++P[0]]=i;
            if(1LL*i*i<=Nmax)
                for(int j=i*i;j<=Nmax;j+=i)v[j]=1;
        }
}

long long Putere(long long N,long long P)
{

    long long m=1;
    while (P!=1)
    if(P % 2==0)
    {
        N=(N*N) % X;
        P/=2;
    }
    else
    {
        m=(m*N) % X;
        P--;
    }
    return (m*N)% X;
}

inline void Querry()
{
    f>>N; D=1,S=1;
    for(int i=1;i<=P[0] && 1LL*P[i]*P[i]<=N;++i)
        if(N % P[i]==0)
        {
            int k=0;
            while(N % P[i]==0)N/=P[i],++k;
            D*=(k+1);
            int x1=(Putere(P[i],k+1)-1) % X;
            int x2=Putere(P[i]-1,X-2);
            S=(( (S*x1) % X ) * x2) % X;
        }
    if(N>1)
    {
        D*=2;
        S=(1LL*S*(N+1) % X);
    }
    g<<D<<' '<<S<<'\n';
}

int main()
{
    Ciur();
    f>>T;
    for(int i=1;i<=T;++i)Querry();
    f.close();g.close();
    return 0;
}