/* 3 0 0 5 M= 0 -1 0 0 (M矩阵) 2 0 0 0 M.chead(colhead) 1 2 3 4 M.rhead(rowhead) 1 1 1 3 1 4 5 2 2 2 -1 3 3 1 -2 4 */ #include<stdio.h> #include<stdlib.h> //宏定义 #define OK 1 #define ERROR 0 #define OVERFLOW -2 typedef int Status; typedef struct OLNode //替换一个结构体的名字 { int i,j,e; //矩阵三元组i代表行 j代表列 e代表当前位置的数据 struct OLNode *right,*down; //指针域 右指针 下指针 }OLNode,*OLink; //结构体的新名 typedef struct { OLink *rhead,*chead; //行和列链表头指针由CreateSMatrix分配 int mu,nu,tu; //稀疏矩阵的行数,列数和非零元的个数 }CrossList; Status CreateMatrix_OL(CrossList *M); void main() { CrossList M; int i; OLNode *p; CreateMatrix_OL(&M); printf("M:\n---------------------\ni\tj\te\n---------------------\n"); for (i=1;i<=M.nu;i++) { if (NULL!=M.chead[i]) { p=M.chead[i]; while (NULL!=p) { printf("%d\t%d\t%d\n",p->i,p->j,p->e); p=p->down; } } } getchar(); getchar(); } Status CreateMatrix_OL(CrossList *M) //创建系数矩阵M,采用十字链表存储表示 { int m,n,t; //定义矩阵的行,列,非零元的数量 int i,j,e; //定义中间变量 OLNode *p,*q; //定义中间变量 // if(M) //{ // free(M); //如果矩阵已经存在,那么释放矩阵 //} scanf("%d%d%d",&m,&n,&t); //输入矩阵的行列及非零元的数量 M->mu=m; M->nu=n; M->tu=t; //初始化矩阵的行列及非零元的数量 if(!(M->rhead=(OLink*)malloc((m+1)*sizeof(OLink)))) { exit(OVERFLOW); //初始化矩阵的行链表 } if(!(M->chead=(OLink*)malloc((n+1)*sizeof(OLink)))) { exit(OVERFLOW); //初始化矩阵的列链表 } for(i=1;i<=m;i++) { M->rhead[i]=NULL; //初始化行 } for(j=1;j<=n;j++) { M->chead[j]=NULL; //初始化列 } for(scanf("%d%d%d",&i,&j,&e);0!=i;scanf("%d%d%d",&i,&j,&e)) //输入三元组 直到行为0结束 { if(!(p=(OLNode*)malloc(sizeof(OLNode)))) { exit(OVERFLOW); //动态生成p } p->i=i; p->j=j; p->e=e; //初始化p if(NULL==M->rhead[i]||M->rhead[i]->j>j) { p->right=M->rhead[i]; M->rhead[i]=p; } else { for(q=M->rhead[i];(q->right)&&q->j>j;q=q->right); p->right=q->right; q->right=p; } if(NULL==M->chead[j]||M->chead[j]->i>i) { p->down=M->chead[j]; M->chead[j]=p; } else { for (q=M->chead[j];(q->down );q=q->down); p->down=q->down; q->down=p; } } return OK; }