Среднее геометрическое для n простых чисел на c++

#include <iostream>
#include <math.h>
 
 
unsigned int prostoe(unsigned int N) //Описание функции.
{
    //Функция определяет, является ли число простым.
    int i, pr;
    for (pr = 1, i = 2; i <= N / 2; i++)
        if (N % i == 0) { pr = 0; break; }
    return pr;
}
 
using namespace std;
int main()
{
    setlocale(LC_CTYPE, "russian");
 
    unsigned int i, N, X, S, kp, ks;
    long int P;
    cout << "N="; cin >> N;
    for (kp = ks = S = 0, P = 1, i = 1; i <= N; i++)
    {
        cout << "X="; cin >> X; //Вводится элемент последовательности.
        if (prostoe(X)) // X — простое число.
        {
            kp++; //Счётчик простых чисел.
            P *= X; //Произведение простых чисел.
        }
        
    }
    if (kp > 0) //Если счётчик простых чисел больше нуля,
        //считаем среднее геометрическое и выводим его,
        cout << "Среднее геометрическое=" << pow(P, (float)1 / kp) << endl;
    else //в противном случае –– сообщение об отсутствии простых чисел.
        cout << "Нет простых чисел\n";
    
    return 0;
}

#include <iostream>
#include <vector>
#include <cmath>
 
void get_primes(std::vector<int> &dest, int n)
{
    dest.clear();
    std::vector<bool> sieve(n + 1, true);
    for(uint64_t i = 2; i <= n; ++i)
        if(sieve[i])
        {
            dest.push_back(i);
            for(uint64_t j = i * i; j <= n; j += i)
                sieve[j] = false;
        }
}
 
int main()
{
    int n = 50'000'000;
    std::vector<int> primes;
    get_primes(primes, n);
 
    double res { 1. };
    double pwr { 1. / primes.size() };
    for(auto i : primes)
        res *= pow(i, pwr);
 
    std::cout << res;
 
    return 0;
}