数组和广义表新幻灯片.ppt
《数组和广义表新幻灯片.ppt》由会员分享,可在线阅读,更多相关《数组和广义表新幻灯片.ppt(120页珍藏版)》请在淘文阁 - 分享文档赚钱的网站上搜索。
1、数组和广义表新第1页,共120页,编辑于2022年,星期六5.1 数组的定义数组的定义5.3 矩阵的压缩存储矩阵的压缩存储 5.2 数组的顺序表示和实现数组的顺序表示和实现5.4 广义表的定义广义表的定义5.5 广义表的存储结构广义表的存储结构5.7 广义表操作的递归函数广义表操作的递归函数5.6 m m元多项式的表示元多项式的表示本章主要内容:本章主要内容:第2页,共120页,编辑于2022年,星期六5.1 数组的定义数组的定义ADT Array 数据对象数据对象:Daj1,j2,.,ji,jn|ji=0,.,bi-1,i=1,2,.,n n(0)称为数组的维数,bi是数组第i维的长度,ji
2、是数组元素的第i维下标,数据关系数据关系:RR1,R2,.,Rn Ri|0 jk bk-1,1 k n 且k i,0 ji bi-2,i=2,.,n 基本操作基本操作:ADT Array 第3页,共120页,编辑于2022年,星期六二维数组的定义二维数组的定义:数据对象数据对象:D=aij|0ib1-1,0 jb2-1数据关系数据关系:R=ROW,COL ROW=|0ib1-2,0jb2-1 COL=|0ib1-1,0 jb2-2基本操作基本操作:InitArray(&A,n,bound1,.,boundn)DestroyArray(&A)Value(A,&e,index1,.,indexn)
3、Assign(&A,e,index1,.,indexn)第4页,共120页,编辑于2022年,星期六 InitArray(&A,n,bound1,.,boundn 操作结果:操作结果:若维数 n 和各维长度合法,则构造相应的数组A,并返回OK。DestroyArray(&A)操作结果:操作结果:销毁数组A。Value(A,&e,index1,.,indexn)初始条件:初始条件:A是n维数组,e为元素变量,随后是n 个下标值。操作结果:操作结果:若各下标不超界,则e赋值为所指定的A 的元素值,并返回OK。第5页,共120页,编辑于2022年,星期六 Assign(&A,e,index1,.,i
4、ndexn)初始条件:初始条件:A是n维数组,e为元素变量,随后是n 个下标值。操作结果:操作结果:若下标不超界,则将e的值赋给所指定的A的元素,并返回OK。第6页,共120页,编辑于2022年,星期六5.2 数组的顺序表示和实现数组的顺序表示和实现 类型特点类型特点:1)只有引用型操作,没有加工型操作;2)数组是多维的结构,而存储空间是一个一维的结构。有两种顺序映象的方式有两种顺序映象的方式:1)以行序为主序(低下标优先);2)以列序为主序(高下标优先)。第7页,共120页,编辑于2022年,星期六例如:例如:称为基地址基地址或基址以以“行序为主序行序为主序”的存储映象的存储映象二维数组A中
5、任一元素ai,j 的存储位置 LOC(i,j)=LOC(0,0)+(b2ij)La0,1a0,0a0,2a1,0a1,1a1,2二维数组存储结构a0,1a0,0a0,2a1,0a1,1a1,2L第8页,共120页,编辑于2022年,星期六 推广到一般情况,可得到推广到一般情况,可得到 n 维数组数据元素存储位维数组数据元素存储位置的映象关系置的映象关系 称为 n 维数组的映象函数。数组元素的存储位置是数组元素的存储位置是其下标的线性函数。其下标的线性函数。其中 cn=L,ci-1=bi ci,1 i n。LOC(j1,j2,.,jn)=LOC(0,0,.,0)+ci ji i=1n第9页,共1
6、20页,编辑于2022年,星期六/数组的顺序存储表示数组的顺序存储表示#include#include /标准头文件,提供宏标准头文件,提供宏va_startva_start、va_argva_arg和和va_end,va_end,用于存取用于存取变长参数表变长参数表#define MAX_ARRAY_DIM 8#define MAX_ARRAY_DIM 8 /假设数组维数的最大值为假设数组维数的最大值为8 8typedef structtypedef struct ElemType*base;ElemType*base;/数组元素地址,由数组元素地址,由InitArrayInitArray分
7、配分配 int dim;int dim;/数组维数数组维数 int*bounds;int*bounds;/数组维界基址,由数组维界基址,由InitArrayInitArray分配分配 int*constants;int*constants;/数组映像函数常量基址,由数组映像函数常量基址,由InitArrayInitArray分配分配Array;Array;第10页,共120页,编辑于2022年,星期六/基本操作的函数原型说明基本操作的函数原型说明Status InitArray(Array&A,int dimStatus InitArray(Array&A,int dim););/若维数若维数
8、dimdim和随后的各维长度合法,则构造相应的数组和随后的各维长度合法,则构造相应的数组A A,并返回,并返回OKOK。int DestroyArray(Array&A);int DestroyArray(Array&A);/销毁数组销毁数组A A int StrCompare(Array&A,ElemType&e,int StrCompare(Array&A,ElemType&e,)/A/A是是n n维数组,维数组,e e为元素变量,随后是为元素变量,随后是n n各下标值。各下标值。/若各下标不超界,则将若各下标不超界,则将e e的值赋给所指定的的值赋给所指定的A A的元素,并返的元素,并返
9、回回OKOK。第11页,共120页,编辑于2022年,星期六/基本操作的算法描述基本操作的算法描述Status InitArray(Array&A,int dim,);/若维数若维数dimdim和随后的各维长度合法,则构造相应的数组和随后的各维长度合法,则构造相应的数组A A,并返回,并返回OKOK。if(dimMAX.ARRAY_DIM)return ERROR;A.dim=dim;A.bounds=(int*)malloc(dim*sizeof(int);if(!A.bounds)exit(OVERFLOW);/若各维长度合法,则存入若各维长度合法,则存入A.boundsA.bounds,
10、并求出,并求出A A的元素总数的元素总数elemtotalelemtotal elemtotal=1;va_start(ap,dim);/ap为为va_list类型,是存放变长参数表信息的数组类型,是存放变长参数表信息的数组 for(i=0;idim;+i)A.boundsi=va.arg(ap,int);if(A.bounds i=0;-i)A.constantsi=A.boundsi+1 A.constantsi+1;return OK;第13页,共120页,编辑于2022年,星期六Status DestroyArray(Array&A)/销毁数组A if(!A.base)return E
11、RROR;free(A.base);A.base=NULL;if(!A.bounds)return ERROR;free(A.bounds);A.bounds=NULL;if(!A.constants)return ERROR;free(A.constants);A.constants=NULL;return OK;第14页,共120页,编辑于2022年,星期六 Status Locate(Array A,va_list ap,int&off)/若ap指示的各下标值合法,则求出该元素 在A中相对地址off off=0;for(i=0;iA.dim;+i)ind=va_arg(ap,int);i
12、f(ind=A.boundsi)return OVREFLOW;off+=A.constantsi*ind;return OK;第15页,共120页,编辑于2022年,星期六 Status Value(Array A,ElemType&e,)/A/A是是n n维数组,维数组,e e为元素变量,随后是为元素变量,随后是n n各下标值。各下标值。/若各下标不超界,则若各下标不超界,则e e赋值为所指定的赋值为所指定的A A的元素的元素 值,并返回值,并返回OKOK va_start(ap,e);if(result=Loate(A,ap,off)=0)return result;e=*(A.base
13、+off);return OK;第16页,共120页,编辑于2022年,星期六 Status Assign(Array&A,ElemType e,)/A/A是是n n维数组,维数组,e e为元素变量,随后是为元素变量,随后是n n各下标值。各下标值。/若各下标不超界,则若各下标不超界,则e e赋值为所指定的赋值为所指定的A A的元素值,并返回的元素值,并返回OKOK va_start(ap,e);if(result=Loate(A,ap,off)=0)return result;*(A.base+off)=e;return OK;第17页,共120页,编辑于2022年,星期六5.3 5.3 矩
14、阵的压缩存储矩阵的压缩存储 在数值分析中经常出现一些阶数很高的矩阵,同时在矩阵中有许多值相同的元素或零元素。有时为节省空间,可以对这类矩阵进行压缩存储压缩存储。假设值相同的元素或零元素在矩阵中的分布有一定规律,则称此类矩阵为特殊矩阵特殊矩阵;反之,称为稀疏矩阵稀疏矩阵。第18页,共120页,编辑于2022年,星期六 以常规方法,即以二维数组表示高阶的稀疏矩阵时产生的问题问题:1)零值元素占了很大空间;2)计算中进行了很多和零值的运算,遇除法,还需判别除数是否为零。1)尽可能少存或不存零值元素;解决问题的原则解决问题的原则:2)尽可能减少没有实际意义的运算;3)操作方便。即:能尽可能快地找到与下
15、标值(i,j)对应的元素,能尽可能快地找到同一行或同一列的非零值元。第19页,共120页,编辑于2022年,星期六 若n阶矩阵A中的元素满足下述性质 aij=aji 1=i,j=j 对角线下方对角线下方j(j-1)2+i-1 当当ij 对角线上方对角线上方 对于任意给定一组下标(i,j),均可在sa中找到矩阵元素aij,反之,对所有的k=0,1,2,(n(n+1)/2-1,都能确定sak 中的元在矩阵中的位置(i,j)。由此称san(n+1)/2 为n阶对称矩阵A的压缩存储第21页,共120页,编辑于2022年,星期六 这种压缩存储的方法同样也适应于三角矩阵和对角矩阵。5.3.2 5.3.2
16、稀疏矩阵稀疏矩阵 假设 m 行 n 列的矩阵含 t 个非零元素,则称 为稀疏因子稀疏因子。通常认为通常认为 0.05 的矩阵为稀疏矩阵。的矩阵为稀疏矩阵。何谓稀疏矩阵?何谓稀疏矩阵?nmt=第22页,共120页,编辑于2022年,星期六/稀疏矩阵的三元组顺序表存储表示稀疏矩阵的三元组顺序表存储表示#define MAXSIZE 12500/最大非零元个数 typedef struct int i,j;/该非零元的行下标和列下标 ElemType e;/该非零元的值 Triple;/三元组类型三元组类型一、三元组顺序表一、三元组顺序表typedef union Triple dataMAXSIZ
17、E+1;/非零元三元组表 int mu,nu,tu;/矩阵的行数、列数和非零元个数 TSMatrix;/稀疏矩阵类型稀疏矩阵类型第23页,共120页,编辑于2022年,星期六如何求稀疏矩阵的转置矩阵?如何求稀疏矩阵的转置矩阵?-028003600070500140-005280000007143600转置例如,稀疏矩阵第24页,共120页,编辑于2022年,星期六用常规的二维数组表示时的算法用常规的二维数组表示时的算法 其中,其中,M为原稀疏矩阵,为原稀疏矩阵,T为转置后的稀疏矩阵,为转置后的稀疏矩阵,其时间复杂度为其时间复杂度为:O(munu)for(col=1;col=nu;+col)fo
18、r(row=1;row=mu;+row)Tcolrow=Mrowcol;第25页,共120页,编辑于2022年,星期六用用“三元组三元组”表示时如何实现?表示时如何实现?1 2 141 5 -52 2 -73 1 363 4 282 1 145 1 -52 2 -71 3 364 3 28a.datab.data第26页,共120页,编辑于2022年,星期六 假设a和b是TSMatrix型的变量,分别表示矩阵M和T。可以有两种处理方法:(1)按照b.data中三元组的次序依次在a.data中找到相应的三元组进行转置。换句话说,按照矩阵M的列序来进行转置。为了找到M的每一列中所有的非零元素,需要
19、对其三元组表a.data 从第一行起整个扫描一遍,由于a.data 是以M的行序为主序来存放每个非零元的,由此得到的恰是b.data应有顺序。具体算法如具体算法如5.1所示所示第27页,共120页,编辑于2022年,星期六Status TransposeSMatrix(TSMatrix M,TSMatrix&T)/采用三元组表存储表示,求稀疏矩阵M的转置矩阵T T.mu=M.nu;T.nu=M.mu;T.tu=M.tu;if(T.tu)q=1 for(col=1;col=M.nu;+col)for(p=1;p=M.tu;+p)if(m.datap.j=col)T.dataq.i=M.datap
20、.j;T.dataq.j=M.datap.i;T.dataq.e=M.datap.e;+q;/if return OK;/TransposeSMatrix算法算法 5.15.1第28页,共120页,编辑于2022年,星期六 该算法主要工作是在p和col的两重循环中完成的,故算法的时间复杂度为O(nu*tu)O(nu*tu)。当非零元的个数tu和mu*nu同数量级时,算法5.1的时间复杂度就为O(nu*nuO(nu*nu2 2)了。(2)按照a.data中三元组的次序进行转置,并将转置后的三元组置入b中恰当的位置。如果能预先确定矩阵M中每一列(即T中每一行)的第一个非零元在b.data中的准确位
21、置,则对a.data中的三元组依次转置时,便可直接放到b.data中恰当的位置上。第29页,共120页,编辑于2022年,星期六 需附设num和cpot两个向量。numcol表示矩阵M中第col列中非零元的个数,cpotcol指示M中第col列的第一个非零元在b.data中的恰当位置。cpot1=1cpot1=1;cpotcol=cpotcol-1+numcol-1cpotcol=cpotcol-1+numcol-1第30页,共120页,编辑于2022年,星期六 例如:例如:下面稀疏矩阵的三元组表示对应各向量值如下:cpot1=1;for(col=2;col=M.nu;+col)cpotcol
22、=cpotcol-1+numcol-1;col12345Numcol12011Cpotcol12445该转置方法称为快速转置,如算法5.2所示。第31页,共120页,编辑于2022年,星期六Status FastTransposeSMatrix(TSMatrix M,TSMatrix&T)T.mu=M.nu;T.nu=M.mu;T.tu=M.tu;if(T.tu)for(col=1;col=M.nu;+col)numcol=0;for(t=1;t=M.tu;+t)+numM.datat.j;/求M中每一 cpot1=1;/列非零元个数 for(col=2;col=M.nu;+col)cpotc
23、ol=cpotcol-1+numcol-1;for(p=1;p=M.tu;+p)/if return OK;/FastTransposeSMatrix 转置矩阵元素算法算法 5.25.2第32页,共120页,编辑于2022年,星期六Col=M.datap.j;q=cpotcol;T.dataq.i=M.datap.j;T.dataq.j=M.datap.i;T.dataq.e=M.datap.e;+cpotcol 转置矩阵元素:第33页,共120页,编辑于2022年,星期六 分析算法FastTransposeSMatrix的时间复杂度:时间复杂度为时间复杂度为:O(M.nu+M.tu)for(
24、col=1;col=M.nu;+col)for(t=1;t=M.tu;+t)for(col=2;col=M.nu;+col)for(p=1;p=M.tu;+p)第34页,共120页,编辑于2022年,星期六 三元组顺序表又称有序的双下标法有序的双下标法,它的特点是,非零元在表中按行序有序存储,因此便于进行依行顺序处便于进行依行顺序处理的矩阵运算理的矩阵运算。然而,若需随机存取某一行中的非零元,则需从头开始进行查找。为便于随机存取某一行中的非零元,则需知道每一行的第一个非零元在三元组表中的位置。为此,可将上节快速转置矩阵算法中建立的,指示行信息的辅助数组cpot固定在稀疏矩阵的存储结构中。称这种
25、带行链接信息的三元组表为行逻辑链接的顺序表。二、行逻辑联接的顺序表二、行逻辑联接的顺序表第35页,共120页,编辑于2022年,星期六#define MAXMN 500 typedef struct Triple dataMAXSIZE+1;/非零元三元组表 int rposMAXRC+1;/各行第一个非零元的位置表 int mu,nu,tu;/矩阵的行数、列数和非零元个数 RLSMatrix;/行逻辑链接顺序表类型行逻辑链接的顺序表的类型描述如下:第36页,共120页,编辑于2022年,星期六例如:给定一组下标,求矩阵的元素值ElemType value(RLSMatrix M,int r,
- 配套讲稿:
如PPT文件的首页显示word图标,表示该PPT已包含配套word讲稿。双击word图标可打开word文档。
- 特殊限制:
部分文档作品中含有的国旗、国徽等图片,仅作为作品整体效果示例展示,禁止商用。设计者仅对作品中独创性部分享有著作权。
- 关 键 词:
- 数组 广义 幻灯片
限制150内