Pagini recente » Cod sursa (job #1641678) | Cod sursa (job #1641660) | Cod sursa (job #1548626) | Cod sursa (job #2819479) | Cod sursa (job #1937775)
#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;
}