Выбрать 3 точки множества, составляющих треугольник наибольшего периметра

/*
Выбрать 3 различные точки заданного на плоскости множества точек,
составляющих треугольник наибольшего периметра.
 
Пусть вершины A(xA; yA; zA), B(xB; yB; zB), C(xC; yC; zC). 
Тогда периметр треугольника найдется как сумма длин сторон ВА, BC, CA:
P = |ВА| + |BC| + |CA| =корень((xA – xB)^2 + (yA – yB)^2 + (zA – zB)^2) +
+ корень( (xC – xB)^2  + (yC – yB)^2  + (zC – zB)^2 ) + 
+ корень( (xA – xC)^2  + (yA – yC)^2  + (zA – zC)^2 );
Пусть множество заданно необходимым и достаточным образом, 
то есть любой элемент множества встречаеться в описании один 
и только один раз, а так же существует правило по которому мы можем
перебрать все элементы множества не повторяясь.
Тогда, решение можно найти полным перебором всех возможных
комбинаций A, B, C 
В частности, если множество задано массивом элементами которого являются
координаты:
(ужжос как пафосно... всё намного проще =) )
*/
 
#include<math.h>
//только для использования sqrt;
 
#include<iostream>
using namespace std;
 
double d(double *p1, double *p2)// возвращает расстояние между точками p1 p2
{
    double dx2,dy2,dz2; //dx2 = (Xa - Xb)^2
    dx2 = (p1[0]-p2[0])*(p1[0]-p2[0]);
    dy2 = (p1[1]-p2[1])*(p1[1]-p2[1]);
    dz2 = (p1[2]-p2[2])*(p1[2]-p2[2]);
    return sqrt(dx2+dy2+dz2);
}
 
int main()
{
    //инициализация массива случайной размерности:
    srand( time(NULL) );
    int dimm = rand()%100+3;//размерность не меньше трёх.
    cout<<"\ndimm = "<<dimm<<endl;
    double** m = new double* [dimm];
    for (int i = 0;i<dimm;i++) m[i] = new double [3];
    
    //Заполним его случайными координатами в диапазоне от -1 до 1:
    for (int i = 0;i<dimm;i++)
    {
        m[i][0] = double(rand()%200 - 100)/100; //x from -1 to +1;
        m[i][1] = double(rand()%200 - 100)/100; //y from -1 to +1;
        m[i][2] = double(rand()%200 - 100)/100; //z from -1 to +1;          
    }
    // ниже идёт описание вспомогательных переменных:
    //p = perimetr ab,bc,ca = distance; a,b,c = points;
    double p = 0;
    double ab,bc,ca; ab = bc = ca = 0;
    int a,b,c; a = b = c = 0;
    
    // Все приготовления завершены, можно приступить к поиску решения:
    for(int ia=0;ia<dimm;ia++)
    {
        for(int ib=0;ib<dimm;ib++)
        {
            for(int ic=0;ic<dimm;ic++)
            {
                ab = d(m[ia],m[ib]);
                bc = d(m[ib],m[ic]);
                ca = d(m[ic],m[ia]);
                if(ab+bc+ca > p) // если верно, то найден новый максимальный периметр
                {
                    p = ab+bc+ca;
                    cout<<"new max p = "<<p<<endl;//можно закомментировать строку.
                    a = ia;
                    b = ib;
                    c = ic;
                }
            }
        }
    }
    // Решение найдено, осталось только вывести его куда-нибудь.
    
    cout<<"\n\nmax P = "<<p<<endl;
    cout<<"\t"<<"\tx"<<"\t\ty"<<"\t\tz"<<endl;
    cout<<" points A\t"<<m[a][0]<<"\t\t"<<m[a][1]<<"\t\t"<<m[a][2]<<"\t\t"<<endl;
    cout<<" points B\t"<<m[b][0]<<"\t\t"<<m[b][1]<<"\t\t"<<m[b][2]<<"\t\t"<<endl;
    cout<<" points C\t"<<m[c][0]<<"\t\t"<<m[c][1]<<"\t\t"<<m[c][2]<<"\t\t"<<endl;       
            
 
    system("pause");
        //Освобождение памяти:
    for (int i = 0;i<dimm;i++) delete m[i];
    delete m;            
    return 0;
}