Cod sursa(job #480380)

Utilizator paul_gabryelPaul Buda paul_gabryel Data 27 august 2010 16:21:54
Problema Suma si numarul divizorilor Scor 0
Compilator cpp Status done
Runda Arhiva educationala Marime 0.82 kb

#include<cstdio>
#include<fstream>

using namespace std;

#define nn 999984
#define nm 78499
#define mod 9973

int c[nn],v[nm];
int nrp,t;
long long n;

int ridica (int baz,int exp)
{
baz%=mod;
int b[50],nrb=0;
for(;exp;exp/=2)
b[++nrb]=exp%2;
int rez=1;
for(;nrb;--nrb){
rez=(rez*rez)%mod;
if(b[nrb])
rez=(rez*baz)%mod;
}

return rez;
}

int main ()
{

for(int i=2;i<nn;++i)
if(c[i]==0){
v[++nrp]=i;
for(int j=2*i;j<nn;j+=i)
c[j]=1;
}

ifstream in ("ssnd.in");
freopen("ssnd.out","w",stdout);
for(in>>t;t;--t){
in>>n;
int i=1,s=1,nrd=1;
while(n>1){
int p=0;
while(n%v[i]==0){
++p;
n/=v[i];
}
s=(s*(ridica(v[i],p+1)-1)/(v[i]-1))%mod;
nrd*=(p+1);
++i;
if(v[i]*v[i]>n&&v[i]&&n%v[i]!=0){
nrd<<1;
s=((s*n*n)%mod);
n=1;
}
}
printf("%d %d\n",nrd,s);
}
in.close();


return 0;}