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();
}