第6章串数组：特殊矩阵的压缩存储 - 详解

6.6 特殊矩阵的压缩存储

在计算机科学中，用二维数组来表示数学中的矩阵，是最自然的方法。

所谓的压缩存储，从而可以提高存储空间利用率。就是有一类矩阵，其非零元素或者零元素的分布有一定规律，比如：对称矩阵、三角矩阵、对角矩阵等。为了节省存储空间，可以利用这类矩阵中元素的分布规律，为多个值相等的元素只分配一个存储空间，对零不分配空间。这就

1. 对称矩阵

对于 $n$ 阶方阵 $A=(a_{ij}{)}_{n\times n}$，位置索引值$i=j$的那些元素，构成了矩阵的主对角线（Main Diagonal）。

• $i=j$ ，$a_{ij}$构成了主对角线元素；
• $i<j$ ，$a_{ij}$构成了上三角部分元素；
• $i>j$ ，$a_{ij}$构成了下三角部分元素；

$$\left[\begin{array}{cccc}{a}_{11}& a_{12}& \cdots & a_{1n}\\ a_{21}& a_{22}& \cdots & a_{2n}\\ ⋮& ⋮& \ddots & ⋮\\ a_{n1}& a_{n2}& \cdots & a_{nn}\end{array}\right]$$

以主对角线为对称轴，两侧元素对称分布的对角矩阵，即$a_{ij}=a_{ji},i\ne j$，这样的矩阵称为对称矩阵（Symmetric Matrix），即上三角部分和下三角部分中的对应元素相等。

如果将对称矩阵存储到一维的存储结构中，为了讲解和操作方便，用一维数组表示一维的存储结构，在 6.5 节中已经讲解过，假设将一维数组保存到内存中，该数组的元素和内存中的存储单元一一对应。所以，能够假设一维数组$B$ 作为 $n$ 阶对称矩阵 $A$的存储结构，$A$ 中的元素 $a_{ij}$ 存储在 $B$ 中的 $b_k$ 处。

在存储对称矩阵时，可以只存储主对角线和上三角部分的元素，或者是主对角线和下三角部分的元素，让对称的两个元素共享一个存储空间。

不失一般性，根据对称矩阵的特点，采用以行序为主序，存储主对角线和下三角部分的元素，如图 6.6.1 所示。

图 6.6.1 压缩对称矩阵

下面探讨 $k$ 与 $i,j$ 的关系（$k$是从 0 开始；$i,j$都是从 1 开始）。

（1）若 $a_{ij}$ 是 $A$中主对角线或下三角部分的元素，则有$i\ge j$（以行序为主序）。$a_{ij}$ 处于第 $i$行，那么在第$1$ 行到第 $i-1$行中，必须保存如下元素：

• 第 1 行有 1 个元素：$a_{11}$

• 第 2 行有 2 个元素：$a_{21},a_{22}$

• $\cdots$

• 第 $i-1$ 行有 $i-1$ 个元素：$a_{(i-1)1},a_{(i-1)2},\cdots ,a_{(i-1)(i-1)}$

以上共计元素数量：$1+2+\cdots +(i-1)$

• 在第 $i$ 行的元素 $a_{ij}(j\le i)$：就是的前面的元素$a_{i1},a_{i2},\cdots ,a_{i(j-1)}$ ，共计有 $j-1$ 个元素

综上，在 $a_{ij}$：就是之前的元素总数量
$$[1+2+\cdots +(i-1)]+(j-1)=\frac{i(i-1)}{2}+(j-1)$$
由于 $k$是从 0 开始，$a_{ij}$ 存储到 $b_k$ 处，在 $b_k$ 之前，数组 $B$ 中共有 $k$个元素，对应矩阵$A$ 中 $a_{ij}$之前的元素，所以则有：
$$k=\frac{i(i-1)}{2}+(j-1),(i\ge j)\text{\hspace{1em}(1)}$$
（2）若 $a_{ij}$ 是 $A$中的上三角部分的元素，即$i<j$ ，其值与 $a_{ji}$相等。根据（1）式，元素$a_{ji}$ 在数组 $B$中对应的存储位置$b_k$ ，则此时 $k$ 与 $i,j$ 的关系是：
$$k=\frac{j(j-1)}{2}+(i-1),(i<j)\text{\hspace{1em}(2)}$$
将（1）式和（2）式综合起来，得到$k$ 与 $i,j$的关系如下：
$$k=\{\begin{array}{l}\frac{i(i-1)}{2}+(j-1),(i\ge j)\\ \frac{j(j-1)}{2}+(i-1),(i<j)\end{array}\text{\hspace{1em}(3)}$$
作为存储对称矩阵的一维数组$B$，其空间大小应该是能够将主对角线和下三角部分的所有元素都存储于其中，即从第 1 行到第$n$行，总计存储的元素数量是：
$$1+2+\cdots +n=\frac{n(n+1)}{2}$$
数组 $B$的分量从 0 开始编号，即$b_0$是第一个元素，则第$n$个元素的编号是$\frac{n(n+1)}{2}-1$ ，即 $b_{n(n+1)/2-1}$是最后一个数组元素。

**重要方法归纳：**在计算矩阵$A$ 中元素 $a_{ij}$ 在数组 $B$ 中存储位置 $k$时，计算步骤如下：

1. 求出元素 $a_{ij}$前面共存放多少个元素，假设为$m$ 个；
2. 注意数组 $B$存放元素的下标开始值，设下标开始值是$s$（一般是从 0 开始，有的从 1 开始）。

则 $k=m+s$ 。

2. 三角矩阵

如果 $n$阶矩阵的下三角部分的元素均为常数$c$或零，那么该矩阵称为上三角矩阵（Upper Triangular Matrix）。同样，矩阵的上三角部分的元素均为常数$c$或零，该矩阵称为下三角矩阵（Lower Triangular Matrix）。

对三角矩阵进行压缩存储，除了存储上三角或下三角部分以及主对角线的元素之外，另外还要存储常数$c$。于是，可能采用存储对称矩阵的思想，但是要对数组$B$在末端增加一个位置，用以存储常数$c$ 。

（1）下三角矩阵

假如将下三角部分和主对角线的元素，存储到一维数组$B$中，由对称矩阵的压缩存储可知，数组$B$中的元素个数是$\frac{n(n+1)}{2}$ ，如果 $B$的下标是从 0 开始的，则最后一个元素的下标是$\frac{n(n+1)}{2}-1$；再加上另一部分的常数$c$，则可知，若将三角矩阵压缩到一维数组$B$ 中，需要 $\frac{n(n+1)}{2}+1$个位置，即数组$B$：就是的下标范围$0\sim \frac{n(n+1)}{2}$（含两端的值），如图 6.6.2 所示。

图 6.6.2 下三角矩阵

下面分析，以行序为主序，将$A$ 中的元素 $a_{ij}$ 存储到 $B$ 的元素 $b_k$ 中，$k$ 和 $i,j$的关系，其中，$k$是从 0 开始，$i,j$是从 1 开始。

• 若 $a_{ij}$ 是 $A$中主对角线或下三角部分的元素，则有$i\ge j$。以行序为主序。

  $a_{ij}$ 处于第 $i$行，那么在第$1$ 行到第 $i-1$行中，需要保存如下元素：

  • 第 1 行有 1 个元素：$a_{11}$

  • 第 2 行有 2 个元素：$a_{21},a_{22}$

  • $\cdots$

  • 第 $i-1$ 行有 $i-1$ 个元素：$a_{(i-1)1},a_{(i-1)2},\cdots ,a_{(i-1)(i-1)}$

  • 在第 $i$ 行的元素 $a_{ij}(j\le i)$的前面的元素是：$a_{i1},a_{i2},\cdots ,a_{i(j-1)}$ ，共计有 $j-1$ 个元素

  综上，在 $a_{ij}$之前的元素总数量：
  $$[1+2+\cdots +(i-1)]+(j-1)=\frac{i(i-1)}{2}+(j-1)$$
  由于 $k$是从 0 开始，$a_{ij}$ 存储到 $b_k$ 处，在 $b_k$ 之前，数组 $B$ 中共有 $k$个元素，对应矩阵$A$ 中 $a_{ij}$之前的元素，所以则有：
  $$k=\frac{i(i-1)}{2}+(j-1),(i\ge j)$$

• 上三角部分的元素为常数$c$ ，即 $i<j$ 时。用 $B$中最后一个位置保存，即下标为$n(n+1)/2$ 的元素。

将以上两种情况合起来，得到$k$ 与 $i,j$ 的关系：
$$k=\{\begin{array}{llc}\frac{i(i-1)}{2}+(j-1),& & (i\ge j)\\ \frac{n(n+1)}{2},& & (i<j)\end{array}\text{\hspace{1em}(4)}$$
（2）上三角矩阵

对于上三角矩阵，也用类似的推导方法，得到$k$ 与 $i,j$ 的关系。

总结：当三角矩阵$A$中的元素下标编号是从 1 开始的时候，即$i,j$从 1 开始计数，但数组$B$ 的 $k$从 0 开始（数组的一般规定），那么数组就是仍然$B$ 的下标 $k$与三角矩阵的下标$i,j$的对应关系是：
$$\begin{array}{rl}& \text{上三角矩阵}\\ & k=\{\begin{array}{llc}\frac{(i-1)(2n-i+2)}{2}+(j-i),& & (i\le j)\\ \frac{n(n+1)}{2},& & (i>j)\end{array}\\ & \text{下三角矩阵}\\ & k=\{\begin{array}{llc}\frac{i(i-1)}{2}+j-1,& & (i\ge j)\\ \frac{n(n+1)}{2},& & (i<j)\end{array}\end{array}\text{\hspace{1em}(5)}$$
3. 对角矩阵

若一个 $n$ 阶方阵 $A$满足其所有非零元素都集中在以主对角线为中心的带状区域中，则称此矩阵为$n$阶对角矩阵（Diagonal Matrix）。

对角矩阵主对角线上、下方各有$b$条非零元素构成的次对角线，称$b$矩阵半带宽，$(2b+1)$为矩阵的带宽。

对于半带宽为$b(0\le b\le \frac{n-1}{2})$的对角矩阵，其$|i-j|\le b$ 的元素 $a_{ij}$不为零，其余元素为零，即带宽之外的元素皆为零（如图 6.6.3 所示）。

图 6.6.3 对角矩阵

如果 $b=1$，则得到如下所示的对角矩阵（为了显现清楚，这里以$5\times 5$矩阵为例），该矩阵有一条主对角线和两条次对角线，故称为三对角矩阵。
$$A=\left[\begin{array}{ccccc}{a}_{11}& a_{12}& 0& 0& 0\\ a_{21}& a_{22}& a_{23}& 0& 0\\ 0& a_{32}& a_{33}& a_{34}& 0\\ 0& 0& a_{43}& a_{44}& a_{45}\\ 0& 0& 0& a_{54}& a_{55}\end{array}\right]$$
以三对角矩阵为例，对其以行序列主序进行压缩存储，只需要将其中的非零元素存储到一维数组$B$ 中，$a_{ij}$的对应存储位置是$b_k$ ，$k$从 0 开始。

观察三角矩阵$A$，可以得知三角矩阵的如下特点：

• 行下标为 1 的行和行下标为$n$的行，都只有两个非零元素。
• 其他行各有三个非零元素。

对于行下标不为 1 的非零元素$a_{ij}$，在其之前的非零元素包括：

• 行下标为 1 的行有 2 个非零元素；

• 从行下标 $2$ 到行下标 $i-1$ 的行，共计 $i-2$行，每行 3 个非零元素，共计$3(i-2)$ 个元素

• $a_{ij}$所在行的行下标为$i$ ，$a_{ij}$之前的元素可能是：

  1. 若 $a_{ij}$是本行的第 1 个非零元素，则$a_{ij}$之前的所有元素数量是：$2+3(i-2)=3i-4$。并且此时，必然有$j=i-1$ 。所以 $k$ 与 $i,j$ 的关系为：
     $$k=3i-4=2i+i-1-3=2i+j-3$$

  2. 若 $a_{ij}$是本行的第 2 个非零元素，则$a_{ij}$：就是之前的所有元素数量$2+3(i-2)+1=3i-3$。并且此时，必然有$i=j$ 。所以 $k$ 与 $i,j$ 的关系为：
     $$k=3i-3=2i+i-3=2i+j-3$$

  3. 若 $a_{ij}$是本行的第 3 个非零元素，则$a_{ij}$之前的所有元素数量是：$2+3(i-2)+2=3i-2$。并且此时，必然有$j=i+1$ 。所以 $k$ 与 $i,j$ 的关系为：
     $$k=3i-2=2i+i+1-3=2i+j-3$$

综上可知，对于三对角矩阵而言，压缩存储的时候，数组中的位置$k$ 和矩阵的 $i,j$ 关系是（$i,j$是从 1 开始，$k$从 0 开始）：
$$k=2i+j-3\text{\hspace{1em}(6)}$$

例 6.6.1若对 n 阶对称矩阵$A$ 以行序为主序方式将其下三角形的元素（包括主对角线上所有元素）依次存放于一维数组 B[1..(n(n+1))/2] 中，则在 B 中确定 $a_{ij}(i<j)$的位置 k 的关系为（ ）。

A. $i\times (i-1)/2+j$\qquad B. $j\times (j-1)/2+i$\qquad C. $i\times (i+1)/2+j$\qquad D. $j\times (j+1)/2+i$

【解】

矩阵 $A$ 是 $n$阶对称矩阵，且元素$a_{ij}$中的下标关系是$i<j$，说明此元素是上三角矩阵的一个元素。题目中要求以行序为主序的方式存储下三角的元素，则与$a_{ij}$ 对称的是 $a_{ji}$ ，$a_{ij}=a_{ji}$。因而，本题转换为：确定$a_{ji},(i\ge j)$ 的位置 $k$ 的关系。

根据以行序为主序的存储方式，存储对称矩阵的下三角（含对角线）不分，则第 1 行存储$a_{11}$；第 2 行存储$a_{21},a_{22}$ ；$\cdots$ ；第 $(j-1)$ 行存储 $a_{(j-1)1},\cdots ,a_{(j-1)(j-1)}$ ；第 $j$ 行存储 $a_{j1},\cdots ,a_{ji}$。将各行所存储的元素个数求和：
$$1+2+3+\cdots +(j-1)+i=\frac{j\times (j-1)}{2}+i$$
所以，$k=\frac{j(j-1)}{2}+i$ 。

本题答案：B