Pagini recente » Cod sursa (job #564720) | Cod sursa (job #2554909) | Cod sursa (job #1080490) | Cod sursa (job #2575391) | Cod sursa (job #978754)
Cod sursa(job #978754)
#include<fstream>
#include<math.h>
#include<cstring>
#define maxd 1000005
#define maxp 80000
#define mod 9973
using namespace std;
typedef long long ll;
ifstream fin("ssnd.in");
ofstream fout("ssnd.out");
int t,p,nrd,sum;
int prime[maxp];
ll n;
bool sieve[maxd];
void preproc()
{
int i,lim=maxd-5;
prime[++p]=2;
for(i=3;i*i<=lim;i+=2)
if(!sieve[i])
{
prime[++p]=i;
for(int j=i;i*j<=lim;j++)
sieve[i*j]=true;
}
for(;i<=lim;i+=2)
if(!sieve[i])
prime[++p]=i;
}
void gcd(int a,int b,int &x,int &y)
{
if(b==0) {x=1; y=0; return;}
int x0,y0;
gcd(b,a%b,x0,y0);
x=y0%mod;
y=(x0-a/b*y0)%mod;
}
int modular(int a,int b)
{
int x,y;
gcd(a,b,x,y);
while(x<0) x+=mod;
return x;
}
void update(int nr,int prod,int x)
{
nrd*=(nr+1);
sum=(sum*((prod-1)%mod))%mod;
sum=(sum*x)%mod;
}
void solve()
{
int nr,prod,inv;
ll sq;
fin>>t;
for(int i=1;i<=t;i++)
{
fin>>n;
sq=sqrt(n); nrd=1; sum=1;
for(int j=1;prime[j]<=sq;j++)
if(n%prime[j]==0)
{
inv=modular(prime[j]-1,mod);
nr=0; prod=prime[j]%mod;
while(n%prime[j]==0)
{
nr++;
prod=(prod*(prime[j]%mod))%mod;
n/=prime[j];
}
update(nr,prod,inv);
}
if(n>1)
{
inv=modular(n-1,mod);
prod=( (n%mod)*(n%mod) )%mod;
update(1,prod,inv);
}
fout<<nrd<<" "<<sum<<'\n';
}
}
int main()
{
preproc();
solve();
fclose(stdin);
fclose(stdout);
return 0;
}