对称矩阵及对称矩阵的压缩存储

1.对称矩阵：设一个N*N的方阵A，A中任意元素Aij，当且仅当Aij == Aji(0 <= i <= N-1 && 0 <= j <= N-1)，则矩阵A是对称矩阵。以矩阵的对角线为分隔，分为上三角和下三角。

2.压缩存储称矩阵存储时只需要存储上三角/下三角的数据，所以最多存储n(n+1)/2个数据

（根据等差数列的前n项和可求得）。

3.对称矩阵和压缩存储的对应关系：下三角存储i>=j, SymmetricMatrix[i][j] == Array[i*(i+1)/2+j]。

4.这就是一个对称矩阵。

template<class T>
class SymmetricMatrix
{
public:
	SymmetricMatrix(T* a,size_t n)    //i>=j
		:_n(n)
		
	{
		_a=new T[n*(n+1)/2];
		size_t k=0;
		for(size_t i=0;i<n;i++)
		{
			for(size_t j=0;j<n;j++)
			{
				if(i>=j)
				{
					  _a[k++]=a[i*_n+j];
				}
			}
		}
	}
	~SymmetricMatrix()
	{
		if(_a!=NULL)
		{
			delete[] _a;
			_n=0;
		}
	}
	T& Access(size_t i,size_t j)
	{
		if(i<j)
		{
			swap(i,j);
		}
		return _a[i*(i+1)/2+j];
	}
	void Display()
	{
		for(size_t i=0;i<_n;i++)
		{
			for(size_t j=0;j<_n;j++)
			{
				cout<<Access(i,j)<<" ";
			}
			cout<<endl;
		}
	}

private:
	T* _a;   //数组
	size_t _n;   //行和列数
};
void TestSymmetricMatrix()
{
	int a[5][5]=
	{
		{0,1,2,3,4},
		{1,0,1,2,3},
		{2,1,0,1,2},
		{3,2,1,0,1},
		{4,3,2,1,0},
	};
	SymmetricMatrix<int> sm((int*)a,5);
    //cout<<sm.Access (4,0)<<endl;
	sm.Display();
}