数据结构实验报告之三元组顺序表存储的稀疏矩阵练习

一、实验内容

1.实验目的

为了更加深入了解三元组以及稀疏矩阵，学会快速转置稀疏矩阵的算法设计并以此为基础思考其他相关问题。本次实验主要以解决某特定6阶稀疏矩阵三对角线元素之和以及转置该6阶矩阵为目的而进行算法设计以及代码实现。

2.实验内容

(1)构建一个三元组顺序表存储的稀疏6阶方阵
(2)求解该方阵中三条对角线上的元素之和
(3)实现该方阵的快速转置
(4)实现代码并输出相关数据                                                                              

---

一．要实现（1）只需要给与指定的行、列以及元素值即可。本次实验原稀疏矩阵设置为S，转置矩阵设置为T，下面展示的是S以及S对应的三元组：

S
0	0	0	0	0	0
1	0	0	0	0	0
0	0	0	2	0	0
0	0	0	0	0	0
0	0	3	0	0	0
0	0	0	0	0	4

S.mu=6	S.nu=6	S.tu=4
S.data	i	j	e
0
1	2	1	1
2	3	4	2
3	5	3	3
4	6	6	4

---

即转置后的T方阵以及对应的三元组应当如下所示：

T
0	1	0	0	0	0
0	0	0	0	0	0
0	0	0	0	3	0
0	0	2	0	0	0
0	0	0	0	0	0
0	0	0	0	0	4

T.mu=6	T.nu=6	T.tu=4
T.data	i	j	e
0
1	1	2	1
2	3	5	3
3	4	3	2
4	6	6	4

---

二.为实现（2）中的三对角线总和，本此实验设计了一个sumDiagnoal3函数以求和。

三.快速转置算法为课上所学FastTransposeTSMatrix函数，并给与了更加详细的描述。

四.实验输出见下文（三、实验结果）。

二、实现过程

本次实验环境：macOS Mojave, Xcode, gcc, C99

根据实验要求，本次实验依次设计了两个函数以完成三对角线求和方阵转置，并给与较为清晰的输出。下面是实验过程：

---

1.编译预处理以及使用类别名

编译预处理以及使用类别名，并没有使用到过大的矩阵但还是将MAXSIZE设置的足够大，创建了三元组(Triple)结构体以及稀疏矩阵(TSMatrix)结构体如下所示:

#include <stdio.h>
#define MAXSIZE 12500
#define OK 1
#define ERROR 0
typedef int Status, ElemType;
typedef struct {
   
   int i, j;//该非零元素的行下标和列下标
   ElemType e;
} Triple;
typedef struct {
   
   Triple data[MAXSIZE+1];//非零元三元组表，data[0]未用
   int mu, nu, tu;//矩阵的行数，列数和非零元个数
} TSMatrix;//稀疏矩阵

---

2.sumDiagnoal3()函数

计算三条对角线元素之和的函数(sumDiagnoal3)，如下所示：
初实现求和算法时并没有仔细考虑，但是突然产生一个念想，要是这次实验要我们完成的是五条对角线呢？于是带着这样的想法，我决定吧算法单独作为一个小函数拿出来。

int sumDiagonal3(TSMatrix S) {
   //求三条对角线上元素之和
   int i, sum=0;
   for(i=1; i<=S.tu; i++) {
   
      if((S.data[i].i-S.data[i].j)<=1 && (S.data[i].j-S.data[i].i)<=1)
         sum += S.data[i].e;
   }
   return sum;
}

---

3.FastTransposeTSMatrix()函数

快速转置算法(FastTransposeTSMatrix)，如下所示：
由于该函数实现机制过于巧妙，故在编写代码时我添加了很多的注释以帮助自己读懂。不妨下面我再解释几行较为难懂的代码，一遍日后复习时需要使用：

Status FastTransposeTSMatrix(TSMatrix M, TSMatrix *T) {
   
   //采用三元组顺序表存储表示，求稀疏矩阵M的转置矩阵T。
   T->mu=M.nu; T->nu=M.mu;T->tu=M.tu;
   int col,t;//col指示列,t指示元素个数
   int num[M.mu], cpot[M.mu];//分别表示M矩阵中第col列非零元素个数 和
   if(T