Cod sursa(job #481579)

Utilizator Mishu91Andrei Misarca Mishu91 Data 31 august 2010 21:28:51
Problema Suma si numarul divizorilor Scor 70
Compilator cpp Status done
Runda Arhiva educationala Marime 1.1 kb
#include <fstream>
#include <bitset>

using namespace std;

const int MAX_N = 1000005;
const int MOD = 9973;

ifstream fin ("ssnd.in");
ofstream fout ("ssnd.out");

long long N;
int T, K, P[MAX_N];
bitset <MAX_N> viz;

void ciur() {
	for(int i = 2; i < MAX_N; ++i) {
		if(viz[i] == 0) {
			P[++K] = i;

			for(int j = i+i; j < MAX_N; j += i) {
				viz[j] = 1;
			}
		}
	}
}

inline int pow(int x, int p) {
	long long ans = 1;
    for (int i = 30; i >= 0; -- i)
        ans = ((p >> i) & 1) ? (ans * ans * x) % MOD : (ans * ans) % MOD;
    return (int) ans;  
}
void solve() {
	fin >> N;

	int	nd = 1, sd = 1;

	for(int i = 1; i <= K && 1LL * P[i] * P[i] <= N; ++i) {
		if(N % P[i]) continue;
		int p = 0, s = 1;
		
		while(N % P[i] == 0) {
			N /= P[i];
			s *= P[i];
			s %= MOD;
			++p;
		}
		
		nd *= (p+1);
		int p1 = (s*P[i] - 1) % MOD;
		if(p1 < 0) p1 += MOD;
		int p2 = pow(P[i]-1, MOD-2) % MOD;
		if(p2 < 0) p2 += MOD;
		sd = (((sd * p1) % MOD) * p2) % MOD;
	}

	if(N > 1) {
		nd *= 2;
		sd = (1LL*sd*(N+1) % MOD);
	}

	fout << nd << " " << sd << "\n";
}

int main() {
	ciur();

	for(fin >> T; T; --T)
		solve();
}