#include <iostream>
#include <cstdlib>

using namespace std;

int main()
{
    setlocale(0, "");
    int m1, n1, m2, n2;
    cout << "Введите размер первой матрицы: " << endl;;
    cin >> m1;
    n1=m1;
    m2=n2=n1;
    int a[m1][n1];
    int b[m2][n2];
    cout << endl;
    cout << "Матрица А:" << endl;
    for (int i=0; i<m1; i++)
    {
        for (int j=0; j<n1; j++)
        {
            a[i][j]=rand()%10+1;
            cout << a[i][j] << "\t";
        }
        cout << endl;
    }
    cout << endl << "Матрица В:" << endl;
    for (int i=0; i<m2; i++)
    {
        for (int j=0; j<n2; j++)
        {
            b[i][j]=rand()%10+1;
            cout << b[i][j] << "\t";
        }
        cout << endl;
    }
    int c[m1][n2];
    int b3[m1][m1];
    for (int i=0; i<m2; i++) //3B
    {
        for (int j=0; j<n2; j++)
        {
            b3[i][j]=3*b[i][j];
        }
    }

    for (int row = 0; row < m1; row++) //A*B
    {
        for (int col = 0; col < n2; col++)
        {
            c[row][col]=0;
            for (int inner = 0; inner < n1; inner++)
            {
                c[row][col] += a[row][inner] * b[inner][col];
            }
        }
    }
    cout << endl;
    int d[m1][n1];
    for (int i=0; i<m1; i++) //3B*A*B
    {
        for (int j=0; j<n1; j++)
        {
            d[i][j]=0;
            for (int k=0; k<n2; k++)
            {
                d[i][j]+=c[i][k]*a[k][j];
            }
        }
    }
    int e[m1][n1];
    for (int i=0; i<m1; i++) //B*B
    {
        for (int j=0; j<m1; j++)
        {
            e[i][j]=0;
            for (int k=0; k<m1; k++)
            {
                e[i][j]+=b[i][k]*b[k][j];
            }
        }
    }
    for (int i=0; i<m1; i++) //2A
    {
        for (int j=0; j<m1; j++)
        {
            a[i][j]*=2;
        }
    }
    for (int i=0; i<m1; i++) //2A*(B*B)
    {
        for (int j=0; j<m1; j++)
        {
            b[i][j]=0;
            for (int k=0; k<m1; k++)
            {
                b[i][j]+=a[i][k]*e[k][j];
            }
        }
    }
    int f[m1][m1];
    for (int i=0; i<m1; i++)
    {
        for (int j=0; j<m1; j++)
        {
            f[i][j]=0;
        }
    }
    cout << "Матрица F=3B*A*B-2A*(B*B)" << endl;
    for (int i=0; i<m1; i++) //3B*A*B-2A*(B*B)
    {
        for (int j=0; j<m1; j++)
        {
            f[i][j]=d[i][j]-b[i][j];
            cout << f[i][i] << "\t";
        }
        cout << endl;
    }
    cin.get();
    cin.get();
    return 0;
}