Cod sursa(job #1851753)

Utilizator Radu_FilipescuFilipescu Radu Radu_Filipescu Data 20 ianuarie 2017 00:33:21
Problema Suma si numarul divizorilor Scor 70
Compilator cpp Status done
Runda Arhiva educationala Marime 1.36 kb
#include <fstream>
#include <cmath>

#define upper 1000001

using namespace std;

ifstream fin("ssnd.in");
ofstream fout("ssnd.out");

bool prim[upper];
int prime[100001];

void Eratostene()
{
   int i,j;

   prim[2]=true;
   for(i=3; i<=upper; i=i+2) prim[i]=true;
   for(i=3; i*i<=upper; i=i+2)
    if(prim[i]==true)
    for(j=3*i; j<=upper; j=j+2*i) prim[j]=false;

}

int main()
{
  int t,i,j,k,exp,nrprime=0;
  long long n,aux,nrdiv=1,sumdiv=1,sum;

  Eratostene();

  for(i=2; i<=upper; i++)
    if(prim[i]==true) {
                        nrprime++;
                        prime[nrprime]=i;
                      }

  fin>>t;
  for(i=1; i<=t; i++)
  {
     fin>>n;

     if(n<upper && prim[n]==true) fout<<'2'<<' '<<(n+1)%9973<<'\n';
     else
     {
       nrdiv=1;
       sumdiv=1;

       j=1;
       while(n>1)
         if(nrdiv==1 && j*j>n) fout<<'2'<<' '<<(n+1)%9973<<'\n';
         else if(n%prime[j]==0)
       {
         exp=1;
         while(n%prime[j]==0) { exp++; n=n/prime[j]; }

         //fout<<j<<' '<<exp<<'\n';
         nrdiv=nrdiv*exp;

         aux=1;
         for(k=1; k<=exp; k++) aux=aux*prime[j];
         //fout<<aux<<'\n';

         sumdiv=sumdiv*((aux-1)/(prime[j]-1));
         j++;
       }
         else j++;

       fout<<nrdiv<<' '<<sumdiv%9973<<'\n';
     }
  }


    return 0;
}