對稱矩陣的壓縮儲存講解

2020-01-26 13:31:23

字體：大中小

來源：轉(zhuǎn)載

供稿：網(wǎng)友

一、存儲矩陣用一個二維數(shù)組即可；

二、什么是對稱矩陣：

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

三、對稱矩陣的壓縮儲存：

壓縮存儲稱矩陣存儲時只需要存儲上三角/下三角的數(shù)據(jù)，所以最多存儲n(n+1)/2個數(shù)據(jù)（相當于1+2+…+n,即等差數(shù)列求和）。

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

四、代碼實現(xiàn)

#include<iostream>using namespace std;template<class T>class CompressionMatrix{public:  CompressionMatrix(T* arr,int sz)    :_data(new T[sz*(sz+1)/2])    ,_size(sz)  {    int index=0;    //壓縮儲存過程    for(int i=0;i<sz;++i)    {      for(int j=0;j<sz;++j)      {        if (i>=j)//_data中儲存下三角的數(shù)據(jù)        {          _data[index]=arr[i*sz+j];          index++;        }        else          break;      }    }  }  //獲取某個坐標的數(shù)據(jù)，i和j代表該數(shù)據(jù)在矩陣中的橫縱坐標  T GetDate(int i,int j)  {    if (i>=j)//下三角數(shù)據(jù)    {      return _data[i*(i+1)/2+j];    }    else//上三角數(shù)據(jù)    {      std::swap(i,j);//將橫坐標和從坐標值交換；      return _data[i*(i+1)/2+j];    }  }    //打印矩陣的數(shù)據(jù)  void PrintfMatrix()  {    for (int i=0;i<_size;++i)    {      for (int j=0;j<_size;++j)      {        cout<<GetDate(i,j)<<" ";      }      cout<<endl;    }  }  ~CompressionMatrix()  {    if (_data!=NULL)    {      delete[] _data;      _data=NULL;      _size=0;    }  }protected:  T* _data;//儲存數(shù)據(jù)的數(shù)組  int _size;//儲存原始對稱矩陣的行數(shù)（或列數(shù)）};

測試代碼：

int main(){  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},  };  CompressionMatrix<int> cm((int*)a,5);//將二維數(shù)組強制轉(zhuǎn)換為一維數(shù)組指針，是問題更簡單  cm.PrintfMatrix();  return 0;}

五、運行結(jié)果