Cod sursa(job #2863580)

Utilizator BlueLuca888Girbovan Robert Luca BlueLuca888 Data 6 martie 2022 21:46:15
Problema Suma si numarul divizorilor Scor 100
Compilator cpp-64 Status done
Runda Arhiva educationala Marime 1.6 kb
#include <bits/stdc++.h>
#define ull unsigned long long

using namespace std;

ifstream fin  ("ssnd.in");
ofstream fout ("ssnd.out");

const ull MOD = 9973;
const ull LIM = 1e6;
bitset <LIM+5> ciur;
ull p, prim[80000], expo;

void mark(ull num){
    for(ull i=num*num; i<=LIM; i+=num)
        ciur[i] = true;
}

ull n, dcnt, dsum;

ull lgput(ull a, ull b){
    ull answer = 1;
    while(b != 0){
        if(b&1)
            answer *= a;
        a *= a;
        b >>= 1;
    }
    return answer;
}

signed main (){
    ios_base::sync_with_stdio(false);
    fin.tie(nullptr);
    fout.tie(nullptr);

    ciur[0] = ciur[1] = true;
    mark(2);
    mark(3);
    for(ull i=5; i<=LIM/i; i+=6){
        if(!ciur[i]) mark(i);
        if(!ciur[i+2]) mark(i+2);
    }
    prim[++p] = 2;
    prim[++p] = 3;
    for(ull i=5; i<=LIM; i+=6){
        if(!ciur[i]) prim[++p] = i;
        if(!ciur[i+2] && i+2<=LIM) prim[++p] = i+2;
    }

    ull teste;
    fin>>teste;
    while(teste--){
        fin>>n;
        dcnt = dsum = 1;

        for(ull i=1; i<=p; i++){
            if(prim[i]*prim[i] > n) break;
            if(n % prim[i] != 0) continue;

            expo = 0;
            while(n % prim[i] == 0){
                expo++;
                n /= prim[i];
            }

            dcnt *= (expo+1);
            dsum = (ull)dsum * ((lgput(prim[i], expo+1) - 1) / (prim[i] - 1)) % MOD;
        }
        if(n > 1){
            dcnt *= 2;
            dsum = (ull)dsum * ((n*n-1) / (n-1)) % MOD;
        }
        fout<<dcnt<<" "<<dsum<<"\n";
    }
    return 0;
}