Cod sursa(job #1366401)

Utilizator sulzandreiandrei sulzandrei Data 1 martie 2015 00:38:58
Problema Suma si numarul divizorilor Scor 0
Compilator cpp Status done
Runda Arhiva educationala Marime 1.62 kb
#include <iostream>
#include <fstream>
#include <bitset>
#include <cmath>
using namespace std;
ifstream in("ssnd.in");
ofstream out("ssnd.out");
#define dim 1000001
#define mod 9973
#define ull unsigned long long int
bitset<dim> v;
ull n,i,j,t,cop,Ndiv,nr,Sdiv,tes;
ull p[1000000],d[1000000];
void Sieve()
{
    for(i=2;i<=dim;i++)
        v[i]=1;
    for(i=2;i*i<=dim;i++)
        if(v[i])
            for(j=i*i;j<=dim;j+=i)
                v[j]=0;
}
void precalcul()
{
    nr=0;
    cop = n;
    for(i=2;i*i<=n;i++)
        if ((v[i]==1) && (n%i==0))
        {
            p[nr++] = i;
            d[nr-1] = 1;
            while(cop%i==0)
            {
                cop /=i;
                d[nr-1]++;
            }
        }
}
ull D()
{
    if (cop==n)
    {
        if (cop==1)
            return 1;
        else
            return 2;
    }
    Ndiv = 1;
    for(i=0;i<nr;i++)
        Ndiv *=d[i];
    return Ndiv;
}
ull put(ull base,ull k)
{
    ull b=base,res=1;
    while(k)
    {
        if(k&1)
            res =((res%mod)*(base%mod))%mod;
        base %=mod;
        base *=base;
        base %=mod;
        k>>=1;
    }
    res = (res-1)/(b-1);
    return res;
}
ull S()
{
    Sdiv = 1;
    if (nr==0)
    {
        if (n==1)
            return 1;
        else
            return n+1;
    }
    for(i=0;i<nr;i++)
        Sdiv = ((Sdiv%mod)*(put(p[i],d[i])%mod))%mod;
    return Sdiv%mod;
}
int main()
{
    in>>t;
    Sieve();
    for(tes=0;tes<t;tes++)
    {
        in>>n;
        precalcul();
        out<<D()<<" "<<S()<<endl;
    }
    return 0;
}