Cod sursa(job #446571)

Utilizator andunhillMacarescu Sebastian andunhill Data 26 aprilie 2010 09:18:56
Problema Suma si numarul divizorilor Scor 60
Compilator cpp Status done
Runda Arhiva educationala Marime 0.91 kb
#include<fstream>
#include<bitset>
#include<math.h>
using namespace std;
#define pmax 1000005
#define max 1000005
#define p 9973
ifstream f("ssnd.in");
ofstream g("ssnd.out");
bitset <max+1> v;
long long n,k,x,nr,sum,q,t,N;
long long prime[max+1];
int pov(long a,long b)
{ long long aux=a;
  for(int i=1;i<b;i++)
	a*=aux;  
  return a;
}

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;
  for(long i=0;1LL*i*i<=N;i++)
  {	if(n%prime[i]==0)
	{  	while(n%prime[i]==0)
		{	q++;
			n/=prime[i];
		}
		nr*=(q+1);
		sum*=(pov(prime[i],q+1)-1)/(prime[i]-1);
		q=0;
	}
  }
  if(n>1)
  { nr*=2;
    sum*=(1LL*n*n-1)/(n-1);
  }
  g<<nr<<" "<<sum%p<<'\n';
}
int main()
{ f>>t;
  ciur();
  for(int i=1;i<=t;i++)
	solve();
  f.close();
  g.close();
  return 0;
}