稀疏矩阵的计算概论

#include

#include

#include

typedef int ElemType;// 稀疏矩阵的十字链表存储表示

typedef struct OLNode

{

int i,j; // 该非零元的行和列下标

ElemType e; // 非零元素值

struct OLNode *right,*down; // 该非零元所在行表和列表的后继链域}OLNode, *OLink;

typedef struct// 行和列链表头指针向量基址,由CreatSMatrix_OL()分配{

OLink *rhead, *chead;

int mu, nu, tu; // 稀疏矩阵的行数、列数和非零元个数

}CrossList;

// 初始化M(CrossList类型的变量必须初始化,否则创建、复制矩阵将出错) 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

{

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++)