稀疏矩阵的存储和快速转置

稀疏矩阵的存储和快速转置

一、目的：

掌握稀疏矩阵的三元组表的压缩存储方法，并且利用该存储方法实现稀疏矩阵的转置与快速转置。

二、要求：

输入一个稀疏矩阵A，由程序将其转换成三元组表存储；转置后的三元组表，由程序将其转换成矩阵形式后输出。

三、实验内容

1、设计程序。

2、调试程序，并设计输入数据。

四、实验报告要求

写出程序和实验结果。

#define MAXSIZE 50

#include "Stdio.h"

#define X 10

#define Y 10

typedef struct{

int i,j;

int e;

}Triple;

typedef struct

{ Triple data[MAXSIZE+1];

int mu,nu,tu;

}TSMatrix;

int num[X],cpot[Y];

TSMatrix createsmatrix()

{

TSMatrix M;

int i,j,k=1,x,m,n;

int a[X][Y];

printf("请输入矩阵的行数mu:\n");

scanf("%d",&m);

printf("请输入矩阵的列数nu:\n");

scanf("%d",&n);

printf("请输入矩阵的各个元素:\n");

for(i=0;i<m;i++)

for(j=0;j<n;j++)

scanf("%d",&a[i][j]);

for(i=0;i<m;i++)

{ for(j=0;j<n;j++)

{ x=a[i][j];

if(x)

{M.data[k].i=i+1;

M.data[k].j=j+1;

M.data[k].e=x;

k++;}

} }

M.mu=m;

M.nu=n;

M.tu=k-1;

return(M);

}

TSMatrix fasttransposesmatrix(TSMatrix M)

{ TSMatrix T;

int col,t,p,q;

T.mu=M.nu;

T.nu=M.mu;

T.tu=M.tu;

if(T.tu)

{

for(col=1;col<=M.nu;++col)

num[col]=0;

for(t=1;t<=M.tu;++t)

++num[M.data[t].j];

cpot[1]=1;

for(col=2;col<=M.nu;++col)

cpot[col]=cpot[col-1]+num[col-1];

for(p=1;p<=M.tu;++p)

{

col=M.data[p].j;

q=cpot[col];

T.data[q].i=M.data[p].j;

T.data[q].j=M.data[p].i;

T.data[q].e=M.data[p].e;

++cpot[col];

}

}

return(T);

}

void print(TSMatrix T)

{ int b[X][Y];

int i,j,t,m=T.mu,n=T.nu;

for(i=0;i<m;i++)

for(j=0;j<n;j++)

b[i][j]=0;

for(t=1;t<=T.tu;t++)

b[(T.data[t].i)-1][(T.data[t].j)-1]=T.data[t].e;

printf("转置矩阵是:");

for(i=0;i<m;i++)

{ printf("\n");

for(j=0;j<n;j++)

printf("%5d",b[i][j]);

}

}

void main()

{ int i,j;

TSMatrix M,T;

M=createsmatrix();

T=fasttransposesmatrix(M);

print(T);

printf("\n");

}