稀疏矩阵.docVIP

#includestdio.h #includemalloc.h #define sMax 45 typedef int elemtype; typedef struct lnode //结构体的定义 { int i,j; struct lnode *cp,*rp; union//共用体的定义,输入输出使用同一段存储结构 { struct lnode *next; elemtype v; }mat; }link; int flag=0; //建立稀疏矩阵的函数，返回十字链表头指针 link *creatlinkmat() { link *p,*q,*head,*cp[sMax]; int i,j,k,m,n,t,s; elemtype v; printf(输入行、列，非零元素个数(m,n,t数字间用逗号分隔)); scanf(%d,%d,%d,m,n,t);//输入行、列，非零元素个数 if(mn)s=m; else s=n; head=(link *)malloc(sizeof(link)); //建立十字链表头结点，分配内存空间 head-i=m;head-j=n; cp[0]=head; //cp[]是指针数组，分别指向头结点和行、列表头结点 for(i=1;i=s;i++) //建立头结点循环链表 { p=(link *)malloc(sizeof(link)); p-i=0;p-j=0; p-rp=p;p-cp=p; cp[i]=p; cp[i-1]-mat.next=p; } cp[s]-mat.next=head; for(k=1;k=t;k++) { printf(\t 第%d个元素(行号i 列号j 值v，数字间用空格分隔):,k); scanf(%d %d %d,i,j,v); p=(link *)malloc(sizeof(link)); p-i=i;p-j=j;p-mat.v=v; q=cp[i]; while((q-rp!=cp[i])(q-rp-jj)) q=q-rp; p-rp=q-rp; q-rp=p; q=cp[j]; while((q-cp!=cp[j])(q-cp-ii)) q=q-cp; p-cp=q-cp; q-cp=p; } return head; } //插入结点函数 void insert(int i,int j,int v,link *cp[]) { link *p,*q; p=(link *)malloc(sizeof(link)); p-i=i;p-j=j;p-mat.v=v; //以下是经*p结点插入第i行链表中 q=cp[i]; while((q-rp!=cp[i])(q-rp-jj)) q=q-rp;//在第i行中找第一个列号大于j的结点*(q-rptr) //找不到时，*q是该行表上的尾结点 p-rp=q-rp; q-rp=p;//*p插入在*q之后 //以下是将结点插入第j列链表中 q=cp[j];//取第j列表头结点 while((q-cp!=cp[j])(q-cp-ii)) q=q-cp ;//在第j行中找第一个列号大于i的结点*(q-cptr) //找不到时，*q是该行表上的尾结点 p-cp=q-cp; q-cp=p;//*p插入在*q之后 } //输出十字链表的函数 void print(link *A) { link *p,*q,*r;//p是控制行q是控制列r是控制输出的格式 int k,c,t; c=A-j;//矩阵A的列数 printf(矩阵为：\n); p=A-mat.next;//p指向第一个结点（不是头结点） while(p!=A) { q=p-rp;//p指向这以一行的一个值 if(q==A-cp)break;//行或列处理完成，跳出循环 r=p;//r指向这一行的头结点 while(q!=p) { for(k=1;kq-j-(r-j);k++)//输出同一行上两非零数据间的零 printf( 0); printf(%3d,q-mat.v);//输出那个非零值 q=q-rp;//q指向这一行的下一个元素 r=r-rp;//r指向q前面的一个非零元素 } k=r-j;//k的值是某一行的最后一个非零元的列数 for(t=k;tc