Cod sursa(job #1311889)

Utilizator pitradaPit-Rada Ionel-Vasile pitrada Data 8 ianuarie 2015 18:19:56
Problema Suma si numarul divizorilor Scor 70
Compilator cpp Status done
Runda Arhiva educationala Marime 1.35 kb
#include<fstream>
using namespace std;
ifstream fin ("ssnd.in");
ofstream fout("ssnd.out");

long long i,j,a,s,x,nrd,p[80001],pj;
int t,np,c, z[1000009];

long long pow(long long base, int expo)
{
    long long p=1, a=base;
    while(expo)
    {
        if(expo&1)
        {
            p*=a;
        }
        a*=a;
        expo=(expo>>1);
    }
    return p;
}

int main()
{
    np=0;
    for(i=2;i<=1000;i++)
    {
        if(z[i]==0)
        {
            np++;
            p[np]=i;
            for(j=i*i;j<=1000000;j=j+i)
            {
                z[j]=1;
            }
        }
    }
    for(i=1001;i<=1000000;i++)
    {
        if(z[i]==0)
        {
            np++;
            p[np]=i;
        }
    }
    fin>>t;
    for(i=1;i<=t;i++)
    {
        fin>>a;
        s=1;
        nrd=1;
        j=1;
        pj=p[j];
        while(pj*pj<=a)
        {
            c=0;
            while(a%pj==0)
            {
                c++;
                a/=pj;
            }
            x=pow(pj,c+1);
            nrd*=(1+c);
            s*=((x-1)/(pj-1));
            //s%=9973;
            j++;
            pj=p[j];
        }
        if(a>1)
        {
            nrd=(nrd<<1);
            s*=(a+1);
            //s%=9973;
        }
        fout<<nrd<<" "<<s%9973<<"\n";
    }
    return 0;
}