对称矩阵的压缩存储 c++模板类实现

首先我们有一个矩阵，如：

{ 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 },

这种行列数相同，且以对角线为中心对称的矩阵称为对称矩阵

对称矩阵因其特点，可以进行压缩存储（只存上半部或下半部）以达到节省空间的目的

需要矩阵时可以将其还原

下面给出对称矩阵压缩存储的一种思路(存储下半部分),使用c++模板实现:

.h:

template < class T, size_t N >
class SymmetricMatrix
{
protected:
		T* _matrix;								//用一维数组存储二维矩阵(数组)下半部分
		size_t _n;
public:
		SymmetricMatrix( T* a, size_t n )		//二维数组实现机制，底层以一维数组实现 所以 参数为 T*
			:_n ( n )
		{
			_matrix = new T[ n * ( n + 1 ) / 2 ];//注意T后圆括号与方括号的区别//只需存储 n * (n + 1) / 2个元素（根据等差数列推导）

			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
					{
						break;					//跳过此行，减少程序执行步骤
					}
				}
			}
		}

		T& Access ( size_t i, size_t j )
		{
			if ( i < j )
			{
				swap ( i, j );
			}

			return _matrix[ i * ( i + 1 ) / 2 + j ];//也是根据等差数列，得出i行前有多少个元素 + j即为此时元素
		}
		
		const T & Access ( size_t i, size_t j ) const
		{
			if ( i < j )
			{
				swap ( i, j );
			}

			return _matrix[ i * ( i + 1 ) / 2 + j ];// i * (i + 1) / 2为前i行元素
		}

		void Display (  )
		{
			for ( size_t i = 0; i < _n; ++i )
			{
				for ( size_t j = 0; j < _n; ++j )
				{
					cout << Access ( i , j ) << " ";
				}

				cout << endl;
			}

			cout << endl;
		}

		~SymmetricMatrix (  )
		{
			delete [] _matrix;
			_matrix = NULL;

			_n = 0;
		}
};

.cpp:

#include <iostream>
#include <windows.h>
using namespace std;
#include "yy.h"

void test ( )
{
	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,5> sm( ( int * )a, 5 );	//构建对象

	sm.Display ( );
}

int main ( )
{
	test ( );

	system ( "pause" );

	return 0;
}