В заданной матрице найти максимальную сумму элементов прямоугольной подматрицы среди всех возможных подматриц

#include <ctime>
#include <cstdlib>
#include <iostream>
#include <iomanip>
 
using namespace std;
 
const size_t NRow = 5, NCol = 7;
 
int a[NRow][NCol];
 
bool subMatrixSum(size_t i0, size_t j0, size_t n, size_t m, int &sum)
{
    if (i0+n > NRow || j0+m > NCol) return false;
 
    sum = 0;
    for (size_t i = i0; i < i0+n; i++)
        for (size_t j = j0; j < j0+m; j++)
            sum += a[i][j];
 
    return true;
}
 
void printSubmatrix(size_t i0, size_t j0, size_t n, size_t m)
{
    for (size_t i = i0; i < i0+n; i++)
    {
        for (size_t j = j0; j < j0+m; j++)
            cout << setw(3) << a[i][j] << ' ';
        cout << '\n';
    }
}
 
int main()
{
    srand(time(NULL));
    for (size_t i=0; i<NRow; i++)
        for (size_t j=0; j<NCol; j++)
            a[i][j] = -50 + rand() % 100;
 
    cout << "Generated matrix:\n";
    printSubmatrix(0, 0, NRow, NCol);
 
    size_t iMax=0, jMax=0, nMax=1, mMax=1;
    int sMax=a[0][0];
 
    for (size_t n=1; n<=NRow; n++)
        for (size_t m=1; m<=NCol; m++)
            for (size_t i=0; i<NRow; i++)
                for (size_t j=0; j<NCol; j++)
                {
                    int s;
                    bool fit = subMatrixSum(i, j, n, m, s);
                    if (fit && s > sMax)
                    {
                        iMax = i;
                        jMax = j;
                        nMax = n;
                        mMax = m;
                        sMax = s;
                    }
                }
 
    cout << "\nSubMatrix geometry: { (" << iMax << ", " << jMax << "), "
         << nMax << "x" << mMax << " }\n";
    cout << "Element sum: " << sMax << endl;
    cout << "\nSubmatrix:\n";
    printSubmatrix(iMax, jMax, nMax, mMax);
}