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

#include <iostream>
#include <random>
#include <vector>
#include <algorithm>
#include <iomanip>
 
using namespace std;
 
int main()
{
    setlocale(LC_ALL, "Rus");
    random_device rd;
    mt19937 g{ rd() };
    uniform_real_distribution<> dist(-10., 10.);
    auto random = [&g, &dist]() { return dist(g); };
 
    cout << "Введите сторону матрицы n = ";
    size_t n;
    cin >> n;
    vector<vector<double>> a(n, vector<double>(n));
    cout << "Исходное содержание матрицы" << endl;
    double p = 1.;
    for (auto& line : a)
    {
        for (auto& value : line)
        {
            cout << setw(8u) << fixed << setprecision(3u) << (value = random());
            p *= value;
        }
        cout << endl;
    }
    cout << "Произведение всех элементов матрицы " << p
        << " или " << scientific << p << endl;
    cout << "Преобразованное содержание матрицы" << endl;
    for (size_t row = 0u; row < n; ++row)
    {
        a[row][row] = p;
        a[row][n - row - 1] = p;
        for (const auto& value : a[row])
        {
            cout << setw(15u) << scientific << value;
        }
        cout << endl;
    }
 
    system("pause");
    return 0;
}