Cod sursa(job #2347055)

Utilizator mateibanuBanu Matei Costin mateibanu Data 18 februarie 2019 12:52:44
Problema Suma si numarul divizorilor Scor 40
Compilator cpp-64 Status done
Runda Arhiva educationala Marime 1.35 kb
#include <bits/stdc++.h>

using namespace std;

#define MOD 9973

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

int inv(int x, int n){
    if (n==0) return 1;
    if (n==1) return x;
    int 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,x,nr=0,j;
    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("%d",&x);
        k=0;
        for (i=1;i<=nr&&x>1;i++){
            if (x%prime[i]==0){
                p[++k]=prime[i];
                d[k]=0;
                while (x%prime[i]==0){
                    x/=prime[i];
                    d[k]++;
                }
            }
        }
        if (x>1){
            p[++k]=x;
            d[k]=1;
        }
        sus=1;
        jos=1;
        dv=1;
        for (i=1;i<=k;i++){
            dv=(dv*(d[i]+1))%MOD;
            sus=(sus*((inv(p[i],d[i]+1)-1+MOD)%MOD))%MOD;
            jos=(jos*((p[i]-1+MOD)%MOD))%MOD;
        }
        jos=inv(jos,MOD-2);
        sus=(sus*jos)%MOD;
        printf("%d %d\n",dv,sus);
    }
    return 0;
}