Cod sursa(job #482347)

Utilizator claudiumihailClaudiu Mihail claudiumihail Data 3 septembrie 2010 12:08:51
Problema Suma si numarul divizorilor Scor 0
Compilator cpp Status done
Runda Arhiva educationala Marime 1.8 kb
#include<fstream>
#include<iostream>
#include<bitset>
#include<vector>
using namespace std;

typedef unsigned long long uint64;

const unsigned int MOD = 9973;

template<unsigned long numbits>
struct PrimeSieve
{
	long long num;
	bitset<numbits> bits;

	void GeneratePrimes(const long long n)
	{
		num=1;
		for(int i=1; ((i<<1)+1)*((i<<1)+1)<n; ++i)
		{
			if(!bits[i])
			{
				for(long long j=i; ((i<<1)+1)*((j<<1)+1)<n; ++j)
				{
					const long long index=((i<<1)+1)*((j<<1)+1);
					bits[index>>1]=1;
				}
			}
		}
	}
	
	const bitset<numbits>& GetPrimes() const
	{
		return bits;
	}
};

uint64 powl(uint64 n, int p)
{
	uint64 rez=1;
	while(p)
	{
		if(p&1)
		{
			rez*=n;
		}
		n*=n;
		p>>=1;
	}
	return rez;
}

void Decompose(uint64& n, uint64 div, uint64& p, uint64& d)
{
	d=p=0;
	if(n % div == 0)
	{
		p=div;
		while(n % div == 0)
		{
			d++;
			n/=div;
		}
	}
}

int main()
{
	int t;
	uint64 n,p,d;
	fstream fin("ssnd.in", fstream::in);
	fstream fout("ssnd.out", fstream::out);
	vector<unsigned long> v;
	
	PrimeSieve<1000007> primes;
	primes.GeneratePrimes(500005);
	v.push_back(2);
	for(int i=1; i<500005; ++i)
	{
		if(!(primes.GetPrimes())[i])
		{
			v.push_back((i<<1)+1);
		}
	}
	
	fin>>t;
	//cout<<t<<endl;
	
	
	for(int l=0; l<t; ++l)
	{
		fin>>n;
		//n=12;
		uint64 i=0,aux=n;
		uint64 s=1,numdiv=1;
		while(n>1 && v[i]*v[i]<=n)
		{
			Decompose(n, v[i], p, d);
			if(p)
			{
				//cout<<p<<" - "<<d<<endl;
				numdiv*=(d+1);
				//cout<<powl(p,d)<<endl<<endl;
				s=(s*((powl(p,d+1)-1)/(p-1)))%MOD;
			}
			i++;
			//cout<<(i<<1)+1<<endl;
		}
		//cout<<n<<" -- "<<v[i]<<endl;
		if(n>1)
		{
			s=(s*(n+1))%MOD;
			numdiv=2;
		}
		
		fout<<numdiv<<" "<<s<<"\n";
	}
	
	fin.close();
	fout.close();
	return 0;
}