//TSMatrix.cpp
#include"predefined.h"
#include"TSMatrix.h"
Status TransposeSMatrix(TSMatrix M,TSMatrix *T)
//算法5.1:采用三元组表存储表示,求稀疏矩阵M的转置矩阵T
{
int q,p,col;
(*T).mu=M.mu;
(*T).nu=M.nu;
(*T).tu=M.tu;
if(M.tu)
{
q=1;
for(col=1;col<=M.nu;col++)
for(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++;
}
}
return OK;
}
Status FastTransposeSMatrix(TSMatrix M,TSMatrix *T)
//算法5.2:采用三元顺序表存储表示,求稀疏矩阵M的转置矩阵T
{
int num[10],cpot[10],t,p,col,q;
(*T).mu=M.mu;
(*T).nu=M.nu;
(*T).tu=M.tu;
if(M.tu)
{
for(col=1;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.nu;col++)
cpot[col]=num[col-1]+cpot[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]++;
}
}
return OK;
}
Status CreateSMatrix(TSMatrix *M)
//创建稀疏矩阵
{
int i,m,n;
ElemType e;
printf("请输入矩阵的行、列、非零元个数:");
scanf("%d,%d,%d",&(*M).mu,&(*M).nu,&(*M).tu);
for(i=1;i<=(*M).tu;i++)
{
printf("请按行序顺序输入第%d个非零元素所在的行(1~%d),列(1~%d),元素值:",i,(*M).mu,(*M).nu);
scanf("%d,%d,%d",&m,&n,&e);
(*M).data[i].i=m;
(*M).data[i].j=n;
(*M).data[i].e=e;
}
return OK;
}
void PrintSMatrix(TSMatrix M)
//打印稀疏矩阵
{
int i;
printf("%d行%d列%d个非零元素。\n",M.mu,M.nu,M.tu);
printf("行 列 元素值\n");
for(i=1;i<=M.tu;i++)
printf("%2d%4d%8d\n",M.data[i].i,M.data[i].j,M.data[i].e);
}
//TSMatrix.h
#define MAXSIZE 12500
typedef int ElemType;
typedef struct
{
int i,j;
ElemType e;
}Triple;
typedef struct
{
Triple data[MAXSIZE+1];//date[0]未用
int mu,nu,tu;
}TSMatrix;
Status TransposeSMatrix(TSMatrix M,TSMatrix *T);
Status CreateSMatrix(TSMatrix *M);
void PrintSMatrix(TSMatrix M);
Status FastTransposeSMatrix(TSMatrix M,TSMatrix *T);