Cod sursa(job #532038)

Utilizator andreitheo87Teodorescu Andrei-Marius andreitheo87 Data 10 februarie 2011 19:03:40
Problema Suma si numarul divizorilor Scor 0
Compilator cpp Status done
Runda Arhiva educationala Marime 1.21 kb
#include <stdio.h>
#include <vector>
#include <bitset>

using namespace std;

const int MOD = 9973;
int primes[80000];

void calcNdivSdiv(long long n, int* nndiv, long long* ssdiv) {
  int ndiv = 1;
  long long sdiv = 1;
  int div = 2, pos_div = 0;
  int nr;
  long long divPowNr;
  while ((long long) div * div <= n) {
    if (n%div == 0) {
      nr = 1;
      divPowNr = div;
      n/=div;
      while (n%div == 0) {
	++nr;
	divPowNr *= div;
	n /= div;
      }
      if (nr) {
	ndiv *= (nr + 1);
	sdiv *= (divPowNr*div - 1) / (div-1);
      }
    }
    pos_div++;
    div = primes[pos_div];
  }
  if (n != 1) {
    ndiv *= 2;
    sdiv *= (n+1);
  }
  *nndiv = ndiv;
  *ssdiv = sdiv;
}

void compute_primes(int n) {
  int nprimes = 0;
  int prim[1001001];
  for(int i=2; i<=n; i++)
    if(!prim[i]) {
      for(int j=2*i; j<=n; j+=i) prim[j] = true;
      primes[nprimes++] = i;
    }
}

int main() {
  freopen("ssnd.in", "r", stdin);
  freopen("ssnd.out", "w", stdout);
  int nrt;
  scanf("%d", &nrt);
  compute_primes(1001000);
  while (nrt--) {
    long long n;
    scanf("%lld", &n);
    long long sdiv;
    int ndiv;
    calcNdivSdiv(n, &ndiv, &sdiv);
    printf("%d ", ndiv % MOD);
    printf("%d\n", sdiv % MOD);
  }
  return 0;
}