Cod sursa(job #446579)

Utilizator andunhillMacarescu Sebastian andunhill Data 26 aprilie 2010 09:38:41
Problema Suma si numarul divizorilor Scor 0
Compilator cpp Status done
Runda Arhiva educationala Marime 1.11 kb
#include<fstream>
#include<bitset>
using namespace std;
#define pmax 1000005
#define max 1000005
const int p=9973;
ifstream f("ssnd.in");
ofstream g("ssnd.out");
bitset <max+1> v;
int k,x,nr,sum,q,t;
long long n,N;
int prime[max+1];
inline int pow(int x, int y) 
{ int rez = 1;    
	x %=p; 
	for(; y; y >>= 1)     
	{  
	if(y & 1)      
		rez *= x,    
		rez %= p;    
		x *= x;  
		x %= p;     
	} 
	return rez; 
} 
void ciur()
{ for(long i=2;i<=pmax;i++)
  	if(v[i]==0)
	{	prime[k]=i;
		for(long j=(i<<1);j<pmax;j+=i)
			v[j]=1;
		k++;
	}
}
void solve()
{ f>>n;
  x=0,q=0,nr=1,sum=1,N=n;
  long long s=1;
  for(long i=0;1LL*prime[i]*prime[i]<=N;i++)
  {	if(n%prime[i]==0)
	{  	while(n%prime[i]==0)
		{	q++;
			s*=prime[i];
			n/=prime[i];
		}
		nr*=(q+1);
		sum*=((1LL*s*prime[i]-1)%p);
		sum%=p;
		sum*=pow(prime[i]-1,p-2); 
		sum%=p;
		q=0;
	}
  }
  if(n>1)
  { nr*=2;
    sum*=((1LL*n*n-1)%p);
	sum%=p;
	sum*=pow(n-1,p-2);
	sum%=p;
  }
  g<<nr<<" "<<sum<<'\n';
}
int main()
{ f>>t;
  ciur();
  for(int i=1;i<=t;i++)
	solve();
  f.close();
  g.close();
  return 0;
}