Cod sursa(job #769067)

Utilizator mi5humihai draghici mi5hu Data 18 iulie 2012 10:58:43
Problema Suma si numarul divizorilor Scor 70
Compilator cpp Status done
Runda Arhiva educationala Marime 1.7 kb
#include <stdio.h>
#include <vector>
#define int64 long long int
#define MOD 9973
#define NMAX_P 1000001

using namespace std;

int64 PR[100000];
int64 PR_S;
bool bifat[NMAX_P];

void print_(int64 S, int64 ND) {
     printf("%lld %lld\n", ND, S);          
}

int64 my_pow(int64 b, int64 e) {
     if (e == 0) {
         return 1;   
     } else if (e == 1) {
         return b;       
     } 
     
     int64 aux = my_pow(b, e/2);
     if (e % 2 == 1) {
         return (aux * aux * b);   
     } else {
         return aux * aux;       
     }
}

void solve_(int64 nr, int64 &S, int64 &ND) {
     int64 DIV_NR;
     for (int i = 0; i < PR_S && nr > 1; i++) {
         if (nr % PR[i] == 0) {
            DIV_NR = 0;
            while (nr % PR[i] == 0) {
                DIV_NR++;
                nr = nr / PR[i];      
            }
            ND *=  (DIV_NR + 1);
            S = (S * (my_pow(PR[i], DIV_NR + 1) - 1) / (PR[i] - 1)) % MOD;
         }    
     }          
} 

void read_() {
     int n;
     int64 nr;
     int64 S, P;
     
     scanf("%d", &n);
     for (int i = 0; i < n; i++) {
         scanf("%lld", &nr);
         S = 1; P = 1;
         solve_(nr, S, P);
         print_(S, P);
     }     
}

void marcheaza(int64 i) {
      for (int64 j = i; j <= NMAX_P; j += i) {
          bifat[j] = true;
      }    
}

void getPrimes() {
     for (int64 i = 2; i <= NMAX_P; i++) {
         if (bifat[i] == false) {
             marcheaza(i);             
             PR[PR_S++] = i;
         }
     }
}

int main() {
    freopen("ssnd.in", "r", stdin);
    freopen("ssnd.out", "w", stdout);
    
    getPrimes();
    read_();
    return 0;
}