Pagini recente » Cod sursa (job #599959) | Cod sursa (job #3143365) | Rating Valeriu Cojocari (Valeriu_Cojocari) | Cod sursa (job #327292) | Cod sursa (job #1988938)
#include <iostream>
#include <cmath>
#include <fstream>
using namespace std;
ifstream f("ssnd.in");
ofstream g("ssnd.out");
int vp[78500];
long long x,nr,s;
long long Exp_Log(long long a,long long b)
{
long long rez=1;
while(b)
{
if(b%2)
{
rez*=a;
b--;
}
else
{
a*=a;
b/=2;
}
}
return rez;
}
long long Exp_Log_Mod(long long a,long long b,long long pr)
{
long long rez=1;
while(b)
{
if(b%2)
{
rez=(rez%pr)*(a%pr)%pr;
b--;
}
else
{
a=(a%pr)*(a%pr)%pr;
b/=2;
}
}
return rez;
}
long long Fractie(long long d,long long p,long long pr)
{
long long nmrt,numi;
nmrt=Exp_Log_Mod(d,p,pr)-1;
numi=Exp_Log_Mod(d-1,pr-2,pr);
return (nmrt*numi)%pr;
}
void Nr_Sum_Div()
{
int r,p,i=1;
long long put;
nr=1;
s=1;
r=sqrt(x);
while(vp[i]<=r)
{
p=0;
while(x%vp[i]==0)
{
p++;
x/=vp[i];
}
if(p)
{
nr*=p+1;
s*=Fractie(vp[i],p+1,9973);
// put=Exp_Log(vp[i],p+1);
// s*=(put-1)/(vp[i]-1);
r=sqrt(x);
}
i++;
}
if(x!=1)
{
nr*=2;
s*=x+1;
}
}
void Ciur()
{
int n=0;
int i,j;
bool c[1000000];
for(i=2;i<=1000000;i++)
if(!c[i])
{
vp[++n]=i;
for(j=2*i;j<=1000000;j+=i)
c[j]=1;
}
}
int main()
{
int i,t;
Ciur();
f>>t;
for(i=1;i<=t;i++)
{
f>>x;
Nr_Sum_Div();
g<<nr<<" "<<s<<endl;
}
f.close();
g.close();
return 0;
}