Pagini recente » Cod sursa (job #1489422) | Cod sursa (job #2245469) | Cod sursa (job #2875972) | Cod sursa (job #1954858) | Cod sursa (job #799560)
Cod sursa(job #799560)
#include<fstream>
#define LL long long
#define M 9973
#define NMAX 1000000 // 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 + 8], ciur[NMAX + 8],lng = 0;
inline void eratosthenes()
{
for(int i = 2; i <= NMAX; ++i)
{
if(!semn[i])
{
ciur[lng] = i;
lng++;
for(int j = i; j <= NMAX; j += i)semn[j] = 1;
}
}
}
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;
}