Pagini recente » Cod sursa (job #545863) | Cod sursa (job #2179451) | Cod sursa (job #1648687) | Cod sursa (job #1601284) | Cod sursa (job #2341388)
#include <iostream>
#include <fstream>
#include <cmath>
#define M 9973
using namespace std;
bool prim[1000001];
long primlist[79000];
ofstream g("ssnd.out");
void precalc()
{
int i, j;
primlist[1]=2;
primlist[0]=1;
for (i=2;i<1000001;i+=2)
prim[i]=1;
for (i=3;i<1000001;i+=2)
{
if (prim[i]==0)
primlist[++primlist[0]]=i;
else
continue;
if (i>1000)
continue;
for (j=i*i;j<1000001;j+=i)
prim[j]=1;
}
}
void dioph(long a, long b, long& alpha, long& beta)
{
long r=a%b, alpha_p, beta_p;
if (r==0)
{
alpha=0;
beta=1;
return;
}
else
{
dioph(b, r, alpha_p, beta_p);
alpha=beta_p;
beta=alpha_p-beta_p*(a/b);
return;
}
}
void solve(long long N)
{
long p=1, s=1, lim=sqrt(N);
for (int i=1;i<=primlist[0] && i<=lim;i++)
{
if (N%primlist[i])
continue;
long d=0, pr=primlist[i], prexp=1, alpha, beta;
dioph(pr-1, M, alpha, beta);
while (alpha<0)
alpha+=M;
while (N%pr==0)
{
d++;
N/=pr;
prexp*=pr;
}
p=(p*(d+1))%M;
s=(s*(prexp*pr-1))%M;
s=(s*alpha)%M;
}
if (N>1)
{
p=(2*p)%M;
s=(s*(N+1))%M;
}
g<<p<<' '<<s<<'\n';
}
void read()
{
long T;
long long N;
ifstream f("ssnd.in");
f>>T;
while (T--)
{
f>>N;
solve(N);
}
}
int main()
{
precalc();
read();
}