Написать программу нахождения суммы диагональных элементов двух матриц

#include <iostream>
#include <iomanip>
#include <cstdlib>
#include <ctime>
 
using namespace std;
 
int main(){
    
    srand(time(0));
    const int N = 10;
    int n;
    do{
        cout << "n = ";
        cin >> n;
    }while(n < 1 || N < n);
    cout << endl;
    
    int a[N][N], b[N][N], c[N];
    for (int i = 0; i < n; ++i){
        for (int j = 0; j < n; ++j){
            a[i][j] = rand()%100;
            b[i][j] = rand()%100;
        }
        
        c[i] = a[i][i] + b[i][i];
        
        for (int j = 0; j < n; ++j)
            cout << setw(3) << a[i][j];
        cout << "   |  ";   
        for (int j = 0; j < n; ++j)
            cout << setw(3) << b[i][j];
        
        cout << endl;
    }
    cout << endl;
    
    for (int i = 0; i < n; ++i)
        cout << setw(4) << c[i];
    cout << endl;
}