class Matrix

 {

     int row;

     int col;

     int *p;

 public:

     Matrix();

     Matrix(int x,int y);

     ~Matrix();

     Matrix(Matrix &m);

     int getrows();

     int getcols();

     void setelem(int x,int y,int elem);

     int getelem(int x,int y);

     void add(Matrix &m);

     void subtract(Matrix &m);    

     bool equals(Matrix &m);

 };

 #include "Matrix.hh"

 #include <cassert>

 #include <iostream>

 using namespace std;

 Matrix::Matrix()

 {

     row=;

     col=;

     p=;

 }

 Matrix::Matrix(int x,int y)

 {

     assert(x>=);

     assert(y>=);

     row=x;

     col=y;

     p=new int[row*col];

     for(int i=;i<row*col;i++)

         p[i]=;

 }

 Matrix::Matrix(Matrix &m)

 {

     row=m.row;

     col=m.col;

     p=new int[row*col];

     for(int i=;i<row*col;i++)

     p[i]=m.p[i];

 }

 Matrix::~Matrix()

 {

     delete[] p;

 }

 int Matrix::getrows()

 {

     return row;

 }

 int Matrix::getcols()

 {

     return col;

 }

 int Matrix::getelem(int x,int y)

 {

     assert(x>=);

     assert(y>=);

     return p[x*col+y];

 }

 void Matrix::setelem(int x,int y,int elem)

 {

     assert(x>=);

     assert(y>=);

     p[x*col+y]=elem;

 }

 void Matrix::add(Matrix &m)

 {

     assert(    (row==m.row) || (col==m.col) );

     for(int i=;i<row*col;i++)

     p[i]+=m.p[i];

 }

 void Matrix::subtract(Matrix &m)

 {

     assert( (row==m.row) || (col==m.col) );

     for(int i=;i<row*col;i++)

     p[i]-=m.p[i];

 }

 bool Matrix::equals(Matrix &m)

 {

     bool b=true;

     if( (row!=m.row) || (col!=m.col) )return false;

     for(int i=;i<row*col;i++)

         if(p[i]!=m.p[i])

         {    b=false;return b;}

     return b;

 }