用三元组存储稀疏矩阵并实现转置

基本概念

在学习线性代数的时候，经常用到矩阵。在C语言中，表示矩阵的最直观形式就是二维数组。然而在实际应用中，很多高阶矩阵中的非零元素非常少，这个时候如果继续使用二维数组存储，那么就会浪费很多存储空间。

在数据结构中，我们用三元组存储稀疏矩阵。三元组定义为(i,v,j)，这三个值一次表示矩阵的行、列、值。

有了基本的概念之后，就可以定义数据结构了

定义一个结构体，来表示三元组的基本属性

typedef struct
{
	int row, col;
	int e;
}Triple;

然后再定义一个存储容器，用来存放三元组的

为了简单起见，我们用数组来实现，并定义最大存储单元MAXSIZE为100

typedef struct
{
	Triple data[MAXSIZE];
	Int m,n,len;
}TSMatrix;
//(TSMatrix表示 Triple Sparse Matrix)

实现矩阵的转置
         实现用三元组表示的矩阵的转置，可以直接把行列互换，然后再执行按行序为主的排序过程。为了避免重新排序引起的元素移动，可以采用列序递增转置法。

         具体做法，就是遍历列的下表值，从列数低的值到列数高的值，依次添加到缓存三元组中。很显然，这是一个双重for循环结构，内层循环实现遍历整个表，寻找合适的列。外层循环，则记录要寻找的列数。

//实现转置 
void TransposeTSMatrix(TSMatrix A, TSMatrix* B)
{
	int i,j,k;
	B->m = A.n;
	B->n = A.m;
	B->len = A.len;
	j=0;
	
	for( k=0; k<A.len; ++k)
	{
		for( i=0; i<A.len; ++i)
		{
			if(A.data[i].col == k)
			{
				B->data[j].row = A.data[i].col;
				B->data[j].col = A.data[i].row;
				B->data[j].e = A.data[i].e;
				++j;
			}
		}
	}
	
}

有了上面的基础，就可以写一个带有测试驱动的函数

完整代码

#include <stdio.h>
#define MAXSIZE 100
//三元组的定义 
typedef struct
{
	int row, col;//表示行列 
	int e;		//表示值 
}Triple;

//三元组容器的定义 
typedef struct
{
	Triple data[MAXSIZE];
	int m,n,len;
}TSMatrix;

//实现转置 
void TransposeTSMatrix(TSMatrix A, TSMatrix* B)
{
	int i,j,k;
	B->m = A.n;
	B->n = A.m;
	B->len = A.len;
	j=0;
	
	for( k=0; k<A.len; ++k)
	{
		for( i=0; i<A.len; ++i)
		{
			if(A.data[i].col == k)
			{
				B->data[j].row = A.data[i].col;
				B->data[j].col = A.data[i].row;
				B->data[j].e = A.data[i].e;
				++j;
			}
		}
	}
	
}

//测试驱动函数 
int main()
{
	//将输入重定向到根目录下的data.txt 
	freopen("data.txt", "r", stdin);
	TSMatrix A,B;
	int i,j,e; 
	int k=0;
	printf("请输入三元组：");
	while(scanf("%d%d%d", &i, &j, &e)!=EOF)
	{
		A.data[k].row = i-1;
		A.data[k].col = j-1;
		A.data[k].e = e;
		A.len = ++k;
	}
	printf("\n原始三元组为:\n");
	for(i=0; i<A.len; ++i )
	{
		printf("%3d%3d%3d\n", A.data[i].row+1, A.data[i].col+1, A.data[i].e);
	}
	
	printf("\n转置后:\n");
	TransposeTSMatrix(A, &B);
	for(i=0; i<B.len; ++i )
	{
		printf("%3d%3d%3d\n", B.data[i].row+1, B.data[i].col+1, B.data[i].e);
	}
	return 0;
}

程序截图