#include<iostream>
using namespace std;

typedef float elemType;//方便我們轉換為template形式

class Matrix
{
	//friend宣告不受訪問許可權的影響，放在class開始處
	friend Matrix operator+(const Matrix&, const Matrix&);  //矩陣加法
	friend Matrix operator*(const Matrix&, const Matrix&);  //矩陣乘法

public:
	//宣告建構函式
	Matrix(
 
);
	Matrix(const elemType*);
	Matrix(elemType a11, elemType a12, elemType a13, elemType a14,
		elemType a21, elemType a22, elemType a23, elemType a24,
		elemType a31, elemType a32, elemType a33, elemType a34,
		elemType a41, elemType a42, elemType a43, elemType a44);
	

	int rows() const { return 4;
 
 }
	int cols() const { return 4; }

	ostream& print(ostream& os) const;
	void operator+=(const Matrix &rhs);

	elemType operator()(int row, int column) const //過載函式呼叫運算子()
	{
		return _matrix[row][column];
	}

	elemType &operator()(int row, int column) 
	{
		return _matrix[row][column]
 
;
	}

private:
	elemType _matrix[4][4];

};

//定義建構函式
//預設建構函式
Matrix::Matrix()
{
	_matrix[0][0] = 0.; _matrix[0][1] = 0.; _matrix[0][2] = 0.; _matrix[0][3] = 0.;
	_matrix[1][0] = 0.; _matrix[1][1] = 0.; _matrix[1][2] = 0.; _matrix[1][3] = 0.;
	_matrix[2][0] = 0.; _matrix[2][1] = 0.; _matrix[2][2] = 0.; _matrix[2][3] = 0.;
	_matrix[3][0] = 0.; _matrix[3][1] = 0.; _matrix[3][2] = 0.; _matrix[3][3] = 0.;
}
//有參無預設建構函式，接受一個數組
Matrix::Matrix(const elemType *array)
{
	int array_index = 0;
	for (int ix = 0; ix < 4; ++ix)
	    for (int jx = 0; jx < 4; ++jx)
		    _matrix[ix][jx] = array[array_index++];
}
//有參無預設建構函式，接受16個數據值
Matrix::Matrix(elemType a11, elemType a12, elemType a13, elemType a14,
	elemType a21, elemType a22, elemType a23, elemType a24,
	elemType a31, elemType a32, elemType a33, elemType a34,
	elemType a41, elemType a42, elemType a43, elemType a44)
{
	_matrix[0][0] = a11; _matrix[0][1] =a12; _matrix[0][2] = a13; _matrix[0][3] =a14;
	_matrix[1][0] = a21; _matrix[1][1] = a22; _matrix[1][2] = a23; _matrix[1][3] =a24;
	_matrix[2][0] = a31; _matrix[2][1] = a32; _matrix[2][2] = a33; _matrix[2][3] =a34;
	_matrix[3][0] = a41; _matrix[3][1] = a42; _matrix[3][2] = a43; _matrix[3][3] = a44;
}

inline ostream & operator<<(ostream& os, const Matrix &rhs)
{
	return rhs.print(os);
}
//定義print函式，用於輸出流過載
ostream& Matrix::print(ostream& os) const
{
	int cnt = 0;
	for (int ix = 0; ix < 4; ++ix)
	{
		for (int jx = 0; jx < 4; ++jx)
		{
			if (cnt && !(cnt % 4)) os << endl;
			os << _matrix[ix][jx] << ' ';
			++cnt;
		}
	}
	os << endl;
	return os;
}

//定義+ 過載運算子
Matrix operator+( const Matrix &m1, const Matrix &m2)
{
	Matrix result;
	for (int ix = 0; ix < m1.rows(); ++ix)
	{
		for (int jx = 0; jx < m1.cols(); ++jx)
		{
			result(ix, jx) = 0;
			result(ix, jx) = m1(ix, jx) + m2(ix, jx);
		}
	}
	return result;
}

//定義 * 過載運算子
Matrix operator*(const Matrix &m1, const Matrix &m2)
{
	Matrix result;
	for (int ix = 0; ix < m1.rows(); ++ix)
	{
		for (int jx = 0; jx < m1.cols(); ++jx)
		{
			result(ix, jx) = 0;
			for (int kx = 0; kx < m1.cols();++kx)
			    result(ix, jx) += m1(ix, kx) * m2(kx, jx);
		}
	}
	return result;
}

//定義 += 過載運算子
void Matrix::operator+=(const Matrix &m)
{
	for (int ix = 0; ix < m.rows();++ix)
	    for (int jx = 0; jx < m.cols(); ++jx)
		    _matrix[ix][jx] += m(ix, jx);
}

int main()
{
	Matrix m;
	cout << m << endl;

	elemType ar[16] = { 1, 0, 0, 0, 0, 1, 0, 0,
		0, 0, 1, 0, 0, 0, 0, 1 };
	Matrix m1(ar);
	cout << m1 << endl;

	Matrix m2(m1);
	cout << m2 << endl;

	elemType ar2[16] = {
		1.3, 0.4, 2.6, 8.2, 6.2, 1.7, 1.3, 8.3,
		4.2, 7.4, 2.7, 1.9, 6.3, 8.1, 5.6, 6.6
	};
	Matrix m3(ar2);
	cout << m3 << endl;

	Matrix m4 = m3*m1;
	cout << m4 << endl;

	Matrix m5 = m3 + m4;
	cout << m5 << endl;

	m3 += m4;
	cout << m3 << endl;

	system("pause");
}