Cod sursa(job #2201260)

Utilizator emiemiEmi Necula emiemi Data 4 mai 2018 00:22:52
Problema Suma si numarul divizorilor Scor 70
Compilator cpp Status done
Runda Arhiva educationala Marime 1.39 kb
#include <fstream>
#include <vector>
#include <cmath>

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

#define MOD 9973

vector<int> v;

void ciur() {
	int *a = (int *)calloc(1000000, sizeof(int));
	int i, j;

	v.push_back(2);
	int l = 1000;
	for(i = 3; i <= l; i += 2)
		if (!a[i]) {
			v.push_back(i);
			for (j = i * i; j < 1000000; j += i)
				a[j] = 1;
		}
	for (i = l + 1; i < 1000000; i += 2)
		if (!a[i])
			v.push_back(i);
}

long long putere(long long x, int n) {
	long long t = x;
	long long res = 1;
	while (n) {
		if (n & 1) {
			res *= t;
			res %= MOD;
		}
		t *= t;
		t %= MOD;
		n /= 2;
	}
	return res;
}

int main() {
	long long nrd, sd, p, put, t, n, i, ii;
	f >> t;

	ciur();
	int len = v.size();

	for (ii = 1; ii <= t; ++ii) {
		f >> n;
		if (n == 1) {
			nrd = 1;
			sd = 1;
		} else {
			nrd = 1;
			sd = 1;
			for (i = 0; i < len; ++i)
				if (n != 1 && v[i] * v[i] <= n)
					if (n % v[i] == 0) {
						n /= v[i];
						p = 1;
						while (n % v[i] == 0) {
							n /= v[i];
							++p;
						}
						nrd *= (p + 1);
						nrd %= MOD;
						sd *= (putere(v[i], p + 1) - 1);
						sd %= MOD;
						sd *= (putere(v[i] - 1, MOD - 2));
						sd %= MOD;
					}
			if (n > 1) {
				nrd *= 2;
				nrd %= MOD;
				sd *= (n + 1);
				sd %= MOD;
			}
		}
		g << nrd << ' ' << sd << '\n';
	}

	return 0;
}