Pagini recente » Cod sursa (job #1366760) | Cod sursa (job #1888736) | Cod sursa (job #3037881) | Cod sursa (job #680387) | Cod sursa (job #2863580)
#include <bits/stdc++.h>
#define ull unsigned long long
using namespace std;
ifstream fin ("ssnd.in");
ofstream fout ("ssnd.out");
const ull MOD = 9973;
const ull LIM = 1e6;
bitset <LIM+5> ciur;
ull p, prim[80000], expo;
void mark(ull num){
for(ull i=num*num; i<=LIM; i+=num)
ciur[i] = true;
}
ull n, dcnt, dsum;
ull lgput(ull a, ull b){
ull answer = 1;
while(b != 0){
if(b&1)
answer *= a;
a *= a;
b >>= 1;
}
return answer;
}
signed main (){
ios_base::sync_with_stdio(false);
fin.tie(nullptr);
fout.tie(nullptr);
ciur[0] = ciur[1] = true;
mark(2);
mark(3);
for(ull i=5; i<=LIM/i; i+=6){
if(!ciur[i]) mark(i);
if(!ciur[i+2]) mark(i+2);
}
prim[++p] = 2;
prim[++p] = 3;
for(ull i=5; i<=LIM; i+=6){
if(!ciur[i]) prim[++p] = i;
if(!ciur[i+2] && i+2<=LIM) prim[++p] = i+2;
}
ull teste;
fin>>teste;
while(teste--){
fin>>n;
dcnt = dsum = 1;
for(ull i=1; i<=p; i++){
if(prim[i]*prim[i] > n) break;
if(n % prim[i] != 0) continue;
expo = 0;
while(n % prim[i] == 0){
expo++;
n /= prim[i];
}
dcnt *= (expo+1);
dsum = (ull)dsum * ((lgput(prim[i], expo+1) - 1) / (prim[i] - 1)) % MOD;
}
if(n > 1){
dcnt *= 2;
dsum = (ull)dsum * ((n*n-1) / (n-1)) % MOD;
}
fout<<dcnt<<" "<<dsum<<"\n";
}
return 0;
}