Использование массивов для решения геометрических задач

#include <iostream>
using namespace std;
int main()
{
            setlocale(0, "");
            int i,j,k,n,m;
            double a,b,c;
            n=28;
            n++;
        double *x = new double[n];
        double *y = new double[n];
        int ***check = new int**[n];
    for(i=0; i<n; i++){
    check[i] = new int*[n]; 
    for(j=0; j<n;j++){
        check[i][j] = new int[n]; 
    }
}
            n=n-1;
           // Координаты точек
            for(i=1;i<n;i++)
            {
                x[i]=0.001*rand();
                y[i]=0.001*rand();
                //cout<<x[i]<<" "<<y[i]<<endl;
            }
            //три точки для проверки
            x[7]=0;    y[7]=0;
            x[14]=0.5; y[14]=1;
            x[21]=1;   y[21]=0;
            //Массив контроля 
            for(i=1;i<n;i++)
            for(j=1;j<n;j++)
            for(k=1;k<n;k++) check[i][j][k]=0;
 
            for(i=1;i<n;i++)
            for(j=1;j<n;j++)
            for(k=1;k<n;k++)
             if(!(i==j) && !(k==j) && !(k==i) && check[i][j][k]==0){
                //cout<<i<<" "<<j<<" "<<k<<" "<<endl;
                check[i][j][k]=1;
                check[i][k][j]=1;
                check[j][i][k]=1;
                check[j][k][i]=1;
                check[k][i][j]=1;
                check[k][j][i]=1;
                //стороны треугольника
                a=sqrt(pow((x[i]-x[j]),2.0)+pow((y[i]-y[j]),2.0));
                b=sqrt(pow((x[i]-x[k]),2.0)+pow((y[i]-y[k]),2.0));
                c=sqrt(pow((x[k]-x[j]),2.0)+pow((y[k]-y[j]),2.0));
                if(a==b || a==c || b==c)
                cout<<a<<" "<<b<<" "<<c<<" "<<i<<" "<<j<<" "<<k<<endl;
             }
        delete [] x;
        delete [] y;
for (i = 0; i < n; i++)
        {
 
        for (j = 0; j < n; j++)
            delete check[i][j];
    }
    
    delete check;
            system("pause"); // Только для тех, у кого MS Visual Studio
}