Cod sursa(job #548791)

Utilizator ProcopliucProcopliuc Adrian Procopliuc Data 7 martie 2011 20:06:33
Problema Suma si numarul divizorilor Scor 70
Compilator cpp Status done
Runda Arhiva educationala Marime 1.14 kb
# include <fstream.h>
# include <math.h>
# define mod 9973
ifstream f ("ssnd.in");
ofstream g ("ssnd.out");
long long int m=1000005,sd,pd,i,n,p,q,w,t,j,k;
int a[1000000];
bool v[1000005];
  
  void ciur ()
  {
	  long long int i,j;
	  for (i=2;i<=m/2;i++)
		  if (v[i]==0)
		  {
			  k++;
			  a[k]=i;
			  for (j=i*i;j<=m;j=j+i)
				  v[j]=1;
		  }
  }
  
  int putere (int a,int b)
  {
	  if (b==1)
		  return a%mod;
	  else
		  if (b%2==0)
			  return (int) pow (putere (a,b/2),2)%mod;
		  else
			  return a*putere (a,b-1)%mod;
  }

int main ()
{
	ciur ();
	f>>t;
	
	for (;t>=1;t--)
	{
		f>>n;
		sd=1;
		pd=1;
		for (j=1;j<=k && a[j]<=sqrt (n);j++)
		{
			i=a[j];
			p=0;
			if (n%i==0)
			{	while (n%i==0)
				{
					n=n/i;
					p++;
				}
				
			sd=sd*(p+1);
			
			q=(putere (i%mod,p+1)%mod)-1;
			w=(putere ((i-1)%mod,mod-2))%mod;
			
			pd=(pd*((q*w)%mod))%mod;
		
			}
		}
		if (n!=1)
		{
			p=1;
			i=n;
		    sd=sd*(p+1);
			
			q=(putere (i%mod,p+1)%mod)-1;
			w=(putere ((i-1)%mod,mod-2))%mod;
			
			pd=(pd*((q*w)%mod))%mod;
		}
		
		g<<sd<<" "<<pd<<"\n";
	}
	
	return 0;
}