Найти номера трёх точек, которые являются вершинами треугольника наибольшей площади

#include <iostream>
#include <random>
#include <iterator>
#include <math.h>
 
int main(void)
{
    size_t n;
    std::cin >> n;
    double* x = new double[n];
    double* y = new double[n];
    std::uniform_real_distribution<double> distribution(-10, 10);
    std::default_random_engine generator;
    std::generate(x, x + n, [&](){ return distribution(generator); });
    std::generate(y, y + n, [&](){ return distribution(generator); });
 
    for (size_t i(0); i < n; ++i) std::cout << "(" << x[i] << ", " << y[i] << ")" << (i == n - 1 ? "\n\n" : ", ");
 
    double S(0);
    size_t I, J, K;
    for (size_t i(0); i < n - 2; ++i)
        for (size_t j(i + 1); j < n - 1; ++j)
            for (size_t k(j + 1); k < n; ++k) {
                double a, b, c, p, s;
                a = sqrt(pow(x[i], 2) + pow(y[i], 2));
                b = sqrt(pow(x[j], 2) + pow(y[j], 2));
                c = sqrt(pow(x[k], 2) + pow(y[k], 2));
                p = 0.5 * (a + b + c);
                s = sqrt(p * (p - a) * (p - b) * (p - c));
                if (s > S) S = s, I = i, J = j, K = k;
            }
    std::cout << "S: " << S << '\n' << "p1: " << ++I << '\n' << "p2: " << ++J << '\n' << "p3: " << ++K;
 
    delete[] x;
    delete[] y;
    return 0;
}

#include <iostream>
#include <locale.h>
#include <iomanip>
#include <time.h>
#include <math.h>
#include <stdio.h>
#include <stdlib.h>
using namespace std;
 
struct Point
{
    double x, y;
};
 
const int nmax = 20;
 
double x, y;
Point trg[3];
 
double Plsh(Point a, Point b, Point c)
{
    double f;
    f = abs(((a.x - c.x) * (b.y - c.y) - (b.x - c.x) * (a.y - c.y)) / 2);
    return f;
}
 
int main()
{
    setlocale(LC_ALL, "Russian");
 
    srand(time(NULL));
 
    int n, i, j, k, imx, jmx, kmx;
    double mx;
    Point a[nmax];
    rand();
 
    do {
        cout << "Количество точек от 3 до " << nmax << " n " << endl;
        cin >> n;
    } while (n < 3 || n > nmax);
    for (i = 1; i < n; i++)
    {
        a[i].x = (double)(rand() % 10000 - 5000) / 100;
        a[i].y = (double)(rand() % 10000 - 5000) / 100;
    }
    cout << "Координаты верщин треугольников: " << "  " << endl;
    for (i = 1; i < n; i++)
    {
        cout << i << endl;
        cout << "x: " << endl;
    }
    for (i = 1; i < n; i++)
        cout << a[i].x;
    cout<< endl;
    cout << "y: " << endl;
    for (i = 1; i < n; i++)
        cout << a[i].y << endl;
    mx = Plsh(a[1], a[2], a[3]);
    imx = 1;
    jmx = 2;
    kmx = 3;
    for (i = 1; i < n - 2; i++)
        for (j = i + 1; j < n - 1; j++)
            for (k = j + 1; k < n; k++)
                if (Plsh(a[i], a[j], a[k]) > Plsh(a[imx], a[jmx], a[kmx]))
                {
                    mx = Plsh(a[i], a[j], a[k]);
                    imx = i;
                    jmx = j;
                    kmx = k;
                };
    cout << "Максимальная площадь= " << mx << " у треугольника с координатами" << endl;
    cout << imx << "( " << a[imx].x << " ; " << a[imx].y << " )" << endl;
    cout << jmx << "( " << a[jmx].x << " ; " << a[jmx].y << " )" << endl;
    cout << kmx << "( " << a[kmx].x << " ; " << a[kmx].y << " )" << endl;
    system("pause");
}