Cod sursa(job #2661875)

Utilizator cristi_macoveiMacovei Cristian cristi_macovei Data 22 octombrie 2020 20:23:13
Problema Suma si numarul divizorilor Scor 10
Compilator cpp-64 Status done
Runda Arhiva educationala Marime 1.33 kb
#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;
}