稀疏矩阵概念
个人理解:一个非零元很少且分布没有规律的矩阵
如何存储?
这里暂且只介绍三元组顺序表的方法
我们使用一个三元组(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);
}