《数据结构复习》三元组表求转置矩阵

• 对M的列进行转置
  先对列进行循环，再对矩阵非零元进行，列相等时，则存入转置T

---

• 快速转置
  引入变量num[col]、cpot[col]，表示每列非零元的个数以及第一个非零元的位置
  取出非零元的列和第一个的位置，存入T的(cpot[col]，j)，cpot[col]++
  读取下一个非零元

//数组三元组表存储
//有序的双下标法
//ch1_7
#include<stdio.h>
#include<stdlib.h>
#include<stdarg.h>
#define OK 1
#define ERROR 0
#define maxsize 10000
typedef int datatype;
typedef int status;
typedef struct
{
	int i, j;
	datatype e;
}Triple;
typedef struct
{
	Triple data[maxsize];
	int  m, n, t;
}TSMatrix;

//按M的列进行转置
status tranmatrix(TSMatrix M, TSMatrix* T) {
	T->m = M.n, T->n = M.m, T->t = M.t;
	if (T->t) {
		int q = 1;
		for(int col=1;col<M.n;col++)//对列进行循环
			for (int p = 1; p < M.t; 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++;
				}
			}
	}
	return OK;
}

//快速转置
status fastmatrix(TSMatrix M, TSMatrix* T) {
	T->m = M.n, T->n = M.m, T->t = M.t;
	int num[20];//想用M的列，但是不允许
	int cpot[20];
	if (T->t) {
		for (int col = 1; col < M.n; col++)	//num数组置0
			num[col] = 0;
		for (int t = 1; t < M.t; t++)	//求每列的非零元个数
			num[M.data[t].j]++;//非零元列num+1
		cpot[1] = 1;//第一列第一个非零元存在第一位
		for (int col = 2; col < M.n; col++)
			cpot[col] = cpot[col - 1] + num[col-1];//前一列的个数加上第一个的位置

		//求出两辅助变量的值
		for (int p = 1; p <= M.t; p++) {//非零元循环
			int col = M.data[p].j;//取出非零元的列
			int q = cpot[col];//所在列的第一位存储位置
			T->data[q].i = M.data[p].j;//即存到cpot[col]的位置
			T->data[q].j = M.data[p].i;
			T->data[q].e = M.data[p].e;
			cpot[col]++;
		}
	}
	return OK;
}