Задана матрица размером NxM. Определить количество "особых" элементов матрицы

// requirements: compiler that supports C++11
 
#ifndef MATRIX_H
#define MATRIX_H
 
#include <iostream>
#include <vector>
 
typedef double Item;
 
class Matrix
{
private:
    int m_row;
    int m_col;
    Item ** matrix;
 
public:
    Matrix(int row, int col);
    ~Matrix();
 
    void fill_rand(Item fromNum, Item toNum);
    Item find_sum_of_col(const int col_index);
    Item & at(const int row_index, const int col_index);
    bool is_special(const Item & value, const int col_index);
 
    friend std::ostream & operator<<(std::ostream & os, const Matrix & m);
};
 
#endif // MATRIX_H

#include "Matrix.h"
#include <random>
#include <ctime>
#include <iomanip>
 
Matrix::Matrix(int row, int col)
{
    m_row = row;
    m_col = col;
 
    matrix = new Item * [m_row];
    for (int i = 0; i < m_row; ++i)
        matrix[i] = new Item[m_col];
}
 
Matrix::~Matrix()
{
    for (int i = 0; i < m_row; ++i)
        delete [] matrix[i];
    delete [] matrix;
}
 
void Matrix::fill_rand(Item fromNum, Item toNum)
{
    std::mt19937 gen(time(0));
    std::uniform_real_distribution<double> urd(fromNum, toNum);
    for (int i = 0; i < m_row; ++i)
        for (int j = 0; j < m_col; ++j)
            matrix[i][j] = urd(gen);
}
 
Item Matrix::find_sum_of_col(const int col_index)
{
    Item sum = 0;
    for (int i = 0; i < m_row; ++i)
        sum += matrix[i][col_index];
    return sum;
}
 
Item & Matrix::at(const int row_index, const int col_index)
{
    return matrix[row_index][col_index];
}
 
bool Matrix::is_special(const Item & value, const int col_index)
{
    return value > find_sum_of_col(col_index) - value;
}
 
std::ostream & operator<<(std::ostream & os, const Matrix & m)
{
    for (int i = 0; i < m.m_row; ++i)
    {
        for (int j = 0; j < m.m_col; ++j)
            os << std::setprecision(3) << m.matrix[i][j] << "\t";
        os << "\n";
    }
 
    return os;
}

#include <iostream>
#include "Matrix.h"
 
int main()
{
    const int rsize1 = 3;
    const int csize1 = 4;
    Matrix matrix(rsize1, csize1);
    matrix.fill_rand(0.11, 9.99);
 
    std::cout << "\nThis is the matrix:\n" << matrix;
 
    int special_element = 0;
    for (int i = 0; i < rsize1; ++i)
        for (int j = 0; j < csize1; ++j)
            if (matrix.is_special(matrix.at(i, j), j))
                special_element++;
 
    std::cout << "\nTotal amount of special elements is: " << special_element;
 
    return 0;
}