Cod sursa(job #634424)

Utilizator ContraPunctContrapunct ContraPunct Data 16 noiembrie 2011 11:52:26
Problema Suma si numarul divizorilor Scor 70
Compilator cpp Status done
Runda Arhiva educationala Marime 1.18 kb
#include<fstream>
#include<bitset>
#include<vector>

const long long Nmax = 1000008;
const long long MOD = 9973;

using namespace std;

bitset<Nmax+1> ciur;
vector<long long> prime;
long long t;
long long n;
long long nr=1;
long long suma=1;

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

void Ciur ()
{
	long long i = 0 ,j = 0;
	prime.push_back(2);

	for(i=3; i<=Nmax; i+=2)
		if(!ciur[i]) 
		{
			prime.push_back(i);
			for(j = i*i; j<= Nmax; j += (i+i))
				if(!ciur[j])
					ciur[j]=1;
		}
}

void Sol()
{
	long long i,j;
	nr = 1;
	long long suma = 1;
	long long aux = n;
	for( i = 0; prime[i]*prime[i] <= aux; ++i)
	{
		if( n == prime[i] && nr == 1 && suma == 1)
		{
			fout<<2<<" "<<(1+n)%MOD<<"\n";
			return;
		}
		if( n % prime[i] == 0)
		{
			long long exp; 
			long long n2 = n;
			for (exp=0; n%prime[i]==0; n/=prime[i])
			{
				++exp;
			}
			nr = nr * (exp+1);
			suma = suma * ((n2/n * prime[i]) - 1 ) / ( prime[i] - 1) % MOD;
			suma = suma  % MOD;
		}
	}
	if( n > 1)
	{
		nr*=2;
		suma=suma*(n*n-1)/(n-1) % MOD;
	}
	fout<<nr<<" "<<suma<<"\n";
}

int main()
{
	Ciur();
	long long i;
	fin>>t;
	for( i=0;i<t;++i)
	{
		fin>>n;
		Sol();
	}
	fout.close();
	fin.close();
	return 0;
}