Cod sursa(job #931872)

Utilizator narcis_vsGemene Narcis - Gabriel narcis_vs Data 28 martie 2013 15:41:43
Problema Suma si numarul divizorilor Scor 70
Compilator cpp Status done
Runda Arhiva educationala Marime 1.67 kb
#include<fstream>
#define In "ssnd.in"
#define Out "ssnd.out"
#define N 1000000LL
#define mod 9973
using namespace std;
bool p[N+10];
int prim[78500],nrp=0;
void Ciur()
{
    prim[++nrp]=2;
    long long i,j;
    for(i=3LL;i<=N;i+=2LL)
        if(p[i]==false)
        {
            prim[++nrp] = i;
            for(j=i*i;j<=N;j+=2LL*i)
                p[j] = true;
        }
}

long long Pow_log(int n,int k)
{
    long long p = 1;
    while(k)
    {
        if(k&1)//k este impar
        {
            k--;
            p = ( p * n % mod) % mod;
        }
        k>>=1;
        n = (n %mod *n %mod)%mod;
    }
    return p;
}

int main()
{
    long long T,n,fact,putere,pas,inv;
    long long sum,nrdiv,x;
    ifstream fin(In);
    ofstream fout(Out);
    Ciur();
    fin>>T;
    while(T--)
    {
        fin>>n;
        pas = 1;
        sum = 1;
        nrdiv = 1;
        while(pas<=nrp && prim[pas]*prim[pas]<=n)
        {
            fact = prim[pas]%mod;
            if(n%fact==0)
            {
                putere = 0;
                while(n%fact==0)
                {
                    putere++;
                    n/=fact;
                }
                putere++;
                x = Pow_log(fact,putere)-1;
                inv = Pow_log(fact-1,mod-2);
                x = (x*inv) % mod;
                sum = (sum*x)%mod;
                nrdiv =(nrdiv % mod * putere % mod)%mod;
            }
            pas++;
        }
        if(n>1)
        {
            nrdiv*=2;
            nrdiv%=mod;
            sum = (sum%mod *(n+1)%mod)%mod;
        }
        fout<<nrdiv<<" "<<sum<<"\n";
    }
    return 0;
}