Перестановкой строк и столбцов получить все возможные варианты матрицы

#include <iostream>
 
using namespace std;
 
int my_permutation(int a[], int n, int k){
    //âûäàåò k-þ ïåðåñòàíîâêó áåç ïîâòîðåíèé, âñåãî n! øòóê
    
    int fact = 1;
    for (int i = 2; i <= n; ++i)
        fact *= i;
        
    int p;
    if (k >= fact || k < 0)
        p = 0;
    else
        p = 1;
    
    int temp[n];
    for (int i = 0; i < n; ++i)
        temp[i] = i;
    
    for (int i = 0; i < n - 1; ++i){
        fact /= (n - i);
        a[i] = temp[k / fact];
        
        //óäàëÿþ èç temp ýëåìåíò ñäâèãîì îñòàëüíûõ ýë-â âïåðåä
        for (int j = k / fact; j < n - i - 1; ++j)
            temp[j] = temp[j + 1];
        
        k = k % fact;
    }
    
    a[n - 1] = temp[0];
    
    return p;
}
 
int main(){
    
    int a[4];
    for (int i = 0; i < 4 * 3 * 2; ++i){
        my_permutation(a, 4, i);
        for (int j = 0; j < 4; ++j)
            cout << a[j] << ' ';
        cout << endl;
    }   
}

#include <iostream>
#include <iomanip>
 
using namespace std;
 
int my_permutation(int a[], int n, int k){
    //âûäàåò k-þ ïåðåñòàíîâêó áåç ïîâòîðåíèé, âñåãî n! øòóê
    
    int fact = 1;
    for (int i = 2; i <= n; ++i)
        fact *= i;
        
    int p;
    if (k >= fact || k < 0)
        p = 0;
    else
        p = 1;
    
    int temp[n];
    for (int i = 0; i < n; ++i)
        temp[i] = i;
    
    for (int i = 0; i < n - 1; ++i){
        fact /= (n - i);
        a[i] = temp[k / fact];
        
        //óäàëÿþ èç temp ýëåìåíò ñäâèãîì îñòàëüíûõ ýë-â âïåðåä
        for (int j = k / fact; j < n - i - 1; ++j)
            temp[j] = temp[j + 1];
        
        k = k % fact;
    }
    
    a[n - 1] = temp[0];
    
    return p;
}
 
int main(){
    
    int n, m;
    cin >> n >> m;
    
    cout << "-----------------------------------------" << endl;
    int matrix[10][10];
    for (int i = 0; i < n; ++i){
        for (int j = 0; j < m; ++j){
            matrix[i][j] = i * m + j;
            cout << setw(3) << matrix[i][j];
        }
        cout << endl;
    }
    cout << "-----------------------------------------" << endl;
    cout << "-----------------------------------------" << endl;
        
    int f_n = 1;
    for (int i = 2; i <= n; ++i)
        f_n *= i;
        
    int f_m = 1;
    for (int i = 2; i <= m; ++i)
        f_m *= i;
        
    for (int i = 0; i < f_n; ++i){
        int a[n];
        my_permutation(a, n, i);
        for (int j = 0; j < f_m; ++j){
            int b[m];
            my_permutation(b, m, j);
            
            for (int i1 = 0; i1 < n; ++i1){
                for (int j1 = 0; j1 < m; ++j1)
                    cout << setw(3) << matrix[a[i1]][b[j1]];
                cout << endl;
            }
            cout << "-----------------------------------------" << endl;
        }       
    }       
}