c语言版数据结构（奇迹冬瓜）-数组和广义表（稀疏矩阵的转置算法一）

//稀疏矩阵
/*
  |0   12   9   0   0   0   0|     |0   0  -3   0   0  15| 
  |0   0    0   0   0   0   0|     |12  0   0   0  18   0|
M=|-3  0    0   0   0  14   0|   T=|9   0   0  24   0   0|
  |0   0    24  0   0   0   0|     |0   0   0   0   0  -7|
  |0   18   0   0   0   0   0|     |0   0   0   0   0   0| 
  |15  0    0   -7  0   0   0|     |0   0  14   0   0   0|
                                   |0   0   0   0   0   0|
M的转置矩阵为T
*/
#include<stdio.h>
#include<stdlib.h>

#define MAXSIZE 12500   //假设非零元个数的最大值为12500
#define OK 1
#define OVERFLOW -2

typedef struct 
{
	int i,j;           //非零元行的行下标和列下标
	int e;              //非零元的值
}Triple;
typedef struct 
{
	Triple data[MAXSIZE+1];          //非零元三元组表,data[0]未用
	int mu,nu,tu;                    //矩阵的行数,列数和非零个数
}TSMatrix;
typedef int Status;

Status TransposeSMatrix(TSMatrix M,TSMatrix *T)
{
	int q,p,col;
	T->mu=M.mu;
	T->nu=M.nu;
	T->tu=M.tu;
	if(T->tu)
	{
		q=1;
		for(col=1;col<=M.mu;++col)
		{
			for (p=1;p<=M.nu;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;
				}
			}
		}
	}
	return OK;
}

Status FastTransposeSMatrix(TSMatrix M,TSMatrix *T)
{
	T->mu=M.mu;
	T->nu=M.nu;
	T->tu=M.tu;
	int col,*num,t,*cpot,p,q;
	if(!(num=(int *)malloc((M.nu+1)*sizeof(int))))
	{
		exit(OVERFLOW);
	}
	if (!(cpot=(int *)malloc((M.nu+1)*sizeof(int))))
	{
		exit(OVERFLOW);
	}
	if (T->tu)
	{
		for (col = 0; col <= M.nu; col++)
		{
			num[col]=0;
		}
		for (t = 1; t <= M.tu; t++)
		{
			++num[M.data[t].j];
		}
		cpot[1]=1;
		for ( col = 2; col <= M.mu; col++)
		{
			cpot[col]=cpot[col-1]+num[col-1];
		}
		for (p = 1; p <= M.tu ; ++p)
		{
			col=M.data[p].j;
			q=cpot[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;
			++cpot[col];
		}
	}
}

void main()
{
	TSMatrix M,T;
	int i=1;
	M.data[1].i=1;
	M.data[1].j=2;
	M.data[1].e=12;
	M.data[2].i=1;
	M.data[2].j=3;
	M.data[2].e=9;
	M.data[3].i=3;
	M.data[3].j=1;
	M.data[3].e=-3;
	M.data[4].i=3;
	M.data[4].j=6;
	M.data[4].e=14;
	M.mu=6;
	M.nu=7;
	M.tu=4;
	FastTransposeSMatrix(M,&T);
	for(i= 1; i<=T.tu;i++)
	{
		printf("%d %d %d\n",T.data[i].i,T.data[i].j,T.data[i].e);
	}
	getchar();
}