#include <iostream>
#include <stdlib.h>


using namespace std;

class matr
{
	int **ms;
	int x, y;
	public:
		matr (int);
		friend void matr_read (matr &);
		friend void matr_print (matr &);
		matr operator& (matr);
		matr operator+ (matr);
        matr (matr &);
};

matr :: matr (int std_size=10)
{
    ms = new int*[std_size];
    for (int i=0; i<std_size; i++)
        ms[i] = new int[std_size];
}

matr matr :: operator& (matr mtr2)
{
    int tmp;
    matr temp;
	if (y == mtr2.x)
	{
		temp.x = mtr2.x;
		temp.y = mtr2.y;
		for (int i = 0; i < x; i++)
			for (int j = 0; j < y; j++)
			{
				tmp = 0;
				for (int k = 0; k < y; k++)
					tmp += (ms[i][k]) * (mtr2.ms[k][j]);
				temp.ms[i][j] = tmp;
			}
	}
	else return temp;
    return temp;
}

matr matr :: operator+ (matr mtr2)
{
    matr temp;
    if (x != mtr2.x || y != mtr2.y)
		return temp;
	temp.x = x;
	temp.y = y;
	for (int i = 0; i < x; i ++)
		for (int j = 0; j < y; j ++)
			temp.ms[i][j] = ms[i][j] + mtr2.ms[i][j];

    return temp;
}

matr :: matr (matr &obj)
{
    ms = new int*[10];
    x=obj.x;
    y=obj.y;
    for (int i=0; i<x; i++)
    {
        ms[i] = new int[10];
        for (int j=0; j<y; j++)
            ms[i][j] = obj.ms[i][j];
    }
}

void matr_print (matr &mtr)
{
	for (int i = 0; i < mtr.x; i ++)
	{
		cout << endl;
		for (int j = 0; j < mtr.y; j ++)
			cout << mtr.ms[i][j]<< " ";
	}
	cout << endl;
}

void matr_read (matr &mtr)
{
	cout << "Enter height: " ;
	cin >> mtr.x;
	cout << "Enter width: " ;
	cin >> mtr.y;
	cout << "Enter data :"<< endl;
	for (int i = 0; i < mtr.x; i ++)
		for (int j = 0; j < mtr.y; j ++)
			cin >> mtr.ms[i][j];
}

void main()
{
	matr m1, m2, m3;
	cout << "Enter matrix 1" << endl;
    matr_read(m1);
    cout << "Enter matrix 2" << endl;
	matr_read(m2);
    m3=m1+m2;
    cout << "Matrix 1 + matrix 2 =";
	matr_print(m3);
    cout << "Matrix 1 * matrix 2 =";
	m3=m1&m2;
	matr_print(m3);
    system ("pause");
}