Поиск седловой точки матрицы

#include "stdafx.h"
#include <iostream>
#include <limits>
 
using std::endl;
using std::cout;
using std::pair;
using std::make_pair;
using std::numeric_limits;
 
template<typename T> class Matrix {
private:
    T * *matrix;
    size_t N, M;
public:
    Matrix(size_t _N, size_t _M) : N(_N), M(_M) {
        this->matrix = new T*[this->N];
        for (size_t i(0); i < this->N; i++) {
            this->matrix[i] = new T[this->M];
        }
    };
    T* operator[](const size_t index) { return this->matrix[index]; }
 
    friend std::ostream& operator<<(std::ostream &os, Matrix &obj) {
        for (size_t i(0); i < obj.N; i++) {
            for (size_t j(0); j < obj.M; j++) {
                os << obj[i][j];
                if (j + 1 != obj.M)
                    os << " ";
            }
            os << std::endl;
        }
        return os;
    }
 
    void getSaddlePoint() {
        T* min = new T[N];
        T* max = new T[M];
        for (size_t i = 0; i < N; i++)
        {
            min[i] = numeric_limits<T>::max();
            for (size_t j = 0; j < M; j++)
            {
                if (matrix[i][j] < min[i])
                {
                    min[i] = matrix[i][j];
                }
            }
        }
        for (size_t j = 0; j < M; j++)
        {
            max[j] = numeric_limits<T>::min();
            for (size_t i = 0; i < N; i++)
            {
                if (matrix[i][j] > max[j])
                {
                    max[j] = matrix[i][j];
                }
            }
        }
        cout << "Saddle point indexes" << endl;
        size_t k = 0;
        for (size_t i = 0; i < N; i++)
        {
            for (size_t j = 0; j < M; j++)
            {
                if (min[i] == max[j])
                {
                    cout << i << " " << j << endl;
                    k++;
                }
            }
        }
        if (!(k))
            cout << "not Saddle points!" << endl;
 
        delete[] min;
        delete[] max;
    }
 
 
    ~Matrix(void) {
        for (size_t i(0); i < this->N; i++)
            delete this->matrix[i];
        delete this->matrix;
    };
};
/*-------------------------------------------------------------------------*/
 
int main() {
    cout << "TEST 1";
    const size_t N(2), M(3);
    Matrix<int> matrix(N, M);
    //in
    matrix[0][0] = 3;
    matrix[0][1] = 2;
    matrix[0][2] = 1;
    matrix[1][0] = 1;
    matrix[1][1] = 3;
    matrix[1][2] = 4;
    //out
    cout << endl << matrix << endl;
    //
    cout << "->Correct answer : not points";
    cout << endl << "->CHECK:" << endl;
    matrix.getSaddlePoint();
 
    cout << endl << "TEST 2";
    const size_t N2(3), M2(4);
    Matrix<int> matrix2(N2, M2);
    //in
    matrix2[0][0] = 2;
    matrix2[0][1] = 3;
    matrix2[0][2] = 5;
    matrix2[0][3] = 2;
    matrix2[1][0] = 2;
    matrix2[1][1] = 4;
    matrix2[1][2] = 6;
    matrix2[1][3] = 2;
    matrix2[2][0] = -2;
    matrix2[2][1] = 7;
    matrix2[2][2] = 2;
    matrix2[2][3] = 0;
    //out
    cout << endl << matrix2 << endl;
    //
    cout << "->Correct answer : four  saddle point: (0, 0), (0, 3), (1, 0), (1, 3)";
    cout << endl << "->CHECK:" << endl;
    matrix2.getSaddlePoint();
 
    system("pause");
    return 0;
}