特殊矩阵的压缩存储

---

数组的存储结构

一个数组的所有元素在内存中是占用一段连续的存储空间的。

<1>一维数组

以一维数组$A[0...n-1]$为例，其存储结构关系式为：
$$LOC(a_i)=LOC(a_0)+(i)\times L(0\le i＜n)$$
其中，$L$是每个数组元素所占的存储单元。

例如，对一个一维数组$A[n]$，其存储的首地址为100，每个元素占一个存储单元，求$A[14]$的存储位置？
$A[14]$的存储位置为：$100+14\times 1=114$

<2>多维数组

对于多维数组，有两种映射方法：按行优先和按列优先。

以二维数组为例，按行优先存储的基本思想是：先行后列。
先存储行号较小的元素，行号相等先存储列号较小的元素。
设二维数组行下标与列下标的范围分别为$[l1,h1]$，$[l2,h2]$，则存储结构关系式为：$$LOC(a_{i,j})=LOC(a_{l_1,l_2})+[(i-l_1)\times (h_2-l_2+1)+(j-l_2)]\times L$$
若 $l_1$ 和 $l_2$ 的值均为0 ，即从$A[0][0]$开始，上式变为
$$LOC(a_{i,j})=LOC(a_{0,0})+[i\times (h_2+1)+j)]\times L$$

以二维数组为例，按列优先存储的基本思想是：先列后行。
先存储列号较小的元素，列号相等先存储行号较小的元素。
设二维数组行下标与列下标的范围分别为$[l1,h1]$，$[l2,h2]$，则存储结构关系式为：$$LOC(a_{i,j})=LOC(a_{l_1,l_2})+[(j-l_2)\times (h_1-l_1+1)+(i-l_1)]\times L$$
若 $l_1$ 和 $l_2$ 的值均为0 ，即从$A[0][0]$开始，上式变为
$$LOC(a_{i,j})=LOC(a_{0,0})+[j\times (h_1+1)+i)]\times L$$

---

特殊矩阵的压缩存储

压缩存储：指为多个值相同的元素只分配一个存储空间，对零元素不分配存储空间。其目的是为了节省存储空间。

特殊矩阵：指具有许多相同矩阵元素或零元素，并且这些相同矩阵元素或零元素的分布有一定规律性的矩阵。常见的特殊矩阵有对称矩阵、上（下）三角矩阵、对角矩阵等。

---

对称矩阵

若对一个 $n$ 阶方阵$A[1...n][1...n]$中的任意一个元素 $a_{i,j}$都有 $a_{i,j}$ = $a_{j,i}(1\le i,j\le n)$，则称其为对称矩阵。

对于$n$阶方阵，其中的元素可以划分为3个部分，即上三角区、主对角线和下三角区。

对于$n$阶对称矩阵，上三角区的所有元素和上三角区的对应元素相同，若仍采用二维数组存放，则会浪费几乎一半的空间，为此，将对称矩阵$A[1...n][1...n]$存放在一维数组$B[n(n+1)/2]$中，即元素$a_{i,j}$存放在$b_k$中。只存放主对角线和下三角区（或上三角区）的元素，采用行优先存储。

那么在数组B中，位于元素$a_{i,j}(i\ge j)$前面的元素个数为：
第1行：($a_{1,1}$) 1个元素
第2行：($a_{2,1}$，$a_{2,2}$) 2个元素
第3行：($a_{3,1}$，$a_{3,2}$，$a_{3,3}$) 3个元素
……
第$i$-1行：($a_{i-1,1}$，$a_{i-1,2}$，…，$a_{i-1,i-1}$) $i$-1个元素
第$i$行：($a_{i,1}$，$a_{i,2}$，…，$a_{i,j-1}$) $j$-1个元素

那么，元素 $a_{i,j}$ 在数组$B$中的下标 $k=1+2+...+(i-1)+(j-1)=i(i-1)/2+j-1$（数组下标从0开始）

因此，元素下标之间的对应关系如下：
$$k=\{\begin{array}{l}\frac{i(i-1)}{2}+j-1,i\ge j(下三角区和主对角线元素)\\ \\ \frac{j(j-1)}{2}+i-1,i＜j(上三角区元素a_{i,j}=a_{j,i})\end{array}$$

若对称矩阵的压缩存储存放上三角区和主对角线的元素（按行优先），那么

那么在数组B中，位于元素$a_{i,j}(i\le j)$前面的元素个数为：
第1行：($a_{1,1}$，$a_{1,2}$，…，$a_{1,n}$) $n$个元素
第2行：($a_{2,2}$，$a_{2,3}$，…，$a_{2,n}$) $n$-1个元素
第3行：($a_{3,3}$，$a_{3,4}$，…，$a_{3,n}$) $n$-2个元素
……
第$i$-1行：($a_{i-1,i-1}$，$a_{i-1,i}$，…，$a_{i-1,n}$) $n$-$i$+2个元素
第$i$行：($a_{i,i}$，$a_{i,i+1}$，…，$a_{i,j-1}$) $j$-$i$个元素

那么，元素 $a_{i,j}$ 在数组$B$中的下标 $k=n+(n-1)+...+(n-i+2)+(j-i)=(i-1)(2n-i+2)/2+j-i$（数组下标从0开始）

因此，元素下标之间的对应关系如下：
$$k=\{\begin{array}{l}\frac{(i-1)(2n-i+2)}{2}+j-i,i\ge j(上三角区和主对角线元素)\\ \\ \frac{(j-1)(2n-j+2)}{2}+i-j,i＜j(下三角区元素a_{i,j}=a_{j,i})\end{array}$$

---

三角矩阵

下三角矩阵中，上三角区的所有元素均为同一常量。
其存储思想与对称矩阵一致，不同之处在于存储完下三角区和主对角线上的元素之后，紧接着存储对角线上方的常量一次，故将下三角矩阵$A[1...n][1...n]$存放在一维数组$B[n(n+1)/2+1]$中（按行优先）。

元素下标之间的对应关系为：
$$k=\{\begin{array}{l}\frac{i(i-1)}{2}+j-1,i\ge j(下三角区和主对角线元素)\\ \\ \frac{n(n+1)}{2},i＜j(上三角区元素)\end{array}$$

上三角矩阵中，下三角区的所有元素均为同一常量。
与下三角矩阵的存储一致，存储完上三角区和主对角线上的元素之后，紧接着存储下三角区的常量一次，故将下三角矩阵$A[1...n][1...n]$存放在一维数组$B[n(n+1)/2+1]$中（按行优先）。

元素下标之间的对应关系为：
$$k=\{\begin{array}{l}\frac{(i-1)(2n-i+2)}{2}+j-i,i\ge j(上三角区和主对角线元素)\\ \\ \frac{n(n+1)}{2},i＜j(下三角区元素)\end{array}$$

---

三对角矩阵

对角矩阵也称为带状矩阵。对角矩阵的概念：
设有一个$n\times n$的矩阵$A$，对于矩阵$A$中的任意元素$a_{i,j}$，当 $|i-j|>1$时有$a_{i,j}=0$，这样的矩阵称为三对角矩阵。
$$\left[\begin{array}{ccccccc}{a}_{1,1}& a_{1,2}& & & & & \\ a_{2,1}& a_{2,2}& a_{2,3}& & & & \\ & a_{3,2}& a_{3,3}& a_{3,4}& & & \\ & & & & & & \\ & & \ddots & \ddots & \ddots & & \\ & & & a_{n-1,n-2}& a_{n-1,n-1}& a_{n-1,n}& \\ & & & & a_{n,n-1}& a_{n,n}& \end{array}\right]$$

三对角矩阵的压缩存储：
在三对角线矩阵中，除了主对角线及在主对角线上下最邻近的两条对角线上的元素外，所有其他元素均为0。为了节省存储空间，只存储主对角线及其上、下两侧次对角线上的元素外，其他零元素一律不存储。

同样需要用一个一维数组$B$来存储三对角矩阵中位于三对角线上的元素。同样要区分两种存储方式：即行优先方式和列优先方式。

若将三对角线矩阵$A$中的元素按行优先方式存放在一维数组$B$中，则

那么在数组B中，位于元素$a_{i,j}(i\le j)$前面的元素个数为：
第1行：($a_{1,1}$，$a_{1,2}$) 2个元素
第2行：($a_{2,1}$，$a_{2,2}$，$a_{2,3}$) 3个元素
……
第$i$-1行：($a_{i-1,i-2}$，$a_{i-1,i-1}$，$a_{i-1,i+1}$) 3个元素
第$i$行：($a_{i,i-1}$，…，$a_{i,j-1}$) $j$-$i$+1个元素

那么，元素 $a_{i,j}$ 在数组$B$中的下标 $k=2+3(i-2)+j-i+1=2i+j-3$（数组下标从0开始），那么就有$$\{\begin{array}{l}i=\lfloor (k+1)/3+1\rfloor \\ \\ j=k-2i+3\end{array}$$

---

稀疏矩阵

若一个矩阵中的非零元素个数远远小于矩阵元素总数，则称该矩阵为稀疏矩阵。

稀疏矩阵的存储：

若用常规二维数组的方法来存储系数矩阵，则相当浪费存储空间，因此对稀疏矩阵仅存储非零元素。

但通常非零元素的分布没有规律，故不仅要存储非零元素的值，还要存储它所在的行和列。因此，可用一个三元组（行标，列标，值）来存储稀疏矩阵。

此外，稀疏矩阵也可用链接存储。采用十字链表表示，它为每个非零元素结点附带了两个链接指针：

• 一个指示同一行下一个非零元素结点（行指针）
• 另一个指示同一列下一个非零元素结点（列指针）

在狮子链表中，每一个三元组结点按照矩阵元素所在行号$i$和列号$j$链接入第$i$个行链表和第$j$个列链表中，即处于所在的行链表和列链表的交汇处。