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

#include <iostream>
#include <iomanip>
#include <locale.h>
#include <cstdlib> 
#include <cmath> 
#include <time.h>
#define USE MATH DEFINES
using namespace std;
int main()
{
    setlocale(0, "");
    srand(time(NULL) * 1000);
    int** matrix = NULL;
    int N;
    int M;
    cout << "Введите количество строк N:";
    cin >> N;
    cout << "Введите количество столбцов M:";
    cin >> M;
    // описание массива и выделение памяти
    matrix = new int* [N];
    int row, col, a, b, k, p;
    for (row = 0; row < N; row++)
        matrix[row] = new int[M];
    // заполнение массива
    for (row = 0; row < N; row++)
    {
        for (col = 0; col < M; col++)
        {
            matrix[row][col] = rand() % 10 + 1;
        }
    }
    // вывод массива на экран
    cout << "Matrix" << endl;
    for (row = 0; row < N; row++)
    {
        for (col = 0; col < M; col++)
        {
            cout << setw(3) << matrix[row][col];
        }
        cout << endl;
    }
    cout << "Введите номера столбцов a , b " << endl;//вводим номера столбцов которые нужно поменять 
    cin >> a >> b;
    a -= 1, b -= 1;
    cout << endl;
    for (row = 0; row < N; row++)
        for (col = 0; col < M; col++)
        {
            int temp = matrix[row][a]; // временному элементу присваиваем  элемент a
            matrix[row][a] = matrix[row][b]; // элементу a элемент b
            matrix[row][b] = temp;        // b временный
        }
 
    for (row = 0; row < N; row++)
        for (col = 0; col < M; col++)
        {
            if (col % M == 0)
                cout << endl;
            cout << setw(3) << matrix[row][col]; // выводим изменённую матрицу
 
        }
    cout << endl;
    col = a;
    k = 1;
    for (row = 0; row < N; row++)
    {
        k = k * matrix[row][col];
    }
    cout << k << endl;
    col = b;
    p = 1;
    for (row = 0; row < N; row++)
    {
        p = p * matrix[row][col];
    }
    cout << p << endl;
    cout << "Сумма произведений a * b= " << k * p << endl;//выводим суму произведения элементов столбцов a,b
    // высвобождение памяти
    for (row = 0; row < N; row++)
    {
        delete[]matrix[row];
        matrix[row] = NULL;
    }
    delete[]matrix;
    matrix = NULL;
    system("PAUSE");
    return 0;
}

#include <iostream>
#include <cstdlib>  
#include <ctime>
using namespace std;
 
int main()
{
    srand((int)time(0));    
    int N, M, a, b; 
    long long mult_a=1LL, mult_b=1LL;
    cout << "N="; cin >> N;
    cout << "M="; cin >> M;    
    
    int** matrix = new int* [N];    
    for (int i = 0; i < N; i++)
        matrix[i] = new int[M];        
    
    cout << "Matrix:\n";
    for (int i = 0; i < N; i++)
    {
        for (int j = 0; j < M; j++)
        {
            matrix[i][j]=rand()%9 + 1;
            cout << matrix[i][j] << " ";
        }
    cout << "\n";    
    }
   
    cout << "Number cols a, b:\n"; 
    cin >> a >> b;
    a--; b--;
    
    for (int i = 0; i < N; i++)
        swap(matrix[i][a],matrix[i][b]);        
       
    cout << "Modified matrix:\n";
    for (int i = 0; i < N; i++)
    {
        mult_a*=matrix[i][a];
        mult_b*=matrix[i][b];
        for (int j = 0; j < M; j++)                    
            cout << matrix[i][j] << " ";             
    cout << "\n";
    }    
 
    cout << "mult_a + mult_b = " << mult_a + mult_b << "\n";
    
    for (int i = 0; i < N; i++)
    delete[]matrix[i];
    delete[]matrix;    
system("pause");
return 0;
}