Pagini recente » Cod sursa (job #2228104) | Cod sursa (job #948212) | Istoria paginii runda/335900/clasament | Cod sursa (job #1141176) | Cod sursa (job #1165454)
#include <cstdio>
#include <bitset>
using namespace std;
#define Nmax 1000005
#define MOD 9973
bitset<Nmax> used;
int N,nrp;
int is_prime[80000];
void precalc()
{
is_prime[nrp++] = 2;
for(int i = 1; ((i<<1)|1) <= Nmax; ++i)
if(!used[((i<<1)|1)])
{
is_prime[ nrp++ ]= ((i<<1)|1);
for(int j = 1;((i<<1)|1)*((j<<1)|1) <= Nmax ; ++j)
used[((i<<1)|1)*((j<<1)|1)] = 1;
}
}
int lg_put(int x, int b)
{
if (b == 0) return 1;
if (b == 1) return x;
int x1=x,x2=1;
while(b != 1){
if(b&1)
x2 = (1LL*x2 * x1)%MOD,b^=1;
x1 = (1LL*x1 * x1)%MOD,b>>=1;
}
return (x1*x2)%MOD;
}
void solve()
{
long long a,b;
int nrd = 1,times = 0,sd = 1,x,y;
scanf("%lld",&a);b = a;
for(int i = 0; is_prime[i] *is_prime[i] <= a; ++i)
{
times = 0;
while(a % is_prime[i] == 0)
{
++times;
a/= is_prime[i];
}
if(times)
{
nrd = nrd * (times + 1);
x = (lg_put(is_prime[i],times + 1) - 1) % MOD ;
y = (lg_put(is_prime[i]-1, MOD-2)) % MOD ;
sd = (1LL*(( 1LL * sd * x )% MOD) * y) % MOD;
}
if(a == 1) break;
}
if(a > 1){nrd *= 2;sd = (1LL*sd*(a + 1))%MOD;}
printf("%d %d\n",nrd,sd);
}
void read()
{
scanf("%d",&N);
while(N--)solve();
}
int main()
{
freopen("ssnd.in","r",stdin);
freopen("ssnd.out","w",stdout);
precalc();
read();
return 0;
}