對稱矩陣及對稱矩陣的壓縮存儲

設一個N*N的方陣A，A中任意元素Aij，當且僅當Aij == Aji(0 <= i <= N-1 && 0 <= j <= N-1)，則矩陣A是對稱矩陣。以矩陣的對角線為分隔，分為上三角和下三角。

壓縮存儲稱矩陣存儲時只需要存儲上三角/下三角的數據，所以最多存儲n(n+1)/2個數據。

對稱矩陣和壓縮存儲的對應關系：下三角存儲i>=j,  SymmetricMatrix[i][j] == Array[i*(i+1)/2+j]

        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

Symmetry.h中
template < class T>
class Symmetry
{
public:
	//構造函數
	Symmetry(T* arr, size_t size)
		:_arr(new T[size*(size+1)/2])
		, _size(size*(size + 1)/2)
	{
		for (int i = 0; i < size; i++)
		{
			for (int j = 0; j < size; j++)
			{
				if (i >= j)
				{
					_arr[i*(i + 1) / 2 + j] = arr[i*size + j];//把對稱矩陣壓縮
				}
			}
		}
	}
	//打印
	void Print(size_t size)
	{
		for (int i = 0; i < size; i++)
		{
			for (int j = 0; j < size; j++)
			{
				int row = i;
				int col = j;
				if (row < col)
				{
					swap(row, col);
				}
					cout << _arr[row*(row+ 1) / 2 + col] << " ";
			}
			cout << endl;
		}
		cout << endl;
	}
protected:
	T *_arr;
	size_t _size;
};
test.cpp中
#include <iostream>
using namespace std;
#include "Symmetry.h"
void Test()
{
	int arr[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}
	};
	Symmetry<int>s((int*)arr, 5);
	s.Print(5);
}
int main()
{
	Test();
	system("pause");
	return 0;
}