Cod sursa(job #2347082)

Utilizator mateibanuBanu Matei Costin mateibanu Data 18 februarie 2019 13:27:29
Problema Suma si numarul divizorilor Scor 70
Compilator cpp-64 Status done
Runda Arhiva educationala Marime 1.46 kb
#include <bits/stdc++.h>

using namespace std;

#define MOD 9973
#define ll long long

int prime[100010],p[100010],d[100010],ciur[1000010];
ll sus, jos, dv;

ll inv(ll x, ll n){
    if (n==0) return 1;
    if (n==1) return x;
    ll y=inv(x,n/2);
    y=(y*y)%MOD;
    if (n%2==1) y=(y*x)%MOD;
    return y;
}

int main()
{
    int i,n,k,nr=0,j;
    ll x;
    freopen("ssnd.in","r",stdin);
    freopen("ssnd.out","w",stdout);
    for (i=2;i<=1000000;i++){
        if (ciur[i]==0){
            prime[++nr]=i;
            for (j=2*i;j<=1000000;j+=i)
                ciur[j]=1;
        }
    }
    scanf("%d",&n);
    while (n){
        n--;
        scanf("%lld",&x);
        k=0;
        sus=1;
        jos=1;
        dv=1;
        for (i=1;1LL*prime[i]*prime[i]<=x&&x>1&&i<=nr;i++){
            if (x%prime[i]==0){
                p[++k]=prime[i];
                d[k]=0;
                while (x%prime[i]==0){
                    x/=prime[i];
                    d[k]++;
                }
                dv=(dv*(d[k]+1))%MOD;
                sus=(sus*((inv(p[k],d[k]+1)-1+MOD)%MOD))%MOD;
                jos=(jos*((p[k]-1+MOD)%MOD))%MOD;
            }
        }

        if (x>1){
            dv=(dv*2)%MOD;
            sus=(sus*((inv(x,2)-1+MOD)%MOD))%MOD;
            jos=(jos*((x-1+MOD)%MOD))%MOD;
        }
        jos=inv(jos,MOD-2);
        sus=(sus*jos)%MOD;
        printf("%lld %lld\n",dv,sus);
    }
    return 0;
}