Pagini recente » Cod sursa (job #519974) | Cod sursa (job #1138104) | Cod sursa (job #1374757) | Cod sursa (job #2470135) | Cod sursa (job #482346)
Cod sursa(job #482346)
#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;
while(n % div == 0)
{
p=div;
d++;
n/=p;
}
}
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(600001);
v.push_back(2);
for(int i=1; i<600001; ++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;
uint64 i=0,aux=n;
uint64 s=1,numdiv=1;
while(n>1 && v[i]*v[i]<=aux)
{
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;
}
if(n>1)
{
s=n+1;
numdiv=2;
}
fout<<numdiv<<" "<<s<<"\n";
}
fin.close();
fout.close();
return 0;
}