稀疏矩阵的十字链表的建立

系数矩阵的十字链表表示法，在我的学习过程中，是个不怎么好掌握的地方，难点就是它的具体建立过程，废话不多说，这是我学习时候的代码，注释很详细，欢迎指正学习和交流

/*
建立稀疏矩阵的十字链表
*/ 
typedef struct node
{ int row, col;
struct node *down , *right;
union v_next
{ datatype v;
struct node *next;
}
} MNode，*MLink;



MLink CreatMLink() /*返回十字链表的头指针*/
{ 
MLink H;//指针 

Mnode *p,*q,*hd[s+1];//s的大小是矩阵行数以及列数的较大值， 
/*p，q为辅助操作指针，hd为辅助数组，存储的是头节点的地址，这样便于行或者列的随机访问*/ 
int i,j,m,n,t;
/*i，j是读取矩阵元素的行和列，m，n是整个矩阵的行数目和列数目，t是非零项的数目*/ 
datatype v;
/*矩阵所存储的值*/ 
scanf(“%d,%,%d”,&m,&n,&t);//读取需要建立的矩阵的基本信息 

H=malloc(sizeof(MNode)); /*申请总头结点*/
H->row=m; H->col=n;

hd[0]=H;//将头节点记录在辅助数组中
 
for(i=1; i<s; i++)//分配s个头节点 
{ 
	p=malloc(sizeof(MNode)); /*申请第i 个头结点*/
	p->row=0; p->col=0;
	p->right=p;
	p->down=p;//初始化行和列的循环链表 
	hd[i]=p;//将头节点地址记录在辅助数组中 
	hd[i-1]->v_next.next=p;//并与上一个头节点的next连接起来 
}

hd[S]->v_next.next=H;//最终，头节点们形成一个循环链表
/*
这里需要注意的是，并不是像结构图上的那样，有行头节点和列头节点，其实行和列公用头节点 
*/ 

for(k=1;k<=t;k++)//输入k个非零的三元组，并插入十字链表中 
{
	scanf (“%d,%d,%d”,&i,&j,&v);
	p=malloc(sizeof(MNode));
	p->row=i ; p->col=j; p->v_next.v=v
/*以下是将*p 插入到第i 行链表中去，且按列号有序*/
	q=hd[i];//从辅助数组中得到所在行的行头节点地址 
	while(q->right!=hd[i] && (q->right->col)<j )//遍历查找，找到要插入的列，终止条件是，起初就是空链，或者找到了比自己更大的第一个列号 
		q=q->right;
	p->right=q->right; //插入操作，拆链 
	q->right=p;//连接 
/*以下是将*p 插入到第j 行链表中去，且按行号有序*/
	q=hd[i];
	while ( q->down!=hd[j] && (q->down->row)<i ) /*按行号找位置*/
	q=q->down;
	p-> down =q-> down; /*插入*/
	q-> down =p;
} /*for k*/
return H;
} /* CreatMLink */