Pagini recente » Cod sursa (job #1144853) | Cod sursa (job #1544645) | Cod sursa (job #2099150) | Cod sursa (job #299045) | Cod sursa (job #2935933)
#include <bits/stdc++.h>
#define NMAX 1000005
#define MOD 9973
using namespace std;
vector <int> primes;
int lp[NMAX];
int logpow(long long n, int p) {
int ans = 1;
n %= MOD;
while (p) {
if (p & 1)
ans = (1LL * ans * n) % MOD;
n = (1LL * n * n) % MOD;
p >>= 1;
}
return ans;
}
int invmod(long long n) {
return logpow(n, MOD - 2);
}
void make_sieve() {
for (int i = 2; i < NMAX; ++i) {
if (lp[i] == 0)
primes.push_back(i), lp[i] = i;
for (int j = 0; j < (int)primes.size() && i * primes[j] < NMAX; ++j) {
lp[i * primes[j]] = primes[j];
if (i % primes[j] == 0)
break;
}
}
}
void solve() {
long long n;
cin >> n;
vector <pair<long long, int>> decomp;
for (int i = 0; i < (int)primes.size() && 1LL * primes[i] * 1LL * primes[i] <= n; ++i) {
if (n % primes[i] == 0) {
int cnt = 0;
while (n % primes[i] == 0)
n /= primes[i], ++cnt;
decomp.push_back({primes[i], cnt});
}
}
if (n != 1)
decomp.push_back({n, 1});
long long cntdiv = 1, sumdiv = 1;
for (auto p : decomp) {
cntdiv *= (p.second + 1);
sumdiv = (sumdiv * (logpow(p.first, p.second + 1) - 1)) % MOD;
sumdiv = (sumdiv * invmod(p.first - 1)) % MOD;
}
cout << cntdiv << ' ' << sumdiv << '\n';
}
int main() {
make_sieve();
freopen("ssnd.in", "r", stdin);
freopen("ssnd.out", "w", stdout);
int t;
cin >> t;
while (t--)
solve();
return 0;
}