Удалить наибольший элементов каждого столбца матрицы

#include "stdafx.h"
#include<iostream>
#include<ctime>
#include<conio.h>
 
using namespace std;
 
int searchmax(int n, int **a, int j)
{
    int max;
    int k;
    max = **a;
    for (int i = 0; i < n + 1;i++)
    {
        if (*(*(a + i) + j) > max)
        {
            max = *(*(a + i) + j);
            k = i;
 
        }
    }
    return(k);
}
void sdvig(int n, int **a, int k, int j)
{
    for (int i = k;i < n + 1;i++)
        *(*(a + i) + j) = *(*(a + i + 1) + j);
}
void main()
{
    srand(time(NULL));
    int n = 0;
    cin >> n;
    int **a = new int*[n + 1];
    for (int i = 0; i < n; i++)
    {
        a[i] = new int[n + 1];
    }
    for (int i = 0; i < n; i++)
    {
        for (int j = 0; j < n; j++)
        {
            *(*(a + i) + j) = rand() % 10;
            cout << *(*(a + i) + j) << " ";
        }
        cout << "\n";
    }
 
    cout << "\n";
    int k = 0;
    for (int j = 0;j < n;j++)
    {
        k = searchmax(n, a, j);
        cout << " max in " << j << " stolbze: " << *(*(a + k) + j);
 
        sdvig(n, a, k, j);
 
    }
    cout << "\n";
    for (int i = 0; i < n; i++)
    {
        for (int j = 0; j < n; j++)
        {
            *(*(a + i) + j) = rand() % 10;
            cout << *(*(a + i) + j) << " ";
        }
        cout << "\n";
    }
    for (int i = 0;i < n;i++)
        delete[] a[i];
    _getch();
}

int max;
...
max = **a;

#include <iostream>
#include <iomanip>
#include <cstdlib>
#include <ctime>
 
int max_column_index(int ** matrix, int rows, int col) {
    int mi = 0;
    
    for ( int i = 1; i < rows; ++i )
        if ( matrix[i][col] > matrix[mi][col] )
            mi = i;
    
    return mi;
}
 
void dump_matrix(int ** matrix, int rows, int columns) {
    for ( int i = 0; i < rows; ++i ) {
        for ( int j = 0; j < columns; ++j ) 
            std::cout << std::setw(3) << matrix[i][j];
        std::cout << std::endl;
    }
}
 
int main() {
    int n;
    
    std::cout << "N = ";
    std::cin >> n;
    
    int ** m1 = new int* [ n + 1 ];
    for ( int i = 0; i < n + 1; ++i )
        m1[i] = new int [ n ];
    
    srand(time(NULL));
    
    for ( int i = 0; i < n + 1; ++i )
        for ( int j = 0; j < n; ++j )
            m1[i][j] = rand() % 100;
    
    int ** m2 = new int* [ n ];
    for ( int i = 0; i < n; ++i )
        m2[i] = new int [ n ];
    
    for ( int j = 0; j < n; ++j ) {
        int mi = max_column_index(m1, n + 1, j);
        int ii = 0;
        
        for ( int i = 0; i < n + 1; ++i ) {
            if ( i == mi )
                continue;
            m2[ii++][j] = m1[i][j];
        }
    }
    
    std::cout << "First matrix:" << std::endl;
    dump_matrix(m1, n + 1, n);
    
    std::cout << "\nSecond matrix:" << std::endl;
    dump_matrix(m2, n, n);
    
    for ( int i = 0; i < n + 1; ++i )
        delete [] m1[i];
    delete [] m1;
    
    for ( int i = 0; i < n; ++i )
        delete [] m2[i];
    delete [] m2;
    
    return 0;
}

#include <iostream>
#include <ctime>
using namespace std;
 
int main()
{       
   srand(time(0)); // генерация случайных чисел
   int N, max=0, imax;
   cout <<"N="; cin >>N;   
   
// объявление двумерных динамических массивов
 
   int **a = new int* [N+1]; 
   for (int i = 0; i < N+1; i++)
        a[i] = new int [N];
        
   int **b = new int* [N]; 
   for (int i = 0; i < N; i++)
        b[i] = new int [N];    
        
    
// заполнение массива случайными числами от 1 до 9 и вывод на экран
   
    cout <<"Matrix 1:\n"; 
    for (int i = 0; i < N+1; i++)
        {
        for (int j = 0; j < N; j++)
        {
            a[i][j] = rand() %9 + 1;
            cout << a[i][j] << " ";            
        }        
            cout <<endl;            
        }
        
// определение последнего max столбца и сдвиг элементов
 
      for (int j = 0; j < N; j++)
        {
        for (int i = 0; i < N+1; i++)            
            if (i==0 ||a[i][j] >= max) 
            {
            max = a[i][j];
            imax=i;
            }
            for (int k = imax; k < N; k++)
            a[k][j]=a[k+1][j];              
        }
        cout <<endl;
        
// заполнение и вывод нового массива
 
        cout <<"Matrix 2:\n";
        for (int i = 0; i < N; i++)
        {
        for (int j = 0; j < N; j++) 
        {
        b[i][j]=a[i][j];
        cout << b[i][j] << " ";
        }
        cout <<endl;            
        }        
        
// освобождение памяти отведённой под массивы:
 
    for (int i = 0; i < N+1; i++)
    delete [] a[i];
    delete [] a;
    for (int i = 0; i < N; i++)
    delete [] b[i];
    delete [] b; 
system("pause");
return 0;
}

int k = 0;
max = **a;