Найти произведение элементов вектора удовлетворяющих неравенству

#include <iostream>
 
int main()
{
    const int N = 10, M = 5;
    int Y[N], mult = 1, counter = 0;
    // инициализируем значениями (помним, что индексация массивов начинается с нуля)
    for (int i = 1; i <= M; i++)
        Y[i-1] = 1 + i*i;
    for (int i = M; i <= 2*M; i++)
        Y[i-1] = i + Y[i-2];
 
    // выводим результат, заодно находим произведение
    for (int i = 0; i < N; i++)
    {
        std::cout << Y[i] << " ";
        // если удобвлетворяет условиям
        if (Y[i] > 0 && Y[i] <= 5)
        {
            mult *= Y[i]; // умножаем результат на значение
            counter++;    // увеличиваем счётчик на 1
        }
    }
 
    std::cout << "\nMult: " << mult << " with counter " << counter << std::endl;
}

// пример обработки - подсчет суммы элементов
int Work(int* M,int n)
{
    int sum=1;
    for(int i=0;i<n;i++)
    {
        if(M[i]>0 && M[i]<5)
          sum*=M[i];
    }
    return sum;
}

void RandomFill(int* M,int n) // можете оставить имя для простоты
{
    for(int i=0;i<n;i++)
    {
        if(i<=5)
           M[i]=1+i*i;
        else
          M[i]=M[i-1]+i;
    }
}