Cod sursa(job #654273)

Utilizator ion824Ion Ureche ion824 Data 29 decembrie 2011 23:38:59
Problema Suma si numarul divizorilor Scor 70
Compilator cpp Status done
Runda Arhiva educationala Marime 1.5 kb
#include<cstdio>
#include<cmath>
#define MOD 9973
bool b[1000010];
int pr[1000010],K; 

void ciur(int n){
     for(int i=2;i<=n;++i)
       if(!b[i]){
         pr[++K]=i;
         for(int j=2;j<=n/i;++j)
           b[i*j]=1; 
           }            
     }
   
/*int pow1(int a,int n){
    int rez=a%MOD;
      for(int i=1;i<n;++i){
        rez*=a;
        rez%=MOD;
        }   
    return rez;    
}     */

inline int pow1(int x, int p){
	int rez = 1; x %= MOD;
	for(; p; p >>= 1) {
		if(p & 1) {
			rez *= x;
			rez %= MOD;
		}
		x *= x;
		x %= MOD;
	}
	return rez;
}

int main(void){
    FILE *fin = fopen("ssnd.in","r");
    FILE *fout = fopen("ssnd.out","w");
    int i,t,sumdiv,nrdiv,d,p,k,p2; long long n=0; 
    fscanf(fin,"%d",&t); 
    ciur(1000005);
    while(t--){
        fscanf(fin,"%llu",&n);      
        nrdiv=1; sumdiv=1; i=1; k=int(sqrt(n));  
        while(n!=1 && i<=K){
                if(n%pr[i]==0){    
                    d=0;  
                    while(n%pr[i]==0 && n!=1){ ++d; n/=pr[i];  }
                    nrdiv*=(d+1); 
                    p=(pow1(pr[i],d+1)-1)%MOD;
                    p2 = pow1(pr[i]-1, MOD-2) % MOD;
                     sumdiv=(((sumdiv*p)%MOD)*p2)%MOD;
                    ++i;
                    }
                else ++i;                         
                }
          if(n>1){ nrdiv*=2; sumdiv=((n+1)*sumdiv)%MOD; }      
        fprintf(fout,"%d %d\n",nrdiv,sumdiv);      
      }     
 return 0;   
}