Pagini recente » Cod sursa (job #1641365) | Cod sursa (job #396304) | Cod sursa (job #1183017) | Cod sursa (job #325135) | Cod sursa (job #799555)
Cod sursa(job #799555)
#include<fstream>
#define LL long long
#define M 9973
#define NMAX 1000002 // un numar are divizori pana la jumatatea lui => n = 10^12 =>d <= 10^6
using namespace std;
LL n, d[50], p[50], nr, sum;
int t, k,semn[NMAX + 1], ciur[NMAX + 1],lng = 0;
inline void eratosthenes()
{
ciur[lng] = 2;
lng = 1;
for(int i = 3; i <= NMAX; i += 2)
{
if(!semn[i])
{
ciur[lng] = i;
lng++;
for(int j = i + i; j + i <= NMAX; j += i)semn[j] = 0;
}
}
}
inline void descompune()
{
int putere = 1;
nr = sum = 1;
LL numar;
int i = 0;
int d = 0;
while(ciur[i] * ciur[i] <= n)
{
if(!(n % ciur[i]))
{
d = 0;
LL s = 1;
putere =1;
while(!(n % ciur[i]))
{
d++;
n /= ciur[i];
putere = putere * ciur[i];
s += putere;
}
sum = (sum * s) % M;
nr = (nr * (1+ d)) % M;
}
i++;
}
if(n > 1)nr = (2 * nr) % M, sum = (sum * (n + 1)) % M;
}
int main()
{
ifstream f("ssnd.in");
ofstream g("ssnd.out");
f>>t;
eratosthenes();
for(int r = 1; r <= t; ++r)
{
f>>n;
nr = 1;
sum = 1;
descompune();
g<<nr<<' '<<sum <<'\n';
}
g.close();
return 0;
}