稀疏矩阵的快速转置

#include<iostream>
#include<cstdio>
#define maxn 111
using namespace std;
int a[maxn][maxn],n,m,cnt=1;
typedef struct
{
int x,y,val;//x为横坐标，y为纵坐标,val为值
}matrix;
void Scan()//输入矩阵
{
for(int i=1;i<=n;i++)
for(int j=1;j<=m;j++)  scanf("%d",&a[i][j]);
}
void Print() //打印转置前的矩阵
{
for(int i=1;i<=n;i++)
{
for(int j=1;j<=m;j++)
printf("%-3d",a[i][j]);
printf("\n");
}
}
void tranTothree(matrix A[]) //将原矩阵转化为三元组表
{
for(int i=1;i<=n;i++)
for(int j=1;j<=m;j++)
if(a[i][j]) A[cnt].x=i,A[cnt].y=j,A[cnt].val=a[i][j],cnt++; //如果元素不为0，则写入三元组表
A[0].x=n,A[0].y=m,A[0].val=cnt-1; //三元组第一结构体存储矩阵的行，列，三元组内容的个数
}
void Printthree(matrix A[])//打印三元组表
{
for(int i=1;i<=A[0].val;i++)
printf("x=%d,y=%d,val=%d\n",A[i].x,A[i].y,A[i].val);
}
void trans(matrix A[],matrix B[]) //矩阵转置
{
int x[maxn],y[maxn];//x数组记录原矩阵每一列的个数（不含零）,y数组记录转置后的矩阵每一行的起始位置
fill(x,x+maxn,0);
fill(y,y+maxn,0);
y[0]=1;//从y[1]开始记录,y[0]空间不用，但也要加入总数
B[0].x=A[0].y,B[0].y=A[0].x,B[0].val=A[0].val; //将转置后的行列数赋给新的三元表
int nx=A[0].x,ny=A[0].y,ncnt=A[0].val,id;
for(int i=1;i<=ncnt;i++)  id=A[i].y,x[id]++; //找转置前每一列元素的个数，也就是转置后每一行的个数
for(int i=1;i<=m;i++)  y[i]=y[i-1]+x[i-1]; //扫一遍整个转置前整个列，记录转置后每一行的起始位置
for(int i=1;i<=ncnt;i++) //扫一遍原三元表，进行转置
{
id=A[i].y;
if(!y[id]) continue; //若个数为零，则不用进行存储
B[y[id]].x=A[i].y,B[y[id]].y=A[i].x,B[y[id]].val=A[i].val;
y[id]++;//起始位置加一
}
}
void Print_trans(matrix B[])//将转置后的三元表打印成矩阵
{
int cnt=1;
for(int i=1;i<=B[0].x;i++)
{
for(int j=1;j<=B[0].y;j++)
{
if(B[cnt].x==i&&B[cnt].y==j) printf("%-3d",B[cnt].val),cnt++;
else printf("%-3d",0);
}
printf("\n");
}
}
int main()
{
matrix A[maxn],B[maxn];
while(~scanf("%d%d",&n,&m))
{
Scan();
tranTothree(A);
//Printthree(A);
trans(A,B);
Printthree(B);
Print_trans(B);
}
return 0;
}