Программа по нахождению всевозможных структур из четырех точек

#include "stdafx.h"
#include <iostream>
using namespace std;
//структура "Точка"
struct Point
{
    double x,y;
};
//структура "Квадрат"
struct Square
{
    Point p1,p2,p3,p4;
};
int fact(int N)
{
    if(N < 0) // если пользователь ввел отрицательное число
        return 0; // возвращаем ноль
    if (N == 0) // если пользователь ввел ноль,
        return 1; // возвращаем факториал от нуля
    else // Во всех остальных случаях
        return N * fact(N - 1); // делаем рекурсию.
}
int _tmain(int argc, _TCHAR* argv[])
{
    //ввод множества точек
    cout<<"Please, enter N:\n";
    int np;
    cin>>np;
    Point *pp=new Point[np];
    cout<<"Please, enter N points:\n";
    for(int i=0;i<np;i++)
        cin>>pp[i].x>>pp[i].y;
    int b =0;
    for(int i=0;i<np;i++)
        cout<<pp[i].x<<"   "<<pp[i].y<<endl;
    b = ((fact(np))/(fact(4)*(fact(np-4))));
    cout<<b<<endl;
    Square *squares = new Square[np];
    
    for (int i =0;i<np;i++){
        for( int j=1;j<i-1;j++){
            for ( int k=2;k<j-2;k++){
                for ( int l = 3;l<k-3;l++){
                    squares[i].p1=pp[i];
                    cout<<squares[1].p1.x<<" " <<squares[i].p1.y<<endl;
                    squares[i].p2=pp[j];
                    squares[i].p3=pp[k];
                    squares[i].p4=pp[l];
                }}}}
    cout<<"test"<<endl;
    for (int i =0;i<np;i++){
        cout<<squares[i].p1.x<<"  "<<squares[i].p1.y<<"     "<<squares[i].p2.x<<" "<<squares[i].p2.y<<"      "<<squares[i].p3.x<<"  "<<squares[i].p3.y<<"     "<<squares[i].p4.x<<" "<<squares[i].p4.y<<endl;}
    return 0;
}