稀疏矩阵的十字链表实现（C语言实现）

又要从头学数据结构了。

默默挨个实现一遍

实现的功能包括两个矩阵的加减 相乘 还有转置 

有点繁琐 

凑合看吧

#include <stdio.h>
#include <malloc.h>
typedef int ElemType;// 稀疏矩阵的十字链表存储表示 

typedef struct OLNode
{
	int i,j;	//该非零元的列和下标 
	ElemType e;		//非零元素值 
	struct OLNode *right,*down;		//该非零元所在行表和列表的后继链域 
}OLNode,*OLink;

typedef struct
{
	OLink *rhead,*chead;
	int mu,nu,tu;		//稀疏矩阵的行数、列数和非零元个数 
}CrossList;

//初始化M  
int InitSMatrix(CrossList *M)
{
	(*M).rhead=(*M).chead=NULL;
	(*M).mu=(*M).nu=(*M).tu=0;
	return 1;
}
//销毁稀疏矩阵M
int DestroySMatrix(CrossList *M)
{
	int i;
	OLNode *p,*q;
	
	for(i=1;i<=(*M).mu;i++)//按行释放节点 
	{
		p=*((*M).rhead+i);
		while(p)
		{
			q=p;
			p=p->right;
			free(q);
		}
	}
	free((*M).rhead);
	free((*M).chead);
	(*M).rhead=(*M).chead=NULL;
	(*M).mu=(*M).nu=(*M).tu=0;
	return 1;
} 

// 创建稀疏矩阵M，采用十字链表储存表示
int CreateSMatrix(CrossList *M)
{
	int i,j,k,m,n,t;
	ElemType e;
	OLNode *p,*q;
	if((*M).rhead)
		DestroySMatrix(M);
	printf("请输入稀疏矩阵的行数 列数 非零元素个数：（space）");
	scanf("%d %d %d",&m,&n,&t);
	(*M).mu=m;
	(*M).nu=n;
	(*M).tu=t;
	//初始化行链表头
	(*M).rhead=(OLink*)malloc((m+1)*sizeof(OLink));
	if(!(*M).rhead)
		exit(0);
	//初始化列链表头
	(*M).chead= (OLink*)malloc((n+1)*sizeof(OLink));
	if(!(*M).chead)
		exit(0);
	for(k=1;k<=m;k++)//初始化行头指针向量 各行链表为空链表 
		(*M).rhead[k]=NULL;
	for(k=1;k<=n;k++)// 初始化列头指针向量;各列链表为空链表 
		(*M).chead[k]=NULL;
	printf("请按任意次序输入%d个非零元的行 列 元素值:（空格）\n",(*M).tu);
	for(k=0;k<t;k++)
     {
         scanf("%d %d %d",&i,&j,&e);
        p=(OLNode*)malloc(sizeof(OLNode));
        if(!p)
            exit(0);
       p->i=i; // 生成结点 
       p->j=j;
       p->e=e;
       if((*M).rhead[i]==NULL||(*M).rhead[i]->j>j)    
       {
           // p插在该行的第一个结点处
           p->right=(*M).rhead[i];
           (*M).rhead[i]=p;
       }
       else // 寻查在行表中的插入位置 
       {
           //从该行的行链表头开始，直到找到
           for(q=(*M).rhead[i]; q->right && q->right->j < j;q = q->right)
               ;
           p->right=q->right; // 完成行插入 
           q->right=p;
       }
       if((*M).chead[j] == NULL || (*M).chead[j]->i > i) 
       {
           // p插在该列的第一个结点处
           p->down = (*M).chead[j];
           (*M).chead[j] = p;
       }
       else // 寻查在列表中的插入位置 
       {
           for(q = (*M).chead[j];q->down && q->down->i < i;q = q->down)
               ;
           p->down=q->down; // 完成列插入 
           q->down=p;
       }
   }
   return 1;
} 

// 按行或按列输出稀疏矩阵M
int PrintSMatrix(CrossList M)
{
    int i,j;
    OLink p;
    printf("%d行%d列%d个非零元素\n",M.mu,M.nu,M.tu);
    printf("请输入选择(1.按行输出 2.按列输出): ");
    scanf("%d",&i);
    switch(i)
    {
    case 1: 
        for(j=1;j<=M.mu;j++)
        {
            p=M.rhead[j];
            while(p)
            {
                printf("%d行%d列值为%d\n",p->i,p->j,p->e);
                p=p->right;
            }
        }
        break;
    case 2: 
        for(j=1;j<=M.nu;j++)
        {
            p=M.chead[j];
            while(p)
            {
                printf("%d行%d列值为%d\n",p->i,p->j,p->e);
                p=p->down;
            }
        }
    }
    return 1;
}

// 由稀疏矩阵M复制得到T
int CopySMatrix(CrossList M,CrossList *T)
{
    int i;
    OLink p,q,q1,q2;
    
    if((*T).rhead)
        DestroySMatrix(T);
    (*T).mu=M.mu;
    (*T).nu=M.nu;
    (*T).tu=M.tu;
    (*T).rhead=(OLink*)malloc((M.mu+1)*sizeof(OLink));
    if(!(*T).rhead)
        exit(0);
    (*T).chead=(OLink*)malloc((M.nu+1)*sizeof(OLink));
    if(!(*T).chead)
        exit(0);
    for(i=1;i<=M.mu;i++) // 初始化矩阵T的行头指针向量;各行链表为空链表 
        (*T).