#include <stdio.h>
#include <stdlib.h>
typedef struct{
int data[100][100]; //存储矩阵
int m,n; //行,列数
}matrix;
typedef int spmatrix[100][3]; //三元数组别名
spmatrix A,B; //创建两个三元数组,A为转置前,B为转置后
int pos[1000],num[1000];
int GetThree(matrix *s,int len1,int len2) // 将非零元素以行,列,数的值的形式存入三元组
{
int i,j;
int sum=1;
for(i=1;i<=s->m;i++)
{
for(j=1;j<=s->n;j++)
if(s->data[i][j])
{
A[sum][0]=i;
A[sum][1]=j;
A[sum++][2]=s->data[i][j];
}
}
A[0][0]=len1;
A[0][1]=len2;
A[0][2]=sum-1;
return sum;
}
void Getposnum(int len,int col)
{
int i,c;
memset(num,0,sizeof(num));
for(i=0;i<len;i++)
{
c=A[i][1];
num[c]++; //对每列含有的元素计数
}
pos[1]=1;
for(i=2;i<=col;i++)
pos[i]=pos[i-1]+num[i-1]; //计算每列非0元素开始下标
}
void Inverse(int len)
{
int i,j,c;
B[0][0]=A[0][1];
B[0][1]=A[0][0];
B[0][2]=A[0][2];
for(i=1;i<len;i++)
{
c=A[i][1];
B[pos[c]][0]=A[i][1]; //开始转置,pos[c]为存入位置下标
B[pos[c]][1]=A[i][0];
B[pos[c]][2]=A[i][2];
pos[c]++; //每段相对位置加1
}
}
int main()
{
int i,j,res,col,k;
matrix s1;
printf("请输入稀疏矩阵的行数和列数:\n");
scanf("%d%d",&res,&col); //输入行,列数
s1.m=res;
s1.n=col;
printf("请输入矩阵\n");
for(i=1;i<=res;i++)
{
for(j=1;j<=col;j++)
scanf("%d",&s1.data[i][j]); //读取矩阵元素
}
int len=GetThree(&s1,res,col); //矩阵非0元素个数
Getposnum(len,col);
Inverse(len);
printf("输出矩阵的转置:\n"); //输出矩阵的转置
for(i=1;i<=col;i++)
{
for(j=1;j<=res;j++)
{
for(k=1;k<len;k++)
{
if(B[k][0]==i&&B[k][1]==j)
{
printf("%-3d",B[k][2]);
break;
}
}
if(k==len)printf("0 "); //表中未查询到则输出0
}
printf("\n");
}
return 0;
}