矩阵的压缩存储.pptVIP

计算机科学与工程系 计算机科学与工程系 第四章 矩阵的压缩存储 多维数组 数组的定义及其基本操作 数组的存储结构 矩阵的压缩存储 特殊矩阵的压缩存储 稀疏矩阵的压缩存储 4.1 数组的定义及其基本操作 3) 数组的基本运算 4.1.2 数组的顺序存储 一个一维数组，一旦第一个元素a0的存储地址Loc(a0)确定，而每个元素所占用的存储空间大小为l, 则第i个元素的地址可以由以下公式计算： 二维数组与高维数组 1）存储方式 计算机的存储结构是一维的，因而多维数组必须按某种次序排成一个线性序列加以存储。 2）任意元素存储地址的计算 以二维数组行优先顺序存储为例：假设每个元素占用l单元。 4.2 特殊矩阵的压缩存储 4.2.1 对称矩阵 4.2.2 三角矩阵 4.3 稀疏矩阵的压缩存储 4.3.1 三元组表 例： 方法二：快速转置算法 思想：即按A.data[]中三元组次序转置，转置结果放入B.data[]中恰当位置。此法关键是要预先确定A中每一列第一个非零元在B.data[]中位置，为确定这些位置，转置前应先求得A的每一列中非零元个数。 增设辅助数组 num 和 cpot 。 num[col]:表示矩阵A 第col列中非0元素的个数。 cpot[col]:表示第col列第一个非0元素在B.data[]中位置 Cpot[1]=1 Cpot[col]=cpot[col-1]+num[col-1] 2=col=A.nu 2) 十字链表 本章小结 数组 数组的顺序存储（行优先，列优先） 逻辑地址 ? 物理地址的计算； 矩阵的压缩存储 对称矩阵、上三角或下三角等特殊矩阵压缩存储时的地址转换公式 稀疏矩阵的压缩存储：存储特点、三元组表表示 十字链表用C语言描述如下： typedef struct OLNode{ int i, j; /*矩阵非零元素的行、列值*/ ElemType e; /*非零元素的值*/ struct OLNode *right, *down; /*该非零元所在行、列的后继链域*/ } OLNode, *Olink; typedef struct{ Olink *rhead, *chead; /*行、列表指针向量基址*/ int mu,nu,tu; /*矩阵的行、列数及非零元数个数*/ }CrossList; 例：用十字链表表示下面的矩阵。 1 1 3 1 4 5 ^ 2 2 -1 3 1 2 ^ ^ ^ ^ M.chead M.rhead ^ ^ 十字链表的程序实现 int CreatLinkMatrix(CrossList M) { if(M) free(M); scanf（m，n，t）; //输入矩阵行、列值和非零元个数 M.mu=m; M.nu=n; M.tu=t; if(!M.rhead=(Olink*)malloc((m+1)*sizeof(Olink)))) return OVERFLOW; /*内存分配失败，返回溢出信息*/ if(!M.chead=(Olink*)malloc((n+1)*sizeof(Olink)))) return OVERFLOW; M.rhead[]=M.chead[]=NULL;//初始化行列头指针向量 for( scanf(i,j,e); i!=0;scanf(i,j,e)) { //按任意次序输入非零元 if(!(p=(OLNode*)malloc(sizeof(OLNode)))) return OVERFLOW; p-i=i; p-j=j; p-e=e; //生成结点 if(M.rhead[i] = =NULL||M.rhead[i]-jj)) { p-right=M.rhead[i]; M.rhead[i]=p;} //生成第一个结点 else { for(q=M.rhead[i]; (q-right) q-right-jj; q=q-right); //寻找插入位置 p-right=q-right; q-right=p;