Pagini recente » Cod sursa (job #3221801) | Cod sursa (job #1939133) | Cod sursa (job #2984240) | Cod sursa (job #2479687) | Cod sursa (job #2661875)
#include <iostream>
#include <fstream>
#include <vector>
const int MOD = 9973;
const int NMAX = 1e6;
bool sieve[1 + NMAX];
std::vector<int64_t> primes;
inline void sieve_of_eratosthenes() {
primes.push_back(2);
for (int64_t i = 3; i <= NMAX; i += 2) {
if (!sieve[i]) {
primes.push_back(i);
for (int64_t j = i * i; j <= NMAX; j += 2 * i)
sieve[j] = true;
}
}
}
inline int64_t pow(int64_t base, int64_t exp) {
if (exp == 0)
return 1;
if (exp % 2 == 1)
return (base * pow(base, exp - 1)) % MOD;
int64_t temp = pow(base, exp / 2);
return (temp * temp) % MOD;
}
int main() {
std::ifstream in("ssnd.in");
std::ofstream out("ssnd.out");
sieve_of_eratosthenes();
int64_t n, a = 1;
int64_t exp, num, sum;
in >> n;
while (in >> a) {
num = sum = 1;
for (int64_t div_ind = 0; div_ind < primes.size() && a != 1; ++div_ind) {
exp = 0;
while (a % primes[div_ind] == 0) {
a /= primes[div_ind];
++ exp;
}
num = (num * (exp + 1)) % MOD;
sum = (int64_t)((sum * ((int64_t)pow(primes[div_ind], exp + 1) - 1) / (primes[div_ind] - 1)) % MOD);
}
if (a != 1) {
exp = 1;
num = (num * (exp + 1)) % MOD;
sum = (sum * ((int64_t)pow(a, exp + 1) - 1) / (a - 1)) % MOD;
}
out << num << ' ' << sum << '\n';
}
return 0;
}