Найти сумму всех диагоналей и сравнить между с собой

#include <iostream>
#include <ctime>
#include <limits>
using namespace std;
#define SIZE 4
int main(){
    srand(time(NULL));
    int a[SIZE][SIZE];
    
    for(int i = 0; i < SIZE; i++) {
        for(int j = 0; j < SIZE; j++) {
            a[i][j] = rand() % 100;
            cout << a[i][j] << "\t";
        }
        cout << endl;
    }
    
    int sumDiags[2 * SIZE - 1] = {};
    int k = 0;
    
    // с левого нижнего края до главной диагонали
    for(int i = 0; i < SIZE; i++) {
        for(int j = i; j >= 0; j--) {
            sumDiags[k] += a[SIZE - i - 1 + j][j];
        }
        ++k;
    }
    
    // после главной диагонали
    for(int i = 0; i < SIZE - 1; i++) {
        for(int j = SIZE - 1; j >= i + 1; j--) {
            sumDiags[k] += a[j - i - 1][j];
        }
        ++k;
    }
    
    int max = std::numeric_limits<int>::min();
    cout << endl << "суммы диагоналей : ";
    for(int i = 0; i < k; i++) {
        if(max < sumDiags[i]) max = sumDiags[i];
        cout << sumDiags[i] << " ";
    }
    cout << endl << "максимум: " << max;
    
    return 0;
}

#define N 4
#define M 4
 
int main()
{
    int matr[N][M] = {
       {1, 3, 4, 5},
       {2, 3, 4, 5},
       {2, 6, 3, 2},
       {4, 2, 4, 1}
    };
    int diagonals[N + M - 1] = {0}, maxDiag = 0;
 
    for( int i = N - 1; i > -1; i-- )
        for( int j = 0; j < M; j++)
            diagonals[N - 1 - i + j] += matr[i][j];
 
    for( int i = 0; i < N + M - 1; i++ )
        if( diagonals[i] > diagonals[maxDiag] )
            maxDiag = i;
 
    cout << "max diag: " << maxDiag + 1 << "\nmax sum : " << diagonals[maxDiag] << endl;
 
    return 0;
}