/**

 * @brief C语言 对称矩阵 压缩 实现

 * @author wid

 * @date 2013-11-03

 *

 * @note 若代码存在 bug 或程序缺陷, 请留言反馈, 谢谢!

 */

 #include <stdio.h>

 #include <stdlib.h>

 #include <assert.h>

 #define MAT_ROW 6

 #define MAT_COL 6

 typedef int ElemType;

 ///方法声明

 ElemType *Compress( ElemType pArr2D[][MAT_COL] );       ///压缩二维对称矩阵

 ElemType Value( ElemType *pComArr, int m, int n );      ///从压缩后的矩阵中取值

 ///方法实现

 /**

 * @brief 压缩对称矩阵

 *

 * @param 待压缩的矩阵

 *

 * @return 返回压缩后线性表的指针

 */

 ElemType *Compress( ElemType pArr2D[][MAT_COL] )

 {

     ElemType *pMat = (ElemType *)malloc( sizeof(ElemType) * ((MAT_COL * (MAT_COL + )) / ) );

     ///映射到线性表中

     int m = , n = , nLen = , nPos = ;

     for( m = ; m < MAT_ROW; ++m )

     {

         for( n = ; n < nLen; ++n )

         {

             pMat[nPos] = pArr2D[m][n];

             ++nPos;

         }

         ++nLen;

     }

     return pMat;

 }

 /**

 * @brief 从压缩后得到的线性表中取值

 *

 * @param pComArr 指向存放压缩矩阵的线性表

 * @param m 矩阵的第一维下标

 * @param n 矩阵的第二维下标

 *

 * @return 返回该下标指向的元素值

 *

 * @note 元素位置由 0 计起

 */

 ElemType Value( ElemType *pComArr, int m, int n )

 {

     ///使用断言

     assert( m >=  && m < MAT_ROW &&  n >=  && n < MAT_COL );

     ///将0索引转化为1索引

     ++m;

     ++n;

     ///根据对称矩阵性质从压缩一维数组中取值

     if( m >= n )

         return pComArr[ (m * (m - )) /  + n -  ];

     else

         return pComArr[ (n * (n - )) /  + m -  ];

 }

 ///测试

 int main()

 {

     ///对称矩阵 arrMat

     ElemType arrMat[MAT_ROW][MAT_COL] = {

         {, , , , , },

         {, , , , , },

         {, , , , , },

         {, , , , , },

         {, , , , , },

         {, , , , , }

     };

     ///对矩阵 arrMat 进行压缩

     ElemType *pCom = Compress( arrMat );

     ///从压缩矩阵输出整个矩阵, 使用 0 索引

     int m = , n = ;

     for( m = ; m < MAT_ROW; ++m )

     {

         for( n = ; n < MAT_COL; ++n )

         {

             printf( "%d ", Value(pCom, m, n) );

         }

         putchar( '\n' );

     }

     return ;

 }