Cod sursa(job #862418)

Utilizator lehman97Dimulescu David lehman97 Data 22 ianuarie 2013 18:09:58
Problema Suma si numarul divizorilor Scor 50
Compilator cpp Status done
Runda Arhiva educationala Marime 1.31 kb
#include <stdio.h>
#include <math.h>


using namespace std;

FILE *f=fopen("ssnd.in","r");
FILE *g=fopen("ssnd.out","w");

const int n=1000000;
bool prim[n];
long long  x,nr2;
unsigned long long  nr,i,a,p,cnt=0;




void ciur()
{
    int i,j;



    for (i=2;i<=n;i++)
        if (!prim[i])
        {
            cnt++;
            for (j = 2*i;j<=n;j+=i)
                prim[j] = 1;
        }
}







void desc()
{
    long long x,dv=1,sum=1;
    long long  i=0,d;
    long long  cnr=nr;

    for(i=2;i<=sqrt(cnr);i++)
    if(!prim[i])
    {

        d=1;
        x=i;
        while(nr%i==0){d++;nr/=i;x*=i;}
        dv=(long long)(d*dv)%9973;
        if(d>1)
        {

            x--;
            x/=(i-1);
            sum=(sum%9973)*(x%9973)%9973;
        }


    }

    if(nr>1)
    {
        dv=(dv*2)%9973;
        sum=((sum%9973)*((nr+1)%9973))%9973;
        }



    fprintf(g,"%lld ",dv);
    fprintf(g,"%lld\n",sum);



}



int main()
{
    ciur();
    fscanf(f,"%d",&a);
    for(i=1;i<=a;i++)
    {
         fscanf(f,"%lld",&nr);
         if(nr==1) fprintf(g,"%d%d\n",1,1);
         else if(nr<1000000 && prim[nr]==0 )
              fprintf(g,"%d%lld\n",2,(nr+1)%9973);
         else desc();
    }


    fclose(g);
    return 0;
}