Pagini recente » Cod sursa (job #654566) | Cod sursa (job #303776) | Cod sursa (job #657990) | Cod sursa (job #642990) | Cod sursa (job #862502)
Cod sursa(job #862502)
#include <stdio.h>
#include <math.h>
#include <bitset>
using namespace std;
FILE *f=fopen("ssnd.in","r");
FILE *g=fopen("ssnd.out","w");
const int n=1000000;
bitset <n> prim;
long long x,nr2;
unsigned long long b,nr,i,a,p,cnt=0,v[100000],cp;
void ciur()
{
int i,j;
for (i=2;i<=n;i++)
if (prim[i]==0)
{
v[++b]=i;
cnt++;
for (j = 2*i;j<=n;j+=i)
prim[j] = 1;
}
}
void ridic(long long nr, long long put)
{
int i;
long long cn=nr,m;
cp=1;
for(i=0;i<32;i++)
{
m=1<<i;
if((m&put)==m)
{
cp=(cn*cp)%9973;
}
cn=(cn*cn)%9973;
}
}
void desc()
{
long long x,dv=1,sum=1;
long long i=0,d;
long long cnr=nr;
for(i=1;i<=b,v[i]<=sqrt(cnr);i++)
if(nr%v[i]==0)
{
d=1;
while(nr%v[i]==0){d++;nr/=v[i];}
ridic(v[i],d);
dv=(long long)(d*dv)%9973;
if(d>1)
{
cp--;
cp/=(v[i]-1);
sum=(sum%9973)*(cp%9973)%9973;
}
}
if(nr>1)
{
dv=(dv*2)%9973;
sum=((sum%9973)*((nr+1)%9973))%9973;
}
fprintf(g,"%lld ",dv);
fprintf(g,"%lld\n",sum);
}
int main()
{
b=0;
ciur();
fscanf(f,"%d",&a);
for(i=1;i<=a;i++)
{
fscanf(f,"%lld",&nr);
if(nr==1) fprintf(g,"%d %d\n",1,1);
else if(nr<1000000 && prim[nr]==0 )
fprintf(g,"%d% lld\n",2,(nr+1)%9973);
else desc();
}
fclose(g);
return 0;
}