最近略忙,今天为大家献上数组一节的拓展知识,我们可以从中加深对数组存储对理解,同时温习数组降低存储量的方法,我们会以矩阵存储为典例。
主干知识,可以对照黑体字来进行自测~
数组
基本概念
一维数组是纯粹的线性表,数组的元素类型就是线性表的元素类型;
二维数组则可以看成“元素是一维数组”的线性表。
其中
二维数组也可以看成是的一个阵列(矩阵)
同理,三维数组也可以看成是二维数组的线性表,N维数组看成元素是N-1维数组的线性表。
数组的主要操作有如下4个:
InitArray(&A,n,bound1,…,boundn)
操作结果:初始化一个n 维数组,各维长度由 boundl,…,boundn 给定。
参数说明:bound1,…,boundn 合法正整数。
DestroyArray(&A)
操作结果:销毁数组 A。
参数说明:数组 A 已经存在。
Value(A,&e, index1,…, indexn)
操作结果:把e赋值为A 的指定下标的元素值。
参数说明:数组A已经存在,下标不越界。
Assign(&A,e, index1,…, indexn)
操作结果:把 e的值赋给 A 的指定下标的元素。
参数说明:数组 A已经存在,下标不越界。
数组的顺序存储表示
由于数组可能是多维的,而计算机的内存只能是一维的连续空间,必须把多维的数组转换成一个一维的序列才可以存储。
大多数语言以“行”为主序存储
假设二维数组中每个元素占用个存储单元,则该数组中任意一个元素在存储映像区的存储地址可以这样确定:
其中是的存储地址,是的存储地址,也就是数组A存储映像区的首地址,称为数组A的基地址。
同样三维数组中任意一个元素在存储映像中地址可以这样确定:
显然,数组元素的地址是其下标的线性函数,存取每个元素所进行的运算及次数都是一样的,称具有这种特点的存储结构为随机存取存储结构。
矩阵的压缩存储
矩阵在科学计算领域有广泛的应用。随计算机应用和科学领域的发展,出现了很多高阶矩阵的计算问题,有的甚至达到几十万阶,上千亿的数组元素。然而这些矩阵很多都是有规律的,或者包含大量的零值元素。
对这类矩阵进行压缩,不仅可以节约存储空间,而且可以避免大量零值元素参与计算。
如果相同值的元素或零值元素在矩阵中的分布具有一定的规律和特征,这类矩阵被称为特殊形状的矩阵。如果没有什么特征,则称为随机稀疏矩阵。
特殊形状矩阵通常有三角形矩阵,条带状矩阵,对称阵等
三角形矩阵是指矩阵中上三角或下三角是全零值元素或同值元素。
条带状矩阵则是只有主对角线附近的若干条对角线上含有非零元素。
利用二维数组存储这些大多是零值元素的矩阵,空间浪费显然巨大,我们知道二维数组在实际存储时也是转换为一维的序列存储的,不管是行主序还是列主序,都可以建立一个元素存储位置和下标的函数关系,利用该函数可以通过下标对数组元素实现随机存取。
那么既然特殊形状矩阵有这么多零值元素,我们能否不存储零值元素,通过修改函数关系而仍然实现对数组元素的随机存取呢?
假设 N 阶方阵(行列相等)A 的元素满足,则称之为对称阵。对称阵中的一对元素可以只存储一个存储空间,于是可以把个元素压缩存储到个存储单元里。假设以来存储这些元素,按照行主序的次序把对称矩阵下半区(包括主对角线)存储到 中。
得到的存储地址满足(L是每个元素占用的存储单元数):
对于三角形状矩阵,如下三角阵,完全可以参照对称阵来给出函数关系,只是 i<j 时元
素值是 0,不再需要以上第二个计算式。
另一种带状特殊矩阵也可以采用类似的思想来存储在一个一维的数组中,下面考虑三
对角矩阵,即只有主对角线及其两侧的对角线上可以有非零值。
带状矩阵存储位置与坐标的函数关系是
在的矩阵中,若非零元素有t个,我们定义
为矩阵的稀疏因子,通常认为时为稀疏矩阵。
由于无特值可循,因此需要记录非零元素的位置。这样就构成了一个三元组。
三元组顺序表存储的实现如下:
#define MAX_SIZE 1000
typedef int ElemType;
typedef struct {
int i, j;
ElemType e;
}Triple;
typedef struct {
Triple data[MAX_SIZE];
int mu, nu, tu;
}TSMatrix;
原先矩阵转置的算法为:
void transpose(ElemType M[][], ElemType T[][], int m, int n) {
for (int i = 0; i < m; i++) {
for (int j = 0; j < n; j++) {
T[j][i] = M[i][j];
}
}
}
很明显这个算法的时间复杂度为
接下来我们使用“按位就座”的思想,有了rpos数组,就可以一次完成所有的非零元的排放。
#define MAXN 100
void createrpos(TSMatrix M) {
for (int col = 0; col < M.nu; col++)num[col] = 0;
for (t = 0; t < M.tu; t++) {
++num[M.data[t].j];//统计各列非零元
}
rpos[0] = 0;
for (int col = 1; col < M.nu; col++) {
rpos[col] = rpos[col - 1] + num[col - 1];//累计num数组获得rpos
}
}
下面给出三元顺序表压缩存储的矩阵M求其转置矩阵T的算法(略去了col的int初始化和q,rpos等元素初始化)
void TRansposeTS(TSMatrix M, TSMatrix& T) {
T.mu = M.nu;
T.nu = M.mu;
T.tu = M.tu;
if (!M.tu) return;
createrpos(M);
for (int p = 0; p < M.tu; p++) {
col = M.data[p].j;
q = rpos[col];
T.data[q].i = M.data[p].j;
T.data[q].j = M.data[p].i;
T.data[q].e = M.data[p].e;
++rpos[col];
}
}
思想:查询当前元素在M中的列,在T中的行
交换i,j坐标,赋元素值
然后rpos指向T中下一个col行的非零元素位置
从上述算法可以看出,包括 createrpos 算法在内,共有 4 个串行工作的单循环,循环次
数分别是2 个 nu 次和 2 个 tu 次,因此总的时间复杂度是 。其效率大大
高于二维数组存储的矩阵转置时间复杂度 。
由于 rpos 数组大大简化了三元组存储矩阵的运用,我们常常在建立三元组表的时候就
建立起这样的数组,表示矩阵“每一行的非零元在三元组顺序表中起始位置”的信息。这种
“带行链接信息”的三元组表被称为逻辑链接的顺序表。
采用三元组存储来表示矩阵时,一般采用顺序表来存储三元组,但当在矩阵运算产生元组的增减时,就需要对顺序表进行元素的插入和删除了,这时也可以考虑使用链表来存储三元组。
对于稀疏矩阵的压缩,除了使用三元组存储外,有时还使用十字链表来存储。十字链表是这样表示的:为一个 m × n 的矩阵创建 m + n 个单链表,m 个行链表和n 个列链表,每个行链表包含该行的所有非零元,同样,每个列链表包含该列的所有非零元。m 个行链表和n个列链表的头指针分别构成一维数组,属于指针数组。这样每个非零元素既属于一个行链表也属于一个列链表。因此,每个非零元素都是两个链表的交叉结点,故称十字链表。十字链表的结点包含 5 个域,具体描述如下:
typedef struct OLNode {
int i, j;
ElemType e;
OLNode* rownext, * colnext;
}OLNode,*OLink;
typedef struct {
OLink* rhead, * chead;
int m, n, t;
}CrossList;
删繁就简三秋树,正确使用数组能够使我们的空间存储大大减小,熟练掌握也有利于我们处理相对较大的数据集,同时降低时间复杂度,提升效率
下一讲,我将带来重要数据结构“树”的一部分内容
今天肝得有点晚,睡觉!
附下一讲链接:【数据结构】五分钟自测主干知识(八)
未完待链接
标签:存储,元素,矩阵,链表,五分钟,数组,自测,数据结构,col From: https://blog.csdn.net/2301_79853895/article/details/136661683