Pagini recente » Cod sursa (job #1598577) | ioit-werebears-01 | Cod sursa (job #1156775) | Cod sursa (job #479838) | Cod sursa (job #2495207)
#include <fstream>
#include <bitset>
#define mod 9973
using namespace std;
ifstream fin("ssnd.in");
ofstream fout("ssnd.out");
long long n,r;
int t,i,j,p[100000],m,sum,nr,cnt,x,d;
bitset <1000001> f;
int pow(int a, int b){
r=1;
a%=mod;
while(b){
if(b&1)
r=(r*a)%mod;
a=(a*a)%mod;
b>>=1;
}
return r;
}
int main(){
p[++m]=2;
for(i=3;i<1000000;i+=2){
if(!f[i]){
p[++m]=i;
for(j=3*i;j<=1000000;j+=2*i)
f[j]=1;
}
}
for(fin>>t;t;t--){
fin>>n;
sum=1; nr=1;
for(i=1;p[i]<=n/p[i] && i<=m;i++){
if(n%p[i]==0){
cnt=1; d=p[i];
while(n%d==0){
cnt++;
n/=d;
}
nr=(nr*cnt)%mod;
x=pow(d,cnt)-1;
if(x<0)
x+=mod;
sum=(sum*x)%mod;
//sum=(sum* ( (pow(p[i],cnt+1)-1<0) ? pow(p[i],cnt+1)-1+mod : pow(p[i],cnt+1)-1 ) )%mod;
sum=(sum*pow((d-1+mod)%mod,mod-2))%mod;
}
}
if(n!=1){
nr*=2;
n%=mod;
x=pow(n,2)-1;
if(x<0)
x+=mod;
sum=(sum*x)%mod;
//sum=(sum* ( (pow(n,2)-1<0) ? pow(n,2)-1+mod : pow(n,2)-1 ) )%mod;
sum=(sum*pow((n-1+mod)%mod,mod-2))%mod;
}
fout<<nr<<" "<<sum<<"\n";
}
return 0;
}