Упорядочить строки матриц по числу элементов кратных заданному числу

#include <iostream>
#include <fstream>
void  m_qsort(int** m, int* cs, int M, int l, int r);
template<class Cmp>
bool  matrix_sort(int** m, int N, int M, Cmp cmp);
int** matrix_input(std::istream& _in, int& N, int& M);
void  matrix_output(std::ostream& _out, int** m, int N, int M);
void  matrix_free(int** m, int N);
 
struct icmp {
    int n;
    icmp(int _n):n(_n){}
 
    bool operator () (int x) const {
        return ((x % n) == 0);
    }
};
 
int main(void){
    int** m, N, M;
/*  ввод из файла
    std::ifstream fp("matrix.txt");
    m = matrix_input(fp, N, M);
    fp.close();
    if(m == NULL)
        return 1;
    matrix_output(std::cout, m, N, M);
    matrix_sort(m, N, M, icmp(3));
    matrix_output(std::cout, m, N, M);
    matrix_free(m, N);
*/
    m = matrix_input(std::cin, N, M);
    if(m == NULL)
        return 1;
    matrix_output(std::cout, m, N, M);
    
    //сортируем строки матрицы по-числу кратному 3-ём
    matrix_sort(m, N, M, icmp(3));
    matrix_output(std::cout, m, N, M);
    matrix_free(m, N);
    return 0;
}
 
//сортировка строк матрицы
template<class Cmp>
bool matrix_sort(int** m, int N, int M, Cmp cmp){
    int* cs = new (std::nothrow) int[N];
    if(cs == NULL)
        return false;
 
    for(int i = 0; i < N; ++i){
        cs[i] = 0;
        for(int j = 0; j < M; ++j){
            if(cmp(m[i][j]))
                ++cs[i];
        }
    }
    m_qsort(m, cs, M, 0, N - 1);
    delete[] cs;
    return true;
}
 
//ввод матрицы из входного потока
int** matrix_input(std::istream& _in, int& N, int& M){
    int** m, i, j;
    if(!(_in >> N >> M) || (N <= 1) || (M <= 1))
        return NULL;
 
    m = new (std::nothrow) int*[N];
    if(m == NULL)
        return NULL;
 
    for(i = 0; i < N; ++i){
        m[i] = new (std::nothrow) int[M];
        if(m[i] == NULL){
            --i;
            goto err;
        }
 
        j = 0;
        while((j < M) && (_in >> m[i][j]) && !_in.fail())
            ++j;
 
        if(j != M)
            goto err;
    }
    return m;
err:
    for(j = i; j >= 0; --j)
        delete[] m[j];
    delete[] m;
    return NULL;
}
 
//вывод матрицы в выходной поток
void matrix_output(std::ostream& _out, int** m, int N, int M){
    for(int i = 0; i < N; ++i){
        for(int j = 0; j < M; ++j)
            _out << m[i][j] << ' ';
        _out << std::endl;
    }
    _out << std::endl;
}
 
//удаление матрицы
void matrix_free(int** m, int N){
    for(int i = 0; i < N; ++i)
        delete[] m[i];
    delete[] m;
}
 
//базовая быстрая сортировка
void m_qsort(int** m, int* cs, int M, int l, int r){
    int p, i = l, j = r;
    if(l >= r)
        return;
    
    p = cs[l + (r - l) / 2];
    while(i <= j){
        while(cs[i] < p)
            ++i;
        while(cs[j] > p)
            --j;
 
        if(i <= j){
            for(int v = 0; v < M; ++v)
                std::swap(m[i][v], m[j][v]);
            std::swap(cs[i], cs[j]);
            ++i;
            --j;
        }
    }
 
    if(j > l)
        m_qsort(m, cs, M, l, j);
    if(r > i)
        m_qsort(m, cs, M, i, r);
}