Cod sursa(job #1166951)

Utilizator alevasluialeHuhurez Marius alevasluiale Data 3 aprilie 2014 23:55:36
Problema Suma si numarul divizorilor Scor 10
Compilator cpp Status done
Runda Arhiva educationala Marime 1.19 kb
#include <cstdio>
#include <bitset>
using namespace std;
const int N=1000002;
bitset <N> v;
int p[78500];
const int mod=9973;
long long putere(long long a,long long p)
{
   long long result=1;
    while(p)
    {
        if(p&1) result=(result*a)%mod;
        a=(a*a)%mod;
        p>>=1;
    }
    return result;
}
void ciur()
{
    int i,j;
    p[++p[0]]=2;
    for(i=3;i*i<N;i+=2)
    {   if(!v[i])
        for(j=i*i;j<N;j+=i) v[j]=1;
    }
    for(i=3;i<N;i+=2) if(!v[i]) p[++p[0]]=i;
}
int main()
{   freopen("ssnd.in","r",stdin);
    freopen("ssnd.out","w",stdout);
   long long n,t,i,j,q,x,nr,ok,a,b;
    long long sum;
    ciur();
    scanf("%d",&t);
    for(i=0;i<t;i++)
    {
        sum=1;
        nr=1;
        scanf("%lld",&n);
        ok=n;
        for(j=1;n>1||p[j]*p[j]<=ok;j++)
        {   q=0;
            x=n;
            a=1;
            while(n%p[j]==0) {q++;a=(a*p[j])%mod;n/=p[j];}
            if(q>0)
            {
            nr*=(q+1);
            a=putere(p[j],q+1)-1;
            b=(a/(p[j]-1))%mod;
            sum*=b%mod;
            }
        }
        if(n>1) {nr*=2;sum+=n;}
        printf("%lld %lld\n",nr,sum%mod);
    }
}