稀疏矩阵的压缩存储及其转置

稀疏矩阵概念

个人理解：一个非零元很少且分布没有规律的矩阵

如何存储？

这里暂且只介绍三元组顺序表的方法

我们使用一个三元组（i，j，k)来确定矩阵中的一个非零元，i表示行，j表示列，k即是第i行第j列的非零元。所以，稀疏矩阵可以用一个三元组表来表示。例如下图：

 我们根据创建的三元组表，打印对应的元素，其余位置打印0，这样一个稀疏矩阵就完成了

稀疏矩阵的转置

矩阵的转置即是将行和列对调，对应到三元组表就是将i和j相互调换，重新排列三元组的次序。

 如何处理呢？

一种方法就是对三元组表1的j进行从1开始的遍历扫描，每次找到一个就转置到表2中。

void TransposeSmatrix(TsMatrix& m, TsMatrix& t)
{//m为待转置的三元组表，t为转置后的
	t.mu = m.nu;//行列数交换
	t.nu = m.mu;
	t.tu = m.tu;
	if (t.tu)
	{
		int q = 1;//t表的下标
		for (int col = 1; col <= m.nu; col++)//对列进行扫描
		{
			for (int p = 1; p <=m.tu; p++)//对m表遍历
			{
				if (m.data[p].j == 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;
					q++;
				}
			}
		}
	}
}

另一种方法是快速转置，我没认真学，以后有机会补上

完整代码

#include<stdio.h>
#include<stdlib.h>

#define Max_size 100
typedef struct Triple  //三元组
{
	int i;//非零元的行下标
	int j;//非零元的列下标
	int e;//该位置对应的元素
}Triple;
typedef struct TsMatrix
{
	Triple data[Max_size + 1];//三元组表
	int mu, nu, tu;//矩阵行数、列数、非零元个数
}TsMatrix;

void InitTsMatrix(TsMatrix&p)
{
	p.mu = rand() % 51 + 6;
	p.nu= rand() % 51 + 8;
	p.tu= rand() % 51 + 6;
	if (p.tu >= p.mu *p. nu)
	{
		p.tu = 6;//对于随机产生的tu如果个数超过了矩阵总元素个数，我们直接初始化为6个
	}
	for (int i = 0; i < p.tu;i++)
	{
		printf("请输入三元组：\n");
		scanf("%d %d %d", &(p.data[i+1].i), &(p.data[i + 1].j), &(p.data[i + 1].e));
	}
}

void print(TsMatrix& p)
{
	int flag = 1;
	for (int row = 1; row <=p.mu; row++)
	{
		for (int col = 1; col <=p.nu; col++)
		{
			if (p.data[flag].i == row && p.data[flag].j == col)
			{
				printf("%3d", p.data[flag].e);
				flag++;
			}
			else {
				printf("%3d", 0);
			}
		}
		printf("\n");
	}
}
void TransposeSmatrix(TsMatrix& m, TsMatrix& t)
{
	t.mu = m.nu;
	t.nu = m.mu;
	t.tu = m.tu;
	if (t.tu)
	{
		int q = 1;
		for (int col = 1; col <= m.nu; col++)
		{
			for (int p = 1; p <=m.tu; p++)
			{
				if (m.data[p].j == 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;
					q++;
				}
			}
		}
	}
}
int main()
{
	TsMatrix a;
	TsMatrix b;
	InitTsMatrix(a);
	print(a);
	TransposeSmatrix(a, b);
	print(b);
}