三元组矩阵转置

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

#define OK 1
#define ERROR 0
#define MAXSIZE 256

typedef int ElemType;
typedef int Status;

typedef struct
{
	int i, j;    //元素在矩阵中的行号，列号
	ElemType v;    //元素值
}Triple;

typedef struct
{
	Triple arr[MAXSIZE];    //非零元素三元组
	int Rows, Cols, Nums;    //矩阵的行数，列数，非零元个数
}SqSMatrix;    //稀疏矩阵的三元组顺序表示

//三元组稀疏矩阵的转置方法——1.简单方法、2.列序遍历法、3.快速转置法

/*1.简单方法：
    把三元组表中行和列的内容互换，然后按照行号对新的三元组表中的各三元
组从小到大进行排序，就可以得到转置后的三元组表，排序算法采用经典的排序
算法，则时间复杂度为O(Nums*Nums)
*/

//2.列序遍历的转置算法
Status TransposeSMatrix(SqSMatrix A, SqSMatrix& B)
{
	//将稀疏矩阵A转置，结果放在稀疏矩阵B中
	B.Rows = A.Cols;    //B的行数等于A的列数
	B.Cols = A.Rows;    //B的列数等于A的行数
	B.Nums = A.Nums;    //B的非零元个数等于A的非零元个数

	if (A.Nums > 0)
	{
		int q = 0;
		int p = 0;
		for (int k = 0; k < A.Cols; k++)
		{
			for (p = 0; p < A.Nums; p++)
			{
				if (A.arr[p].j == k)
				{
					//若第p个三元组中元素的列号为k
					B.arr[q].i = A.arr[p].j;    //新三元组中的行号
					B.arr[q].j = A.arr[p].i;    //新三元组中的列号
					B.arr[q].v = A.arr[p].v;    //新三元组中的值
					q++;    //B.arr的当前位置增加1
				}
			}
		}
	}
	return OK;
}

//3.稀疏矩阵快速转置算法
Status FastTransposeSMatrix(SqSMatrix A, SqSMatrix& B)
{
	B.Rows = A.Cols;
	B.Cols = A.Rows;
	B.Nums = A.Nums;

	int rowStart[MAXSIZE];
	int rowNum[MAXSIZE];
	int k = 0;
	int p = 0;
	if (A.Nums > 0)
	{
		//第一步：统计A中每列非零元素个数
		for (k = 0; k < A.Cols; k++)    rowNum[k] = 0;    //rowNums数组的初始化
		for (p = 0; p < A.Nums; p++)    rowNum[A.arr[p].j]++;    //统计A中每列（即B中每行）非零元素个数
		//第二步：计算B中每行三元组的起始位置
		rowStart[0] = 0;
		for (k = 1; k < A.Cols; k++)
			rowStart[k] = rowStart[k - 1] + rowNum[k - 1];
		for (p = 0; p < A.Nums; p++)
		{
			int q = rowStart[A.arr[p].j];
			B.arr[q].i = A.arr[p].j;
			B.arr[q].j = A.arr[p].i;
			B.arr[q].v = A.arr[p].v;
			rowStart[A.arr[p].j]++;
		}
	}
	return OK;
}