对称矩阵的存储及转置算法

本文要点：

1.对称矩阵与稀疏矩阵

2.两种矩阵的压缩存储

3.代码实现两种矩阵

对称矩阵<SymmetryMatrix>

1.对称矩阵也是一种特殊矩阵，满足Aij = Aji(设矩阵为A,且有0<=i<N-1 && 0<=j<N-1),这种矩阵以对角线分割为上三角和下三角，关于对角线对称的元素相等。

2.对称矩阵的压缩存储

如果把矩阵中的每个元素都存储起来，那么就会显得浪费空间，因为每两个关于对角线对称的元素相等，因此就可以将矩阵压缩存储到一个数组Array中，即：将对称的两个元素存一份，对角线上的元素都存储起来，也就是说只存储上三角或下三角，那么存储的元素的总个数就为n(n+1)/2个，这样，当n特别大的时候，也就最有效。

压缩存储的数组与矩阵之间满足：Array[i*(i+1)/2+j] = Martix[i][j](下三角存储i>=j)

3.代码实现

template<class T>
class SymmetryMatrix
{
public:
	SymmetryMatrix(const T* a,size_t n)
		:_matrix(new T[n*(n+1)/2])
		,_size(n*(n+1)/2)
		,_n(n)
	{
		size_t index = 0;	//表示一维数组的下标
		for (size_t i=0; i<n; ++i)
		{
			for (size_t j=0; j<n; ++j)
			{
				if (i >= j)		//存储下三角
				{
					_matrix[index++] = a[i*n+j];
					//index++;
				}
				else
					continue;
			}
		}
	}

	T& Access(size_t row,size_t col)	//访问数据
	{
		if (row < col)
		{
			std::swap(row,col);
		}
		return _matrix[row*(row+1)/2+col];
	}

	void Display()
	{
		for (size_t i=0; i<_n; i++)
		{
			for (size_t j=0; j<_n; j++)
			{
				cout<<Access(i,j)<<" ";
			}
			cout<<endl;
		}
	}
protected:
	T* _matrix;		//压缩存储的一维数组
	size_t _size;	//可存储的大小
	size_t _n;		//行列的大小		
};

稀疏矩阵<SparseMatrix>

1.稀疏矩阵中，有效数据（非0）的数量远小于非法数据的个数（0为非法数据）