#include <iostream>
#define N 4
using namespace std;
void Square_Matrix_Multiply(int mat[N][N],int m1[N][N],int m2[N][N])
{
for (int i = 0;i < N ;++i)
{
for (int j = 0;j < N;++j)
{
mat[i][j] = 0;
for (int k = 0; k < N; ++ k)
{
mat[i][j] += m1[i][k]*m2[k][j];
}
}
}
}
void Input(int mat[N][N])
{
for( int i = 0; i < N; i++ )
{
for( int j = 0; j < N; j++ )
{
cin>>mat[i][j];
}
}
}void Output(int mat[N][N])
{
for(int i = 0; i < N; i++ )
{
for( int j = 0; j < N; j++ )
{
cout << mat[i][j]<<" ";
}
cout<<"\n";
}
}
//matrix + matrix
void plus( int t[N/2][N/2], int r[N/2][N/2], int s[N/2][N/2] )
{
int i, j;
for( i = 0; i < N / 2; i++ )
{
for( j = 0; j < N / 2; j++ )
{
t[i][j] = r[i][j] + s[i][j];
}
}
}
//matrix - matrix
void minus( int t[N/2][N/2], int r[N/2][N/2], int s[N/2][N/2] )
{
int i, j;
for( i = 0; i < N / 2; i++ )
{
for( j = 0; j < N / 2; j++ )
{
t[i][j] = r[i][j] - s[i][j];
}