用十字链表和一般方法分别实现稀疏矩阵的乘法和加法..doc

#includestdio.h #includemalloc.h #define Size 2501 # define Size1 51 typedef struct { int i; int j; int e;//非零元的值 }triple; //定义三元组 typedef struct { triple data[Size+1];//矩阵中的元素 int rops[Size1+1];// rops[i]为第i行元素中的首非零元在data[]中的序号 int mu;//行数 int nu;//列数 int tu;//非零元数 } juzhen;//定义矩阵 typedef struct node// 定义十字链表元素 { int i,j,e; struct node *right,*down;// 该非零元所在行表和列表的后继元素 }node,*link; typedef struct // 定义十字链表对象结构体 { link *rhead,*chead;//行和列的头指针 int m,n,t;// 系数矩阵的行数，列数，和非零元素个数 }crosslist; void createcross(crosslist M)//建立十字链表 { int i,j,e,k; node *p,*q; printf(输入行，列和非零元数，空格隔开：\n); scanf(%d %d %d,M.m,M.n,M.t); M.rhead=(link *)malloc((M.m+1)*sizeof(link));//给行和列的头指针分配内存 M.chead=(link *)malloc((M.n+1)*sizeof(link)); for(k=1;k=M.m;k++)//初始化行，列的头指针 M.rhead[k]=NULL; for(k=1;k=M.n;k++) M.chead[k]=NULL; printf(输入非零元的行，列和值，空格隔开：\n); for(k=1;k=M.t;k++)//输入非零元 { scanf(%d %d %d,i,j,e); p=(node *)malloc(sizeof(node)); 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; } if(M.chead[j]==NULL||(M.chead[j]-ii))//插入元素所在列无非零元或首非零元的行标大于插入元素的行标 { p-down=M.chead[j]; M.chead[j]=p; } else { for(q=M.chead[j];(q-down)q-down-ii;q=q-down);//空循环找到第一个行标大于或等于插入元素行标的元素 p-down=q-down; q-down=p; } } } void printcross(crosslist A)//输出十字链表 { if(A.m==0) printf(十字链表为空!\n); else { printf(十字链表为:\n); int i,j; for(i=1;i=A.m;i++) { link p=A.rhead[i]; for(j=1;j=A.n;j++) { if((p)(j==p-j)) { printf(%5d,p-e); p=p-right; } else printf(%5d,0); } printf(\n); } } printf(\n); } crosslist addcross() { printf(十字链表加法:\n); crosslist a,b;// 创建两个十字链表对象，并初始化 createcross(a); createcross(b); node *pre,*h[51],*pa,*pb,*q;//定义辅助指针，pa，pb分别为a,b当前比较的元素，pre为pa的前驱元素 int i,