Cod sursa(job #930789)

Utilizator gabriel.badeaGabriel Badea gabriel.badea Data 27 martie 2013 20:16:49
Problema Suma si numarul divizorilor Scor 100
Compilator cpp Status done
Runda Arhiva educationala Marime 1.06 kb
#include <fstream>
#include <bitset>

using namespace std;

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

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

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

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)
{
	int rez = 1;
	x %= MOD;

	for(;p;p>>=1)
	{
		if(p&1)
		{
			rez *= x;
			rez %= MOD;
		}

		x *= x;
		x %= MOD;
	}
	return rez;
}

void sol()
{
	f >> 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;

		while(N % P[i] == 0)
		{
			N /= P[i];
			++p;
		}

		nd *= (p+1);

		int p1 = (pow(P[i], p+1) - 1) % MOD;
		int p2 = pow(P[i]-1, MOD-2) % MOD;

		sd = (((sd * p1) % MOD) * p2) % MOD;
	}

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

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

int main()
{
	ciur();
	for(f>>T;T;--T)
		sol();
}