Из заданного на плоскости множества точек выбрать три различные точки

// epiphanies.cpp: определяет точку входа для консольного приложения.
//
/* каждый пункт в виде шаблона функции
   данные в кач. параметров
   примеры программ, использующих шаблоны для int, float, double
   
   "из заданного на плоскости множества точек выбрать три различные точки так, 
   чтобы разность между площадью круга, ограниченного окружнотью, проходящей через эти три точки, 
   и площадью треугольника с вершинами в этих точках была миимальной"  */
 
/*1) создать массив x,y; создать массив для конечных результатов 1 итерации функции с площадями
2) создать функцию, выбирающую 3 точки и считающую площадь треугольника и площдь круга (полученные результаты записывать в табличку!)
3) создать другую функцию, которая будет выполнять сортировку по площадям, дабы найти минимальную.
+++++
x,y вбиваются пользователем!*/
 
 
#include "stdafx.h"
#include <iostream>
#include <iomanip>
#include <fstream>
#include <cstring>
 
#include <math.h>
//#define _USE_MATH_DEFINE
 
 
 
 
using namespace std;
 
template <class type> void create_massiv (int lstr, i, j, type **a);
template <class type2> void points (int lstr, i, j, type2 **a);
template <class type3> void sorting (int i, j, str, type3 **s);
 
 
 
template <class type> 
void create_massiv (int lstr, int i, int j, type **a) 
{
    
    cout << "Enter matrix: ";
    
    for (i = 0; i < lstr; i++) {
        for (j = 0; j < 2; j++) {
            cin >> a[i][j];
        }
        cout << endl;
    }
    
    cout << "  X  " << " | " << "  Y  " << endl;
    cout << endl;
    
    for (i = 0; i < lstr; i++) {
        for (j = 0; j < 2; j++){ 
            cout << setw(4) << a[i][j] << "  ";
            cout << endl;
        }
    }
    
    
}
 
 
 
template <class type2> 
void points (int lstr, i, j, str, type2 **a, type2 **s) 
{
    type2 const pi = 3.14;
    type2 x1, x2, x3, y1, y2, y3, ma, mb, x, y, r, SC;
    
    for (i = 0; i < lstr - 2; i++){
        for (j = i + 1; j < lstr - 1; j++) {
            for (int k = j + 1; k < lstr; k++) {
                x1 = a[i][0];
                y1 = a[i][1];
 
                x2 = a[i+1][0];
                y2 = a[i+1][0];
  
                x3 = a[i+3][0];
                y3 = a[i+3][0];
 
                cout << " P1(x1;y1) =  ( " << x1 << " ; " << y1 << " )" << endl;
                cout << " P2(x2;y2) =  ( " << x2 << " ; " << y2 << " )" << endl;
                cout << " P3(x3;y3) =  ( " << x3 << " ; " << y3 << " )" << endl;
                
                ma = (y1-y2)/(x2-x1);
                mb = (y3-y2)/(x3-x2);
 
                if (mb - ma == 0) {
                    cout << "P1P2 || P2P3 : these points don't form a circle " << endl;
                    break;
                }
 
                x = (ma*mb*(y1-y3)+mb*(x1+x2)-ma*(x2+x3))/2*(mb-ma);
                y = (-1)*(x - (x1+x2)/2)/ma +(y1+y2)/2;
                //cout << "The center O(x;y)" << x << "  " << y << endl;
                
                r = sqrt( (x1-x)*(x1-x) - (y1-y)*(y1-y) );
                //cout << "The radius " << r << endl;
                
                SC = pi*r*r;
                //cout << "Square of the circle " << SC << endl;
                
                type2 ST = ((x1-x3)*(y2-y3) - (x2-x3)*(y1-y3))/2;
                
                cout << "  №  " << "|" << " S of circle " << "|" << " S of tringle " << "|" << "   SC - ST  " << endl;
                for (int f = 0; f < str; f++){
                    s[f][0] = f;
                    s[f][1] = SC;
                    s[f][3] = ST;
                    s[f][4] = SC - ST;
                }
                
                for (int f = 0; f < str; f++) {
                    for (int g = 0; g < 4: g++) {
                        cout << setw(4) << s[f][g] << "  ";
                        cout << endl;
                    }
                }
            }
        } 
    }
}
 
 
 
 
 
template <class type3> 
void sorting (int i, j, str, type3 **s)
{
    for (j = str-1; j > 0; j--){
        for (i = 0; i < j; i++) {
            if (s[i+1][4] < s[i][4]) {
                type3 min = s[i];
                s[i] = s[i+1];
                s[i+1] = min;
            } 
        }
    }
    
    cout << "sorted data" << endl;
    cout << endl;
    cout << "  №  " << "|" << " S of circle " << "|" << " S of tringle " << "|" << "   SC - ST  " << endl;
    
    for (i = 0; i < str; i++) {
        for (j = 0; j < 4: j++) {
            cout << setw(4) << s[i][j] << "  ";
            cout << endl;
        }
    }
}
 
 
 
 
int main()
{
    int lstr;
    int i, j;
    cout << "Enter lstr strings" << endl;
    cin >> lstr;
    int str = lstr;
    
    double **a = new double*[lstr];
    for (i = 0; i < lstr; i++){ 
        a[i] = new double[2];
    }
    
    double **s = new double*[str];
    for (i = 0; i < str; i++){ 
        s[i] = new double[4];
    }
 
    cout << "With the type DOUBLE" << endl;
    create_massiv <double> (lstr, i, j, **a);
    points <double> (/*lstr, i, j,str, **a,*/ **s);
    sorting <double> (/*i, j, str,*/ **s);
    
    /*float **a = new float*[lstr];
    for (i = 0; i < lstr; i++){ 
        a[i] = new float[2];
    }
    
    float **s = new float*[str];
    for (i = 0; i < str; i++){ 
        s[i] = new float[4];
    }*/
 
    cout << "With the type FLOAT" << endl;
    create_massiv <float> (lstr, i, j, **a);
    points <float> (/*lstr, i, j,str, **a,*/ **s);
    sorting <float> (/*i, j, str, */**s);
    
    /*int **a = new int*[lstr];
    for (i = 0; i < lstr; i++){ 
        a[i] = new int[2];
    }
    
    int **s = new int*[str];
    for (i = 0; i < str; i++){ 
        s[i] = new int[4];
    }*/
 
    cout << "With the type INT" << endl;
    create_massiv <int> (lstr, i, j, **a);
    points <int> (/*lstr, i, j,str, **a,*/ **s);
    sorting <int> (/*i, j, str,*/ **s);
    
    return 0;
}

// happiness.cpp: определяет точку входа для консольного приложения.
//
/* каждый пункт в виде шаблона функции
   данные в кач. параметров
   примеры программ, использующих шаблоны для int, float, double
   
   "из заданного на плоскости множества точек выбрать три различные точки так, 
   чтобы разность между площадью круга, ограниченного окружнотью, проходящей через эти три точки, 
   и площадью треугольника с вершинами в этих точках была миимальной"  */
 
/*1) создать массив x,y; создать массив для конечных результатов 1 итерации функции с площадями
2) создать функцию, выбирающую 3 точки и считающую площадь треугольника и площдь круга (полученные результаты записывать в табличку!)
3) создать другую функцию, которая будет выполнять сортировку по площадям, дабы найти минимальную.
+++++
x,y вбиваются пользователем!*/
 
 
#include "stdafx.h"
#include <iostream>
#include <iomanip>
#include <fstream>
#include <cstring>
#include <stdlib.h>
#include <math.h>
 
using namespace std;
 
 
template <class type> 
void func (int lstr, type **a1, type **s1)
{
    int i, j;
type **s = new type*[lstr];
    for ( i = 0; i < lstr; i++){ 
        s[i] = new type[4];
    }
 
type **a = new type*[lstr];
    for ( i = 0; i < lstr; i++){ 
        a[i] = new type[2];
    }
    
    
    
    cout << "Enter matrix: ";
    
    for ( i = 0; i < lstr; i++) {
        for (j = 0; j < 2; j++) {
            cin >> a[i][j];
        }
        cout << endl;
    }
    
    cout << "  X  " << " | " << "  Y  " << endl;
    cout << endl;
    
    for (i = 0; i < lstr; i++) {
        for (j = 0; j < 2; j++){ 
            cout << setw(4) << a[i][j] << "  ";
            cout << endl;
        }
    }
    
double const pi = 3.14;
    type x1, x2, x3, y1, y2, y3, ma, mb, x, y;
    
    for (i = 0; i < lstr - 2; i++){
        for (j = i + 1; j < lstr - 1; j++) {
            for (int k = j + 1; k < lstr; k++) {
                x1 = a[i][0];
                y1 = a[i][1];
 
                x2 = a[i+1][0];
                y2 = a[i+1][0];
  
                x3 = a[i+3][0];
                y3 = a[i+3][0];
 
                cout << " P1(x1;y1) =  ( " << x1 << " ; " << y1 << " )" << endl;
                cout << " P2(x2;y2) =  ( " << x2 << " ; " << y2 << " )" << endl;
                cout << " P3(x3;y3) =  ( " << x3 << " ; " << y3 << " )" << endl;
                
                ma = (y1-y2)/(x2-x1);
                mb = (y3-y2)/(x3-x2);
 
                if (mb - ma == 0) {
                    cout << "P1P2 || P2P3 : these points don't form a circle " << endl;
                    break;
                }
 
                x = (ma*mb*(y1-y3)+mb*(x1+x2)-ma*(x2+x3))/2*(mb-ma);
                y = (-1)*(x - (x1+x2)/2)/ma +(y1+y2)/2;
                //cout << "The center O(x;y)" << x << "  " << y << endl;
                
                double /*type*/ r = sqrt( (x1-x)*(x1-x) - (y1-y)*(y1-y) );
                //cout << "The radius " << r << endl;
                
                double /*type*/ SC = pi*r*r;
                //cout << "Square of the circle " << SC << endl;
                
                type ST = ((x1-x3)*(y2-y3) - (x2-x3)*(y1-y3))/2;
                
                cout << "  №  " << "|" << " S of circle " << "|" << " S of tringle " << "|" << "   SC - ST  " << endl;
                for (int f = 0; f < lstr; f++){
                    s[f][0] = f;
                    s[f][1] = SC;
                    s[f][3] = ST;
                    s[f][4] = SC - ST;
                }
                
                for (int f = 0; f < lstr; f++) {
                    for (int g = 0; g < 4; g++) {
                        cout << setw(4) << s[f][g] << "  ";
                        cout << endl;
                    }
                }
            }
        } 
    }
 
for (j = lstr-1; j > 0; j--){
        for (i = 0; i < j; i++) {
            if (s[i+1][4] < s[i][4]) {
                type min = s[i];
                s[i] = s[i+1];
                s[i+1] = min;
            } 
        }
    }
    
    cout << "sorted data" << endl;
    cout << endl;
    cout << "  №  " << "|" << " S of circle " << "|" << " S of tringle " << "|" << "   SC - ST  " << endl;
    
    for (i = 0; i < lstr; i++) {
        for (j = 0; j < 4; j++) {
            cout << setw(4) << s[i][j] << "  ";
            cout << endl;
        }
    }
 
 
} 
 
 
int main()
{
    int lstr;
    cout << "Enter lstr strings" << endl;
    cin >> lstr;
    int x;
 
    while ( x < 4 ){
        cout << "Choose the type of array" << endl;
        cout << "1 - int" << endl;
        cout << "2 - double" << endl;
        cout << "3 - float" << endl;
        cin >> x;
        
        switch(x) 
        {
        case 1: int ** s1;
                int ** a1;
                func ( lstr, a1, s1);
                break;
        
        case 2: double **s2;
                double **a2;
                func ( lstr, a2, s2);
                break;
        
        case 3: float **s3;
                float **a3;
                func ( lstr, a3, s3);
                break;
        }
    }
    
    return 0;
}

// happiness.cpp: определяет точку входа для консольного приложения.
//
/* каждый пункт в виде шаблона функции
   данные в кач. параметров
   примеры программ, использующих шаблоны для int, float, double
   
   "из заданного на плоскости множества точек выбрать три различные точки так, 
   чтобы разность между площадью круга, ограниченного окружнотью, проходящей через эти три точки, 
   и площадью треугольника с вершинами в этих точках была миимальной"  */
 
/*1) создать массив x,y; создать массив для конечных результатов 1 итерации функции с площадями
2) создать функцию, выбирающую 3 точки и считающую площадь треугольника и площдь круга (полученные результаты записывать в табличку!)
3) создать другую функцию, которая будет выполнять сортировку по площадям, дабы найти минимальную.
+++++
x,y вбиваются пользователем!*/
 
 
//#include "stdafx.h"
#include <iostream>
#include <iomanip>
#include <fstream>
#include <cstring>
#include <cstdlib>
#include <cmath>
 
using namespace std;
 
 
template <class type> 
void func (int lstr, type **a1, type **s1)
{
    int i, j;
    type **s = new type*[lstr];
    for ( i = 0; i < lstr; i++){ 
        s[i] = new type[4];
    }
 
    type **a = new type*[lstr];
    for ( i = 0; i < lstr; i++){ 
        a[i] = new type[2];
    }
    
    
    
    cout << "Enter matrix: ";
    
    for ( i = 0; i < lstr; i++) {
        for (j = 0; j < 2; j++) {
            cin >> a[i][j];
        }
        cout << endl;
    }
    
    cout << "  X  " << " | " << "  Y  " << endl;
    cout << endl;
    
    for (i = 0; i < lstr; i++) {
        for (j = 0; j < 2; j++){ 
            cout << setw(4) << a[i][j] << "  ";
        }
        cout << endl;
    }
    
    double const pi = M_PI;
    type x1, x2, x3, y1, y2, y3, ma, mb, x, y;
    
    for (i = 0; i < lstr - 2; i++){
        for (j = i + 1; j < lstr - 1; j++) {
            for (int k = j + 1; k < lstr; k++) {
                x1 = a[i][0];
                y1 = a[i][1];
 
                x2 = a[i+1][0];
                y2 = a[i+1][0];
  
                x3 = a[i+3][0];
                y3 = a[i+3][0];
 
                cout << " P1(x1;y1) =  ( " << x1 << " ; " << y1 << " )" << endl;
                cout << " P2(x2;y2) =  ( " << x2 << " ; " << y2 << " )" << endl;
                cout << " P3(x3;y3) =  ( " << x3 << " ; " << y3 << " )" << endl;
                
                ma = (y1-y2)/(x2-x1);
                mb = (y3-y2)/(x3-x2);
 
                if (mb - ma == 0) {
                    cout << "P1P2 || P2P3 : these points don't form a circle " << endl;
                    break;
                }
 
                x = (ma*mb*(y1-y3)+mb*(x1+x2)-ma*(x2+x3))/2*(mb-ma);
                y = (-1)*(x - (x1+x2)/2)/ma +(y1+y2)/2;
                //cout << "The center O(x;y)" << x << "  " << y << endl;
                
                type r =(type) sqrt(double( (x1-x)*(x1-x) - (y1-y)*(y1-y) ));
                //cout << "The radius " << r << endl;
                
                double /*type*/ SC = pi*r*r;
                //cout << "Square of the circle " << SC << endl;
                
                type ST = ((x1-x3)*(y2-y3) - (x2-x3)*(y1-y3))/2.0;
                
                cout << "  №  " << "|" << " S of circle " << "|" << " S of tringle " << "|" << "   SC - ST  " << endl;
                for (int f = 0; f < lstr; f++){
                    s[f][0] = f;
                    s[f][1] = SC;
                    s[f][3] = ST;
                    s[f][4] = SC - ST;
                }
                
                for (int f = 0; f < lstr; f++) {
                    for (int g = 0; g < 4; g++) {
                        cout << setw(4) << s[f][g] << "  ";
                        cout << endl;
                    }
                }
            }
        } 
    }
 
for (j = lstr-1; j > 0; j--){
        for (i = 0; i < j; i++) {
            if (s[i+1][4] < s[i][4]) {
                type* min = s[i];
                s[i] = s[i+1];
                s[i+1] = min;
            } 
        }
    }
    
    cout << "sorted data" << endl;
    cout << endl;
    cout << "  №  " << "|" << " S of circle " << "|" << " S of tringle " << "|" << "   SC - ST  " << endl;
    
    for (i = 0; i < lstr; i++) {
        for (j = 0; j < 4; j++) {
            cout << setw(4) << s[i][j] << "  ";
        }
        cout << endl;
    }
 
 
} 
 
 
int main()
{
    int lstr;
    cout << "Enter lstr strings" << endl;
    cin >> lstr;
    int x;
 
    while ( x < 4 ){
        cout << "Choose the type of array" << endl;
        cout << "1 - int" << endl;
        cout << "2 - double" << endl;
        cout << "3 - float" << endl;
        cin >> x;
        
        switch(x) 
        {
        case 1: int ** s1;
                int ** a1;
                func ( lstr, a1, s1);
                break;
        
        case 2: double **s2;
                double **a2;
                func ( lstr, a2, s2);
                break;
        
        case 3: float **s3;
                float **a3;
                func ( lstr, a3, s3);
                break;
        }
    }
    
    return 0;
}

// happiness.cpp: определяет точку входа для консольного приложения.
//
/* каждый пункт в виде шаблона функции
   данные в кач. параметров
   примеры программ, использующих шаблоны для int, float, double
   
   "из заданного на плоскости множества точек выбрать три различные точки так, 
   чтобы разность между площадью круга, ограниченного окружнотью, проходящей через эти три точки, 
   и площадью треугольника с вершинами в этих точках была миимальной"  */
 
/*1) создать массив x,y; создать массив для конечных результатов 1 итерации функции с площадями
2) создать функцию, выбирающую 3 точки и считающую площадь треугольника и площдь круга (полученные результаты записывать в табличку!)
3) создать другую функцию, которая будет выполнять сортировку по площадям, дабы найти минимальную.
+++++
x,y вбиваются пользователем!*/
 
 
#include "stdafx.h"
#include <iostream>
#include <iomanip>
#include <fstream>
#include <cstring>
#include <cstdlib>
#include <cmath>
 
using namespace std;
 
 
template <class type> 
void func (int lstr, type **a1, type **s1)
{
    int i, j;
    double const pi = 3.14;
    type x1, x2, x3, y1, y2, y3, x, y;
    double ma, mb;
    
    type **s = new type*[lstr];
    for ( i = 0; i < lstr; i++){
        s[i] = new type[4];
    }
   
    type **a = new type*[lstr];
    for ( i = 0; i < lstr; i++){ 
        a[i] = new type[2];
    }
    
    
    
    cout << "Enter matrix: ";
    
    for ( i = 0; i < lstr; i++) {
        for (j = 0; j < 2; j++) {
            cin >> a[i][j];
        }
        cout << endl;
    }
    
    cout << "  X  " << " | " << "  Y  " << endl;
    cout << endl;
    
    for (i = 0; i < lstr; i++) {
        for (j = 0; j < 2; j++){ 
            cout << setw(4) << a[i][j] << "  ";
        }
        cout << endl;
    }
    
    
    
    for (i = 0; i < lstr - 2; i++){
        for (j = i + 1; j < lstr - 1; j++) {
            for (int k = j + 1; k < lstr; k++) {
                x1 = a[i][0];
                y1 = a[i][1];
 
                x2 = a[j][0];
                y2 = a[j][1];
  
                x3 = a[k][0];
                y3 = a[k][1];
 
                cout << " P1(x1;y1) =  ( " << x1 << " ; " << y1 << " )" << endl;
                cout << " P2(x2;y2) =  ( " << x2 << " ; " << y2 << " )" << endl;
                cout << " P3(x3;y3) =  ( " << x3 << " ; " << y3 << " )" << endl;
                
                ma = (y2-y1)/(x2-x1);
                cout << "ma" <<  endl; cout <<ma <<endl;
                mb = (y3-y2)/(x3-x2);
                cout << "mb" << mb << endl;
 
                if (mb - ma == 0) {
                    cout << "P1P2 || P2P3 : these points don't form a circle " << endl;
                    break;
                }
 
                 x =  (ma*mb*(y1-y3)+mb*(x1+x2)-ma*(x2+x3))/2*(mb-ma);
                 y =  ((-1)*(x - (x1+x2)/2))/ma +(y1+y2)/2;
                
                cout << "The center O(x;y)" << x << "  " << y << endl; // начинаются проблемы!!!!!!!!
 
                
                type r =(type) sqrt(double( (x1-x)*(x1-x) - (y1-y)*(y1-y) )); 
                cout << "The radius " << r << endl;
                
                double SC = pi*r*r;
                cout << "Square of the circle " << SC << endl;
                
                double ST = ((x1-x3)*(y2-y3) - (x2-x3)*(y1-y3))/2.0;
                
                cout << "  №  " << "|" << " S of circle " << "|" << " S of tringle " << "|" << "   SC - ST  " << endl;
                
                
                for (int f = 0; f < lstr; f++){
                    s[f][0] = f;
                    s[f][1] = SC;
                    s[f][3] = ST;
                    s[f][4] = SC - ST;
                }
                
                for (int f = 0; f < lstr; f++) {
                    for (int g = 0; g < 4; g++) {
                        cout << setw(4) << s[f][g] << "  ";
                        cout << endl;
                    }
                }
                
                for (j = lstr-1; j > 0; j--){
                    for (i = 0; i < j; i++) {
                        if (s[i+1][4] < s[i][4]) {
                            type* min = s[i];
                            s[i] = s[i+1];
                            s[i+1] = min;
                        } 
                    }
                }
                cout << "sorted data" << endl;
                cout << endl;
                cout << "  №  " << "|" << " S of circle " << "|" << " S of tringle " << "|" << "   SC - ST  " << endl;
                
                for (i = 0; i < lstr; i++) {
                    for (j = 0; j < 4; j++) {
                        cout << setw(4) << s[i][j] << "  ";
                    }
                    cout << endl;
                }
            }
        } 
    }
 
 
 
} 
 
 
int main()
{
    int lstr;
    cout << "Enter lstr strings" << endl;
    cin >> lstr;
    int q;
 
    do {     
        cout << "Choose the type of array" << endl;
        cout << "1 - int" << endl;
        cout << "2 - double" << endl;
        cout << "3 - float" << endl;
        cin >> q;
        
        switch(q) 
        {
        case 1: int ** s1;
                int ** a1;
                func ( lstr, a1, s1);
                break;
        
        case 2: double **s2;
                double **a2;
                func ( lstr, a2, s2);
                break;
        
        case 3: float **s3;
                float **a3;
                func ( lstr, a3, s3);
                break;
        }
    } while (q != 4);
    
    
    return 0;
}

#include "stdafx.h"
#include <iostream>
#include <iomanip>
#include <fstream>
#include <cstring>
#include <cstdlib>
#include <cmath>
 
using namespace std;
 
 
template <class type> 
void func (int lstr, type **a1, type **s1)
{
    int i, j;
    double const pi = 3.14;
    type x1, x2, x3, y1, y2, y3, x, y;
    double ma, mb;
    
    type **s = new type*[lstr];
    for ( i = 0; i < lstr; i++){
        s[i] = new type[4];
    }
   
    type **a = new type*[lstr];
    for ( i = 0; i < lstr; i++){ 
        a[i] = new type[2];
    }
    
    
    
    cout << "Enter matrix: ";
    
    for ( i = 0; i < lstr; i++) {
        for (j = 0; j < 2; j++) {
            cin >> a[i][j];
        }
        cout << endl;
    }
    
    cout << "  X  " << " | " << "  Y  " << endl;
    cout << endl;
    
    for (i = 0; i < lstr; i++) {
        for (j = 0; j < 2; j++){ 
            cout << setw(4) << a[i][j] << "  ";
        }
        cout << endl;
    }
    
    
    
    for (int ii = 0; ii < lstr - 2; ii++){
        for (int jj = ii + 1; jj < lstr - 1; jj++) {
            for (int kk = jj + 1; kk < lstr; kk++) {
                x1 = a[ii][0];
                y1 = a[ii][1];
 
                x2 = a[jj][0];
                y2 = a[jj][1];
  
                x3 = a[kk][0];
                y3 = a[kk][1];
 
                cout << " P1(x1;y1) =  ( " << x1 << " ; " << y1 << " )" << endl;
                cout << " P2(x2;y2) =  ( " << x2 << " ; " << y2 << " )" << endl;
                cout << " P3(x3;y3) =  ( " << x3 << " ; " << y3 << " )" << endl;
                
                /////////здесь будет ошибка в расчетах (бесконечность), если в знаменателе будет 0... а он там будет при равенстве x1 и x2 или x2 и x3
                ma = (y2-y1)/(x2-x1);
                cout << "ma" <<  endl; cout <<ma <<endl;
                mb = (y3-y2)/(x3-x2);
                cout << "mb" << mb << endl;
                //////////////
 
                if (mb - ma == 0) {
                    cout << "P1P2 || P2P3 : these points don't form a circle " << endl;
                    break;
                }
 
                 x =  (ma*mb*(y1-y3)+mb*(x1+x2)-ma*(x2+x3))/2*(mb-ma);
                 y =  ((-1)*(x - (x1+x2)/2))/ma +(y1+y2)/2;
                
                cout << "The center O(x;y)" << x << "  " << y << endl; // начинаются проблемы!!! проблемы, т.к. ma или mb - бесконечность..
 
                
                type r =(type) sqrt(double( (x1-x)*(x1-x) - (y1-y)*(y1-y) )); 
                cout << "The radius " << r << endl;
                
                double SC = pi*r*r;
                cout << "Square of the circle " << SC << endl;
                
                double ST = ((x1-x3)*(y2-y3) - (x2-x3)*(y1-y3))/2.0;
                
                cout << "  №  " << "|" << " S of circle " << "|" << " S of tringle " << "|" << "   SC - ST  " << endl;
                
                
                for (int f = 0; f < lstr; f++){
                    s[f][0] = f;
                    s[f][1] = SC;
                    s[f][2] = ST; //!!!!!!!!! здесь было s[f][3]
                    s[f][3] = SC - ST; //!!!!!!! здесь было s[f][4]
                }
                
                for (int f = 0; f < lstr; f++) {
                    for (int g = 0; g < 4; g++) {
                        cout << setw(4) << s[f][g] << "  ";
                        cout << endl;
                    }
                }
                
                for (j = lstr-1; j > 0; j--){
                    for (i = 0; i < j; i++) {
                        if (s[i+1][4] < s[i][4]) {
                            type* min = s[i];
                            s[i] = s[i+1];
                            s[i+1] = min;
                        } 
                    }
                }
                cout << "sorted data" << endl;
                cout << endl;
                cout << "  №  " << "|" << " S of circle " << "|" << " S of tringle " << "|" << "   SC - ST  " << endl;
                
                for (i = 0; i < lstr; i++) {
                    for (j = 0; j < 4; j++) {
                        cout << setw(4) << s[i][j] << "  ";
                    }
                    cout << endl;
                }
            }
        } 
    }
 
 
 
} 
 
 
int main()
{
    int lstr;
    cout << "Enter lstr strings" << endl;
    cin >> lstr;
    int q;
 
    do {     
        cout << "Choose the type of array" << endl;
        cout << "1 - int" << endl;
        cout << "2 - double" << endl;
        cout << "3 - float" << endl;
        cin >> q;
        
        switch(q) 
        {
        case 1: int ** s1;
                int ** a1;
                func ( lstr, a1, s1);
                break;
        
        case 2: double **s2;
                double **a2;
                func ( lstr, a2, s2);
                break;
        
        case 3: float **s3;
                float **a3;
                func ( lstr, a3, s3);
                break;
        }
    } while (q != 4);
    
    
    return 0;
}

#include "stdafx.h"
#include <iostream>
#include <iomanip>
#include <fstream>
#include <cstring>
#include <cstdlib>
#include <cmath>
 
using namespace std;
 
 
template <class type> 
void func (int lstr, type **a1, type **s1)
{
    int i, j;
    double const pi = 3.14;
    //type x1, x2, x3, y1, y2, y3, x, y, SC, ST;
    double ma, mb, x1, x2, x3, y1, y2, y3, x, y, SC, ST,r;
    
    type **s = new type*[lstr];
    for ( i = 0; i < lstr; i++){
        s[i] = new type[4];
    }
 
/*double **s = new double*[lstr];
    for ( i = 0; i < lstr; i++){
        s[i] = new double[4];
    }*/
   
    type **a = new type*[lstr];
    for ( i = 0; i < lstr; i++){ 
        a[i] = new type[2];
    }
    
    
    
    cout << "Enter matrix: ";
    
    for ( i = 0; i < lstr; i++) {
        for (j = 0; j < 2; j++) {
            cin >> a[i][j];
        }
        cout << endl;
    }
    
    cout << "  X  " << " | " << "  Y  " << endl;
    cout << endl;
    
    for (i = 0; i < lstr; i++) {
        for (j = 0; j < 2; j++){ 
            cout << setw(4) << a[i][j] << "  ";
        }
        cout << endl;
    }
    
    bool flag, flag2, flag3;
  
    
    for (int ii = 0; ii < lstr - 2; ii++){
         if (flag3 == false) {
            break;
        }
        for (int jj = ii + 1; jj < lstr - 1; jj++) {
            if (flag2 == false) {
            flag3 = false;
            break;
        }
            for (int kk = jj + 1; kk < lstr; kk++) { 
                if(flag == false) {
                    flag2 = false;
                    break;
                }
                for (int f = 0; f < lstr ; f++){ 
                    
                x1 = a[ii][0];
                y1 = a[ii][1];
 
                x2 = a[ii+1]/*[jj]*/[0];
                y2 = a[ii+1]/*[jj]*/[1];
  
                x3 = a[ii+2]/*[kk]*/[0];
                y3 = a[ii+2]/*[kk]*/[1];
 
                cout << " P1(x1;y1) =  ( " << x1 << " ; " << y1 << " )" << endl;
                cout << " P2(x2;y2) =  ( " << x2 << " ; " << y2 << " )" << endl;
                cout << " P3(x3;y3) =  ( " << x3 << " ; " << y3 << " )" << endl;
                
                /////////здесь будет ошибка в расчетах (бесконечность), если в знаменателе будет 0... а он там будет при равенстве x1 и x2 или x2 и x3
                if ((x1 != x2) && (x2 != x3)) {
                
                ma = (y2-y1)/(x2-x1);
                cout << "ma" <<  endl; cout <<ma <<endl;
                mb = (y3-y2)/(x3-x2);
                cout << "mb" << mb << endl;
                }
                else {
                    cout << "False! x1=x2 or x2=x3! Please enter new points" << endl;
                    flag = false;
                        break;}
                //////////////
 
                if (mb - ma == 0) {
                    cout << "P1P2 || P2P3 : these points don't form a circle " << endl;
                    break;
                }
            
 
                 x =  (ma*mb*(y1-y3)+mb*(x1+x2)-ma*(x2+x3))/2*(mb-ma);
                 y =  ((-1)*(x - (x1+x2)/2))/ma +(y1+y2)/2;
                
                cout << "The center O(x;y)" << x << "  " << y << endl; 
 
                //double sq = (x1-x)*(x1-x) - (y1-y)*(y1-y);
                /*if (sq < 0) {
                    cout << "Radius incorrect" << endl;
                    break;
                }*/
                //type r =(type) sqrt(double( /*(x1-x)*(x1-x) - (y1-y)*(y1-y)*/sq )); 
                r = sqrt(double (fabs(( (x1-x)*(x1-x) - (y1-y)*(y1-y))) )); 
                cout << "The radius " << r << endl;
                
                if (r == 0) {
                    cout << "these points don't form a circle"<< endl;
                    break;
                }
                
                SC = pi*r*r;
                cout << "Square of the circle " << SC << endl;
                
                ST = fabs(((x1-x3)*(y2-y3) - (x2-x3)*(y1-y3))/2.0);
                cout << "Square of the tringle " << ST << endl;
                
         
                    s[f][0] = ii;
                    s[f][1] = SC;
                    s[f][2] = ST; 
                    s[f][3] = SC - ST;
             
            }
        }
         
     }            
                
    }
    
   // как сделать проверку на НЕНАЛИЧИЕ хотя бы одного решения, дабы не выполнять сортировку?????
      
 cout << "  №  " << "|" << " S of circle " << "|" << " S of tringle " << "|" << "   SC - ST  " << endl;
   for (int f = 0; f < lstr; f++) {
                    for (int g = 0; g < 4; g++) {
                        cout << setw(4) << s[f][g] << "  ";
                       
                    } 
                    cout << endl;
                }
                
                for (j = lstr-1; j > 0; j--){
                    for (i = 0; i < j; i++) {
                        if (s[i+1][4] < s[i][4]) {
                            type* min = s[i];
                            s[i] = s[i+1];
                            s[i+1] = min;
                        } 
                    }
                }
                cout << "sorted data" << endl;
                cout << endl;
                cout << "  №  " << "|" << " S of circle " << "|" << " S of tringle " << "|" << "   SC - ST  " << endl;
                
                for (int iw = 0; iw < lstr; iw++) {
                    for (int jw = 0; jw < 4; jw++) {
                        cout << setw(4) << s[iw][jw] << "  ";
                    }
                    cout << endl;
                }
       } 
   }
  
 
 
int main()
{
    int lstr;
    cout << "Enter lstr > 3 strings" << endl;
    cin >> lstr;
    int q;
 
    do {     
        cout << "Choose the type of array" << endl;
        cout << "1 - int" << endl;
        cout << "2 - double" << endl;
        cout << "3 - float" << endl;
        cin >> q;
        
        switch(q) 
        {
        case 1: int ** s1;
                int ** a1;
                func ( lstr, a1, s1);
                break;
        
        case 2: double **s2;
                double **a2;
                func ( lstr, a2, s2);
                break;
        
        case 3: float double **s3;
                float **a3;
                func ( lstr, a3, s3);
                break;
        }
    } while (q != 4);
    
    
    return 0;
}

#include "stdafx.h"
#include <iostream>
#include <iomanip>
#include <fstream>
#include <cstring>
#include <cstdlib>
#include <cmath>
 
using namespace std;
 
 
template <class type> 
void func (int lstr, type **a1, type **s1)
{
    int i, j;
    double const pi = 3.14;
    //type x1, x2, x3, y1, y2, y3, x, y, SC, ST;
    double ma, mb, x1, x2, x3, y1, y2, y3, x, y, SC, ST,r;
 
    int count_points = 0; // количество нормальных точек
    
    type **s = new type*[lstr];
    for ( i = 0; i < lstr; i++){
        s[i] = new type[4];
    }
 
 
   
    type **a = new type*[lstr];
    for ( i = 0; i < lstr; i++){ 
        a[i] = new type[2];
    }
    
       
    cout << "Enter matrix:\n";
    
    for ( i = 0; i < lstr; i++) {
        for (j = 0; j < 2; j++) {
            cin >> a[i][j];
        }
        cout << endl;
    }
    
    cout << "  X  " << " | " << "  Y  " << endl;
    cout << endl;
    
    for (i = 0; i < lstr; i++) {
        for (j = 0; j < 2; j++){ 
            cout << setw(4) << a[i][j] << "  ";
        }
        cout << endl;
    }
    
    bool flag=true, flag2=true, flag3=true;
  
    
    for (int ii = 0; ii < lstr - 2; ii++){
        for (int jj = ii + 1; jj < lstr - 1; jj++) {
            for (int kk = jj + 1; kk < lstr; kk++) {                    
 
                x1 = a[ii][0];
                y1 = a[ii][1];
 
                x2 = a[jj][0];
                y2 = a[jj][1];
  
                x3 = a[kk][0];
                y3 = a[kk][1];
 
                cout<<endl<<endl;
                cout << " P1(x1;y1) =  ( " << x1 << " ; " << y1 << " )" << endl;
                cout << " P2(x2;y2) =  ( " << x2 << " ; " << y2 << " )" << endl;
                cout << " P3(x3;y3) =  ( " << x3 << " ; " << y3 << " )" << endl;
                cout << endl;
                
                /////////здесь будет ошибка в расчетах (бесконечность), если в знаменателе будет 0... а он там будет при равенстве x1 и x2 или x2 и x3
                if ((x1 != x2) && (x2 != x3)) {            
                    ma = (y2-y1)/(x2-x1);
                    cout << "ma = " <<ma <<endl;
                    mb = (y3-y2)/(x3-x2);
                    cout << "mb = " << mb << endl;
                }
                else {
                    cout << "False! x1=x2 or x2=x3!" << endl;
                    flag = false;
                    continue;
                    }
                //////////////
 
                if (mb - ma == 0) {
                    cout << "P1P2 || P2P3 : these points don't form a circle \n" << endl;
                    continue;
                }
            
                if(ma == 0) {
                    cout<<"Bad points"<<endl<<endl;
                    continue;
                }
 
                 x =  (ma*mb*(y1-y3)+mb*(x1+x2)-ma*(x2+x3))/2*(mb-ma);
                 y =  ((-1)*(x - (x1+x2)/2))/ma +(y1+y2)/2;
                
                 cout << "The center O(x;y)" << x << "  " << y << endl; 
 
 
                 r = sqrt(double (fabs(( (x1-x)*(x1-x) - (y1-y)*(y1-y))) )); 
                 cout << "The radius " << r << endl;
                
                 if (r == 0) {
                     cout << "these points don't form a circle"<< endl;
                     continue;
                 }
                
                 SC = pi*r*r;
                 cout << "Square of the circle " << SC << endl;
                
                 ST = fabs(((x1-x3)*(y2-y3) - (x2-x3)*(y1-y3))/2.0);
                 cout << "Square of the tringle " << ST << endl;
                
         
                 s[count_points][0] = ii;
                 s[count_points][1] = SC;
                 s[count_points][2] = ST; 
                 s[count_points][3] = SC - ST;
                 count_points++;
             
            } 
        }    
    }
    
        
    
   // как сделать проверку на НЕНАЛИЧИЕ хотя бы одного решения, дабы не выполнять сортировку?????
    if (count_points > 0)
    {
        cout<<endl<<endl;
        cout << "unsorted data" << endl;
        cout << "  №  " << "|" << " S of circle " << "|" << " S of tringle " << "|" << "   SC - ST  " << endl;
        for (int f = 0; f < count_points; f++) {
            for (int g = 0; g < 4; g++) {
                cout << setw(4) << s[f][g] << "  ";      
            } 
            cout << endl;
        }
                
        for (j = count_points-1; j > 0; j--){
            for (i = 0; i < j; i++) {
                if (s[i+1][4] < s[i][4]) {
                    type* min = s[i];
                    s[i] = s[i+1];
                    s[i+1] = min;
                } 
            }
        }
        cout<<endl;
 
        cout << "sorted data" << endl;
        cout << endl;
        cout << "  №  " << "|" << " S of circle " << "|" << " S of tringle " << "|" << "   SC - ST  " << endl;
                
        for (int iw = 0; iw < count_points; iw++) {
            for (int jw = 0; jw < 4; jw++) {
                cout << setw(4) << s[iw][jw] << "  ";
            }   
            cout << endl;
        }
 
        cout<<endl<<endl<<endl;
    
    } else cout<<"No good points"<<endl<<endl<<endl;
        
}
   
  
 
 
int main()
{
    int lstr;
    cout << "Enter lstr > 3 strings" << endl;
    cin >> lstr;
    int q;
 
    do {     
        cout << "Choose the type of array" << endl;
        cout << "1 - int" << endl;
        cout << "2 - double" << endl;
        cout << "3 - float" << endl;
        cin >> q;
        
        switch(q) 
        {
        case 1: int ** s1;
                int ** a1;
                func ( lstr, a1, s1);
                break;
        
        case 2: double **s2;
                double **a2;
                func ( lstr, a2, s2);
                break;
        
        case 3: double **s3;
                double **a3;
                func ( lstr, a3, s3);
                break;
        }
    } while (q != 4);
    
    
    return 0;
}