对称矩阵的缩放存储

对称矩阵的缩放存储：重点在于需要达到同一个位置被重复访问。
公式记清是最主要的一步，就是因为我公式没记清，浪费了太多时间

：对称矩阵缩放，目的就是为了节省空间，当矩阵很大时，节省的空间就很可观， 编写该程序的最主要问题在于：考虑清楚，当将下三角创建好之后，如何让上三角得到下三角的值？
：当构建上（下）三角后，一维数组就此固定：因为二维数组的值已经全部在一维数组中，将i，j发生相应的替换就可以得到一维数组相应的坐标，从而可以得到 下（上）三角 想得到的值。
这就利用到公式：

/*对称矩阵里的中任一元素在一维数组中的下标k  与i，j的对应关系为:当i>=j,k=i*(i+1)/2+j；当i<j k =j*(j+1)/2+i;*/


#include<iostream>
//对称矩阵的缩放存储
using namespace std;



template<class T>
class symmetric_matrix
{
public:
    symmetric_matrix(T* a,int n)
        :_a(new T [n*(n+1)/2])
        ,_size(n+(n+1)/2)
    {
        for(int i=0;i<n;i++)
        {
            for(int j=0;j<n;j++)
            {
                if(i>=j)                
                    _a[i*(i+1)/2+j]= a[i*n+j];//构建下三角
                /*else
                _a[j*(j+1)/2+i]=a[j*n+i]//构建上三角
                */
            }
        }

    }
    void Print(int n)
    {
        for(int i=0;i<n;i++)
        {
            for(int j=0;j<n;j++)
            {
                if(i>=j)            
                    cout<<_a[i*(i+1)/2+j]<<" ";
                else
                    cout<<_a[j*(j+1)/2+i]<<" ";//看似不同，实际上指向的是同一块空间
            }
            cout<<endl;
        }
        cout<<endl;
    }

    T& Access(int i, int j )
    {
        //if (i < j)//上三角
        //{
        //  swap(i, j);
        //}

        //return _a[i*(i+1)/2+j];
        if(i<j)
            return _a[j*(j+1)/2+i];
        else 
            return _a[i*(i+1)/2+j];
    }

private:
    T* _a;
    int _size;
};
int main()
{
    int array[][3]=
    {
        {1,0,3},    
        {0,1,0},
        {3,0,1}
    };
    symmetric_matrix<int> sm((int*)array,3);
    sm.Print (3);
    cout<<sm.Access (2,2)<<endl;
    return 0;
}