Переставить строки и столбцы квадратной матрица А, чтобы алгебраические суммы произведений элементов

#include <iostream>
 
using namespace std;
 
void stolb(int **mas,int n,int a,int b);       
void strok(int **mas,int n,int a,int b);        
 
main()
{
    srand(time(0));
    cout<<"Enter n : ";
    int n;
    cin>>n
    int **mas=new int*[n];
    for(int i=0;i<n;i++)
        mas[i]=new int[n];
    for(int i=0;i<n;i++)
    {
        for(int j=0;j<n;j++)
        {
            mas[i][j]=rand()%10+1;
            cout.width(5);
            cout<<mas[i][j]<<" ";
        }
        cout<<endl;
    }
    int max=0;
    int summa;
    int is,js;
    for(int z=0;z<n;z++)
    {
        is=0;
        js=0;
        max=0;
        for(int i=z;i<n;i++)
        {
            for(int j=z;j<n;j++)
            {
                summa=0;
                for(int d=0;d<n;d++)
                {
                    summa=mas[i][d]*mas[d][j]+summa; 
                }
                if(max<summa) 
                    {
                        max=summa;
                        is=i;
                        js=j;
                    }
            }
 
        }
        stolb(mas,n,is,z);  
        strok(mas,n,js,z);  
    }
    cout<<"\n\n\n";
    for(int i=0;i<n;i++)
    {
        for(int j=0;j<n;j++)
        {
            cout.width(5);
            cout<<mas[i][j]<<" ";
        }
        cout<<endl;
    }
    for(int i=0;i<n;i++)
        delete [] mas[i];
    delete [] mas;
    return 0;
}
void stolb(int **mas,int n,int a,int b) 
{
    int zamena;
    for(int i=0;i<n;i++)
    {
        zamena=mas[i][a];
        mas[i][a]=mas[i][b];
        mas[i][b]=zamena;
    }
}
void strok(int **mas,int n,int a,int b) 
{
    int zamena;
    for(int i=0;i<n;i++)
    {
        zamena=mas[a][i];
        mas[a][i]=mas[b][i];
        mas[b][i]=zamena;
    }
}