Сформировать массив по указанному правилу

#include <bits/stdc++.h>
 
using namespace std;
 
int main() {
  int n = 20;
 
  vector<double> a(n);
  for (int i = 1; i <= n; ++i)
    a[i - 1] = abs(i * i * i + 3) / i - i * i + 15 * i;
 
  copy(a.cbegin(), a.cend(), ostream_iterator<double>(cout, " "));
 
  for (int i = 0; i < n; ++i) {
    if (a[i] <= 25 || int(a[i]) % 3 != 0) {
      a.erase(a.cbegin() + i);
 
      --i;
      --n;
    }
  }
 
  cout << '\n';
 
  copy(a.cbegin(), a.cend(), ostream_iterator<double>(cout, " "));
}

srand((int)time(0));
    const int n = 20;
    double a[n];
    int menu, sum = 0;
cout << "\nFilling the array with pseudo-random numbers" << endl;
        cout << "a[n] = ";
        for (int i = 0; i < n; i++)
        {
            a[i] = rand() % 99;
            cout << a[i] << ' ';
        }
 
        for (int i = 0; i < n; i++)
        {
            sum += a[i];
        }
        cout << "\nArithmetic average of array elements: " << sum / n << endl;
 
        int min_index = 0, min_elem = a[0];
        for (int i = 0; i < n; i++)
        {
            if (min_elem > a[i])
            {
                min_elem = a[i];
                min_index = i;
            }
        }
        cout << "The minimum value = " << min_elem << " and index = a[" << min_index << "]" << '\n';
        cout << endl;
        cout << "Array of real numbers: ";
        for (int i = 0; i < n; i++)
        {
            a[i] = abs(pow(i, 3) + 3) / 3 - pow(i, 2) + 15 * i;
            cout << a[i] << ' ';
        }
        cout << endl;
        cout << "\nConvert an array, leaving only elements larger than 25 and multiple 3" << endl;
        for (int i = 0; i < n; i++)
        {
            if ((a[i] > 25) && ((int)a[i] % 3 == 0))
            {
                cout << a[i] << ' ';
            }
        }
        cout << endl;