Cod sursa(job #879636)

Utilizator roots4Irimia Alexandru Gabriel roots4 Data 15 februarie 2013 18:19:35
Problema Suma si numarul divizorilor Scor 70
Compilator cpp Status done
Runda Arhiva educationala Marime 1.08 kb
#include<fstream>
#include<string.h>
#include<math.h>
#define mod 9973
#define max_n 100010
#define max_s 80
using namespace std;

ifstream f("ssnd.in");
ofstream g("ssnd.out");

int T,i,d,k,j;
char V[max_n];
unsigned long long nr,p1,p2,nr_d,s_d;
int Prim[max_n],Nr_d;
double rad;

long long power(int x , int y){
	
	long long p=x,sol=1;
	
	while(y){
		if(y%2==1)
			sol*=p;
		p*=p;;
		y>>=1;
	}
	
	return sol;
}


void ciur(){
	
	Prim[++k]=2;
	
	for(i=3 ; i<=max_n ; i+=2)
		if(V[i]==0){
			Prim[++k]=i;
			for(j=2*i ; j<=max_n ; j+=i)
				V[i]=1;
		}
	
}

int main(){
	
	f>>T;
	
	ciur();
	
	while(T--){
		f>>nr;
		
		nr_d=1;s_d=1;
		
		rad=sqrt(nr);
		
		for(i=1 ; i<=k && Prim[i]<=rad && nr!=1; i++)
			if(nr%Prim[i]==0){
				Nr_d=0;
				while(nr%Prim[i]==0){
					Nr_d++;
					nr/=Prim[i];
				}
				nr_d*=(Nr_d+1);
				s_d=(  (s_d%mod) * (  ( power(Prim[i],Nr_d+1)-1 )/ (Prim[i]-1)  )  ) % mod;
			}
		
		if(nr!=1){
			nr_d*=2;
			s_d=((s_d%mod)*((nr*nr)-1)/(nr-1))%mod;
		}
		
		g<<nr_d<<" "<<s_d<<"\n";
	}
	
	
	return 0;
}