Cod sursa(job #2103286)

Utilizator dragos.galeteanu2001Dragos Iulian dragos.galeteanu2001 Data 9 ianuarie 2018 23:46:33
Problema Suma si numarul divizorilor Scor 100
Compilator cpp Status done
Runda Arhiva educationala Marime 1.17 kb
#include <fstream>
#define MOD 9973
#define limit 1000005

using namespace std;

ifstream f("ssnd.in");
ofstream g("ssnd.out");

long long int n, k = -1, DA[limit];
int t;
bool prim[limit];

void ciur()
{
	for (long long int i = 2; i < limit; ++i) {
		if (prim[i] == 0) {
			DA[++k] = i;

			for (long long int j = i + i; j < limit; j += i) {
				prim[j] = 1;
			}
		}
	}
}

long long int power(long long int x, long long int p)
{
	long long int rez = 1;
	while (p) {
		if (p & 1) rez = (rez*x) % MOD;
		x = (x*x) % MOD;
		p >>= 1;
	}
	return rez % MOD;
}

void solve()
{
	f >> n;
	long long int nr = 1, sum = 1, pw, i, pw1, pw2;
	for (i = 0; i <= k && DA[i] * DA[i] <= n; i++) {
		if (n % DA[i]) continue;
		pw = 0;
		while (n % DA[i] == 0) {
			n /= DA[i];
			pw++;
		}
		nr *= (pw + 1);
		pw1 = (power(DA[i], pw + 1) - 1) % MOD;
		pw2 = (power(DA[i] - 1, MOD - 2)) % MOD;
		sum = (((sum*pw1) % MOD) * pw2) % MOD;
	}

	if (n > 1) {
		nr *= 2;
		sum = (sum * (n + 1)) % MOD;
	}

	g << nr << ' ' << sum << '\n';
}

int main()
{
	ciur();
	f >> t;
	while (t) {
		solve();
		t--;
	}
	f.close();
	g.close();
	return 0;
}