Определить номера точек, которые могут являться вершинами квадрата

#include<iostream>
#include<iomanip>
#include<cmath>
using namespace std;
double dist(double *a, int x, int y){
       return sqrt((a[x]-a[y])*(a[x]-a[y])+(a[x+1]-a[y+1])*(a[x+1]-a[y+1]));
}
bool isosceles(double *a, int x, int y, int z){
     if (dist(a, x, y)==dist(a, x, z) && dist(a, x, y)==dist(a, y, z) && a[x]!=a[y]) return true;
     return false;
}
int main(){
    int n;
    cout<<"n=";
    cin>>n;
    double *a=new double[2*n];
    for (int i=0; i<2*n; i+=2){
        cout<<"\nx"<<i/2+1<<" y"<<i/2+1<<" :\n";
        cin>>a[i]>>a[i+1];
    }
    cout<<"\nisosceles triangles :\n";
    for (int x=0; x<2*n; x+=2)
        for (int y=x+2; y<2*n; y+=2)
            for (int z=y+2; z<2*n; z+=2)
                if (isosceles(a, x, y, z))
                   cout<<setprecision(15)<<"x1="<<a[x]<<"\ty1="<<a[x+1]<<"\tx2="<<a[y]<<"\ty2="<<a[y+1]<<"\tx3="<<a[z]<<"\ty3="<<a[z+1]<<endl;
    system("pause");
    return 0;
}

/*
В одномерном массиве с чётным количеством элементов (2N) находятся 
координаты N точек плоскости. Они располагаются 
в следующем порядке: x1, y1, x2, y2, x3, y3, и т.д. 
Определить номера точек, которые могут являться вершинами квадрата.
*/
#include <algorithm>
#include <complex>
#include <iostream>
#include <iterator>
#include <vector>
///////////////////////////////////////////////////////////////////////////////
int const   SQUARE_VERTICES_TOTAL   =   4;
///////////////////////////////////////////////////////////////////////////////
typedef int                             T_coord;
typedef std::complex    < T_coord   >   T_point;  
typedef std::vector     < T_point   >   T_points;
typedef T_points                        T_square;
typedef std::vector     < T_square  >   T_squares;
typedef std::vector     < T_coord   >   T_norms;
///////////////////////////////////////////////////////////////////////////////
T_coord     get_norm
    (
        T_point  const   &   a,     
        T_point  const   &   b
    )
{
    return  std::norm( a - b );
}
///////////////////////////////////////////////////////////////////////////////
bool    is_square( T_square     const   &   square )
{
    T_norms     norms;
 
    for( size_t  i = 0; i < square.size() - 1; ++i )
    {
        for( size_t  j = i + 1; j < square.size(); ++j )
        {
            norms.push_back
                (
                    get_norm
                        (
                            square[i],
                            square[j]
                        )   
                );
        }//for
    }//for
 
    std::sort
        (
            norms.begin     (),
            norms.end       ()
        );
 
    return      std::distance
                    (
                        norms.begin(),
 
                        std::unique
                            (
                                norms.begin     (),
                                norms.end       ()
                            )
                    )
 
            ==  2;
}
///////////////////////////////////////////////////////////////////////////////
void    set_squares
    (
        T_points    const   &   points,
        T_squares           &   squares
    )
{
    for( size_t  i_0 = 0; i_0 < points.size() - 3; ++i_0 )
    {
        for( size_t  i_1 = i_0 + 1; i_1 < points.size() - 2; ++i_1 )
        {
            for( size_t  i_2 = i_1 + 1; i_2 < points.size() - 1; ++i_2 )
            {
                for( size_t  i_3 = i_2 + 1; i_3 < points.size(); ++i_3 )
                {
                    T_square    square_cur;
 
                    square_cur.push_back( points[i_0] );
                    square_cur.push_back( points[i_1] );
                    square_cur.push_back( points[i_2] );
                    square_cur.push_back( points[i_3] );
 
                    if  (
                            is_square( square_cur )
                        )
                    {
                        squares.push_back( square_cur );
                    }//if
                }//for i_3
            }//for i_2
        }//for i_1
    }//for i_0
}
///////////////////////////////////////////////////////////////////////////////
template< typename  T_cont >
void    print_cont( T_cont  const   &   cont )
{
    std::copy
        (
            cont.begin                                      (),
            cont.end                                        (),
            std::ostream_iterator< T_cont::value_type >     ( std::cout, "\t" )
        );
 
    std::cout   <<  std::endl;
}
///////////////////////////////////////////////////////////////////////////////
int     main()
{
    std::locale::global(std::locale(""));
 
    int     n   =   0;
 
    do
    {
        std::cout   <<  "Введите количество точек >= "
                    <<  SQUARE_VERTICES_TOTAL
                    <<  ": ";
 
        std::cin    >>  n;
    }
    while( n < SQUARE_VERTICES_TOTAL );
 
    std::cout   <<  "\nВведите "
                <<  n
                <<  " различных точек с целыми координатами в виде (2,3): "
                <<  std::endl;
 
    T_points    points;
    T_point     point_cur;
 
    while   (
                    int ( 
                            points.size() 
                        )
 
                <   n
            )
    {
        do
        {
            std::cout   <<  "#"
                        <<  points.size() + 1
                        <<  ": ";
 
            std::cin    >>  point_cur;
        }
        while   (
                        std::find
                            (
                                points.begin    (),
                                points.end      (),
                                point_cur
                            )
 
                    !=  points.end()
                );
 
        points.push_back( point_cur );
    }//while
 
    T_squares   squares;
 
    set_squares
        (
            points,
            squares
        );
 
    if  (
            squares.empty()
        )
    {
        std::cout   <<  "Нет вершин квадрата."
                    <<  std::endl;
    }
    else
    {
        int     ind_cur     =   0;
 
        for     ( 
                    auto    const   &   square_cur 
                    :
                    squares
                )
        {
            std::cout   <<  "\nКвадрат #"
                        <<  ++ind_cur
                        <<  ":\n";
 
            print_cont( square_cur );
        }
    }//else
 
    system("pause");
}