对称矩阵的存储

对称矩阵是以对角线为对称轴，每个元素相对对角线是相等的，对于存储这样的数组，如果对每个元素都进行存储将会造成空间的浪费，因此，一般选择只存储下三角这些元素，存储方法代码如下：
template<class T>
class SymmetricMatrix
{
private:
T *_a;
int _n;
int _size;
public:
SymmetricMatrix(T *a, int size)
:_a(new T[size*(size + 1) / 2])
, _size(size)
, _n(size*(size + 1) / 2)
{
int i = 0;
int j = 0;
int index = 0;
for (i = 0; i < size; i++)
{
for (j = 0; i < size; j++)
{
if (i >= j)
{
_a[index] = a[i*size + j];
}
else
{
break;
}
index++;
}
}
}
~SymmetricMatrix()
{
if (_a != NULL)
delete[]_a;
}
void Display()
{
int i = 0;
int j = 0;
for (i = 0; i <_size; i++)
{
for (j = 0; j < _size; j++)
{
if (i >=j)
{
cout << _a[i*(i + 1) / 2 + j] << " ";
}
else
{
cout << _a[j*(j + 1) / 2 + i] << " ";
}
}
cout << "\n" << endl;
}
}
T & Access(size_t i, size_t j)
{
if (i < j)
{
swap(i, j);
}
return _a[i*(i + 1) / 2 + j];
}
};
void test()
{
int a[3][3] = { { 1, 2, 3 }, { 2, 1, 4 }, { 3, 4, 1 } };
SymmetricMatrix<int> s((int*)a, 3);
s.Display();
}
int main()
{
test();
getchar();
return 0;
}