Не могу реализовать

#include <iostream>
#include <stdio.h>
#include <math.h>
using namespace std;
typedef struct l* pl;
struct l 
{
pl pe;
pl ne;
double x;
int t;
l(){pe=NULL;ne=NULL;};
};
class plf
{
public:
pl cr(){ return new l;};
};
plf fa;
pl ty;//указатель на первый элемент списка
pl ct; //указатель на текущий элемент списка
pl g;
double **a, **pairs3, **tr;
int st,tri;
int fact(int y);
int C(int m,int n);
double** init_d_data(double** k, int t, int y);
void delete_d_data(double** k, int t, int y);
void print_d_data(double** k, int t, int y);
void points3(double** m, double** n, int m1);
void check_data(double** m,  int m1, int m2);
int intersection(double ax1,double ay1,double ax2,double ay2,double bx1,double by1,double bx2,double by2);
double** add_d_data(double** k, int t, int y, int o);
int main(int argc, char* argv[]) 
{ 
int Cnm,i;
int a2,sp21,sp22,sp31,sp32;
st = 0;
// количество точек
cin >> a2;
a = init_d_data(a,a2,2);
// задание точек
for (i=0;i<a2;i++){ cin >> a[i][0]; cin >> a[i][1];}
// количество способов выбрать 3 элемента из a2
Cnm = C(3,a2); printf("C(%d,%d) =  %d\n", 3, a2, Cnm);
sp31 = Cnm; sp32 = 6;
// найти тройки точек (треугольники)
pairs3=init_d_data(pairs3,sp31,sp32); 
points3(a,pairs3,a2);
// напечатать тройки точек (треугольники)
printf("(pairs3):\n");
print_d_data(pairs3,sp31,sp32);
//  проверить на равносторонесть
check_data(pairs3,sp31,sp32);
//  уничтожить данные
delete_d_data(a,a2,2);
delete_d_data(pairs3,sp31,sp32);//
return 0;
}
// факториал
int fact(int y){if (y>1){ return y*fact(y-1);} else {return 1;}}
// количество сочетаний
int C(int m,int n){ return ( fact(n)/(fact(n-m)* fact(m)));}
// создание динамических массивов
double** init_d_data(double** k, int t, int y){int i; k = new double*[t];for (i=0;i<t;i++){k[i]=new double[y];} return k;}
// удаление динамических массивов
void delete_d_data(double** k, int t, int y){ int i; for (i=0;i<t;i++){delete[] k[i];} delete[] k;}
// печать динамических массивов
void print_d_data(double** k, int t, int y){ int i,j; for (i=0;i<t;i++){for (j=0;j<y;j++){printf("%.2f ",k[i][j]);}printf("\n");}}
// добавление o элементов в динамический массив
double** add_d_data(double** k, int t, int y, int o)
{
double** tmp =init_d_data(tmp,t+o,y);
int i,j;
for (i=0;i<t;i++){
for (j=0;j<y;j++){
tmp[i][j]=k[i][j];}}
delete_d_data(k,t,y);
return(tmp);
}
// выбор 3 точек
void points3(double** m, double** n, int m1)
{
int i,j,k,l;k=0; 
for (i=0;i<m1-2;i++){
for (j=i+1;j<m1-1;j++){
for (l=j+1;l<m1;l++){
n[k][0]=m[i][0]; 
n[k][1]=m[i][1];
n[k][2]=m[j][0]; 
n[k][3]=m[j][1];
n[k][4]=m[l][0]; 
n[k][5]=m[l][1];k++;}}}
}
void check_data(double** m,  int m1, int m2)
{
int i,k;
double eps =0.1;
double sc1,sc2,sc3;
for (i=0;i<m1;i++) { k=0;
//найти расстояния между парами точек - если они все равны треугольник равносторонний
sc1= sqrt((m[i][2] - m[i][0])*(m[i][2] - m[i][0]) + (m[i][3] - m[i][1])*(m[i][3] - m[i][1]));
sc2= sqrt((m[i][4] - m[i][0])*(m[i][4] - m[i][0]) + (m[i][5] - m[i][1])*(m[i][5] - m[i][1]));
sc3= sqrt((m[i][4] - m[i][2])*(m[i][4] - m[i][2]) + (m[i][5] - m[i][3])*(m[i][5] - m[i][3]));
if (fabs(sc1-sc2)<eps){k++;} if (fabs(sc1-sc3)<eps){k++;} if (fabs(sc2-sc3)<eps){k++;}
printf("%d %f %d %f %d %f  %f %f %f %d\n",(fabs(sc1-sc2)<eps),fabs(sc1-sc2),(fabs(sc1-sc3)<eps),fabs(sc1-sc3),(fabs(sc2-sc3)<eps),fabs(sc2-sc3), sc1,sc2,sc3,k);
if (k==3) { //треугольник равносторонний
if (st==0) {
// список размеров треугольников и количества с определенным размером
ty = fa.cr();ty->x=sc1; ty->t=1;ct=ty; st=1;
// массив вершин равносторонних треугольников
tr = init_d_data(tr,1,6);tri=1;
tr[0][0] =m[i][0];tr[0][1] =m[i][1];
tr[0][2] =m[i][2];tr[0][3] =m[i][3];
tr[0][4] =m[i][4];tr[0][5] =m[i][5];
printf("tr: %.2f %.2f %.2f %.2f %.2f %.2f\n", tr[0][0],tr[0][1],tr[0][2],tr[0][3],tr[0][4],tr[0][5]);
}
else {
//проходим список с начала - если такой треугольник есть до увеличиваем счетчик -- если нет добавляем в конец
g = ty; 
while (g->ne!=NULL) 
{ if (fabs(g->x - sc1)<eps){ g->t+=1; }; g=g->ne;} 
if (fabs(g->x - sc1)<eps){ g->t+=1; }
else {pl nt = fa.cr(); nt->x = sc1; nt->t=1; nt->pe = ct; (nt->pe)->ne = nt;ct = nt; }
tr=add_d_data(tr,tri,6,1);
tr[tri][0] =m[i][0];tr[tri][1] =m[i][1];
tr[tri][2] =m[i][2];tr[tri][3] =m[i][3];
tr[tri][4] =m[i][4];tr[tri][5] =m[i][5];
printf("tr: %.2f %.2f %.2f %.2f %.2f %.2f\n", tr[tri][0],tr[tri][1],tr[tri][2],tr[tri][3],tr[tri][4],tr[tri][5]);
tri++;printf("%d\n",tri);
}}}
// tr - массив равносторонних треугольников tri - размер
int j;
if (tri>1){
for (i=0;i<tri-1;i++) {
for (j=i+1;j<tri;j++) {
int e,h,p,w;
for (e=0; e<6;e+=2){
for (h=e+2; h<6;h+=2){
for (p=0; p<6;p+=2){
for (w=p+2; w<6;w+=2){
printf ("indexes: e %d  h %d p %d  w %d\n", e,h,p,w);
printf ("treugolnik 1: %.2f %.2f %.2f %.2f\n", tr[i][e],tr[i][e+1],tr[i][h],tr[i][h+1]);
printf ("treugolnik 2: %.2f %.2f %.2f %.2f\n", tr[j][p],tr[j][p+1],tr[j][w],tr[j][w+1]);
//проверка пересечения треугольников qw = 1 если треугольники пересекаются
int qw=intersection(tr[i][e],tr[i][e+1],tr[i][h],tr[i][h+1],tr[j][p],tr[j][p+1],tr[j][w],tr[j][w+1]);
printf("intersection : %d\n",qw);
}}}}
 
}}}
 
}
//проверка пересечения отрезков: [URL]http://acmp.ru/article.asp?id_text=170[/URL]
int intersection(double ax1,double ay1,double ax2,double ay2,double bx1,double by1,double bx2,double by2){
return((((bx2-bx1)*(ay1-by1)-(by2-by1)*(ax1-bx1))*((bx2-bx1)*(ay2-by1)-(by2-by1)*(ax2-bx1)))<=0  &&
   (((ax2-ax1)*(by1-ay1)-(ay2-ay1)*(bx1-ax1))*((ax2-ax1)*(by2-ay1)-(ay2-ay1)*(bx2-ax1)))<=0);}