三元组十字链表下稀疏矩阵的加、转、乘实现

#include<stdio.h>
#include<stdlib.h>
#include<iostream.h>
#include<process.h>
#define OK 1
#define ERROR 0
#define OVERFLOW -2
typedef int ElemType;

//------------------------------------------------------------------------------------------------

struct OLNode
{
 int i,j;                                           //非零元所在行、列
 ElemType e;                                        //非零元值
 OLNode *right,*down;                               //定义向右域和向下域
};

typedef OLNode *OLink;

struct CrossList
{
 OLink *rhead,*chead;                               //行、列表头的头节点
 int mu,nu,tu;                                      //矩阵的行、列和非零元个数
};

//------------------------------------------------------------------------------------------------

void Create(CrossList &M)                             //矩阵的创建
{
 int i,j,k,m,n,t;
 ElemType e;                                       //用于存放非零元值
 OLNode *p,*q;
 printf("Please input the rows 、columns and the values of the matrix：");
 scanf("%d%d%d",&m,&n,&t);                         //输入稀疏矩阵的行、列、非零元
 M.mu=m;                                           //把矩阵的总行数存入M.mu
 M.nu=n;                                           //把矩阵的总列数存入M.nu
 M.tu=t;                                           //把矩阵的非零元个数存入M.tu
 M.rhead=(OLink*)malloc((m+1)*sizeof(OLink));      //给对象M的行分配m+1个结点大小的空间用行头指针M.rhead指向第一个行结点
 if(!M.rhead)                                      //假如头结点分配成功
  exit(OVERFLOW);                               //退出程序，程序返回给操作系统的返回码是OVERFLOW
 M.chead=(OLink*)malloc((n+1)*sizeof(OLink));      //给对象M的列分配n+1个结点大小的空间用列头指针M.chead指向第一个列结点
 if(!M.chead)                                      //假如头结点分配成功
  exit(OVERFLOW);                               //退出程序，程序返回给操作系统的返回码是OVERFLOW
 for(k=1;k<=m;k++)                                 //初始化行头指针
  M.rhead[k]=NULL;
 for(k=1;k<=n;k++)                                 //初始化列头指针
  M.chead[k]=NULL;
 printf("Please input %d data:/n",M.tu);
 for(k=0;k<t;k++)
 {
  scanf("%d%d%d",&i,&j,&e);                     //输入t个非零元的信息
  if(i>m||j>n)                                  //假如输入的元素不在矩阵中
  {
            printf("Please input again:/n");
   exit(OVERFLOW);                           //退出程序，程序返回给操作系统的返回码是OVERFLOW
  }
  else                                          //假如输入的元素在矩阵中
  {
   p=(OLNode*)malloc(sizeof(OLNode));        //为p分配一个结点大小的空间
   if(!p)                                    //假如p不存在
    exit(OVERFLOW);                       //退出程序，程序返回给操作系统的返回码是OVERFLOW
   p->i=i;                                   //把非零元的行值赋给p->i
   p->j=j;                                   //把非零元的列值赋给p->j
   p->e=e;                                   //把非零元的值赋给p->e
   if(M.rhead[i]==NULL||M.rhead[i]->j>j)     //如果第i行的第一个结点为空或者要插入的非零元的列值小于第i行的第一个结点的列值
   {
    //p插入该行第一节点处
    p->right=M.rhead[i];                  //p的向右域指向没有插入前的第i行的第一个结点
    M.rhead[i]=p;                         //p指向没有插入前的第i行的第一个结点
   }
   else
   {
    //否则寻找行表插入位置
    for(q=M.rhead[i];q->right&&q->right->j<j;q=q->right);//这里的;表示对这个循环不执行任何操作
    p->right=q->right;                    //p的右向域指向q的右向域
    q->right=p;                           //p指向q的右向域
   }                                         //完成行插入
   if(M.chead[j]==NULL||M.chead[j]->i>i)     //假如第一个列结点为空或者第j列的第一个结点的行值大于要插入的行值
   {
    //p插入该列第一节点处
    p->down=M.chead[j];                   //p的向下域指向没有插入前的第j列的第一个结点
    M.chead[j]=p;                         //p指向没有插入前的第j列的第一个结点
   }                                         //完成列插入
   else
   {
    //否则寻找列表插入位置
    for(q=M.chead[j];q->down&&q->down->i<i;q=q->down);//这里的;表示对这个循环不执行任何操作
    p->down=q