Pagini recente » Cod sursa (job #1507688) | Profil M@2Te4i | Cod sursa (job #764764) | Rating cora alexandru (coraalexandru) | Cod sursa (job #481618)
Cod sursa(job #481618)
#include <fstream>
#include <bitset>
using namespace std;
const int MAX_N = 1000005;
const int MOD = 9973;
ifstream fin ("ssnd.in");
ofstream fout ("ssnd.out");
long long N;
int T, K, P[MAX_N];
//bitset <MAX_N> viz;
unsigned char viz[1 << 17];
void ciur ( void ) {
P[++K] = 2;
for (int i = 1; ( i * i << 1 ) + ( i << 1 ) <= MAX_N; ++i) {
if ( ( viz[i >> 3] & ( 1 << ( i & 7 ) ) ) ) {
continue;
}
for (int j = ( i * i << 1 ) + ( i << 1 ); ( j << 1 ) + 1 <= MAX_N; j += ( i << 1 ) + 1) {
viz[j >> 3] |= ( 1 << ( j & 7 ) );
}
}
for (int i = 1; ( i << 1 ) + 1 <= MAX_N; ++i) {
if ( ( viz[i >> 3] & ( 1 << ( i & 7 ) ) ) == 0 ) {
P[++K] = ( i << 1 ) + 1 ;
}
}
}
inline int pow ( int N, int P ) {
int sol = 1;
for ( int i = 0; 1 << i <= P; ++i ) {
if ( ( 1 << i & P ) > 0 ) {
sol *= N, sol %= MOD ;
}
N *= N, N %= MOD ;
}
return sol ;
}
inline int pow1(int x, int p) {
int rez = 1;
x %= MOD;
for(; p; p >>= 1) {
if(p & 1) {
rez *= x;
rez %= MOD;
}
x *= x;
x %= MOD;
}
return rez;
}
void solve() {
fin >> N;
int nd = 1, sd = 1;
for(int i = 1; i <= K && 1LL * P[i] * P[i] <= N; ++i) {
if(N % P[i]) continue;
int p = 0;
while(N % P[i] == 0) {
N /= P[i];
++p;
}
nd *= (p+1);
int p1 = (pow(P[i], p+1) - 1) % MOD;
int p2 = pow(P[i]-1, MOD-2) % MOD;
sd = (((sd * p1) % MOD) * p2) % MOD;
}
if(N > 1) {
nd *= 2;
sd = (1LL*sd*(N+1) % MOD);
}
fout << nd << " " << sd << "\n";
}
int main() {
ciur();
for(fin >> T; T; --T)
solve();
}