Pagini recente » Cod sursa (job #1080442) | Cod sursa (job #1595153) | Cod sursa (job #1515550) | Cod sursa (job #2676588) | Cod sursa (job #2610718)
#include <bits/stdc++.h>
using namespace std;
ifstream fin("ssnd.in");
ofstream fout("ssnd.out");
const int NMAX = 1000005;
const int MOD = 9973;
int primes[NMAX], c, d, sum_div, nr_div, power;
bool ciur[NMAX / 2];
long long val;
void init()
{
for (int i = 2; i * i < NMAX; i++)
if (!ciur[i])
{
primes[++c] = i;
for (int j = i * i; j < NMAX; j += i)
ciur[j] = 1;
}
}
int exp_log(int n, int p)
{
int r = 1;
n %= MOD;
while(p)
{
if(p % 2)
r = (r * n) % MOD;
n = n * n;
n %= MOD;
p /= 2;
}
return r;
}
void divizori(long long numar)
{
nr_div = sum_div = 1;
int j = 1;
d = primes[j];
while (numar > 1)
{
power = 0;
while (!(numar % d))
{
numar /= d;
power++;
}
if (power)
{
nr_div *= (power + 1);
sum_div = sum_div * (exp_log(d, power + 1) - 1) / (d - 1);
}
j++;
d = primes[j];
if(d * d > numar && numar > 1)
{
nr_div *= 2;
sum_div = (sum_div * (numar + 1) * 1LL) % MOD;
numar = 1;
}
}
fout << nr_div << ' ' << sum_div << '\n';
}
int main()
{
int nr; fin >> nr;
init();
for (int i = 0; i < nr; i++)
{
fin >> val;
divizori(val);
}
return 0;
}