数据结构实验之数组三：快速转置

Problem Description

 

转置运算是一种最简单的矩阵运算，对于一个m*n的矩阵M( 1 = < m < = 10000,1 = < n < = 10000 )，它的转置矩阵T是一个n*m的矩阵，且T( i , j )=M( j , i )。显然，一个稀疏矩阵的转置仍然是稀疏矩阵。你的任务是对给定一个m*n的稀疏矩阵( m , n < = 10000 )，求该矩阵的转置矩阵并输出。矩阵M和转置后的矩阵T如下图示例所示。
   
   稀疏矩阵M                             稀疏矩阵T

 

Input

连续输入多组数据，每组数据的第一行是三个整数mu, nu, tu(tu <= 50)，分别表示稀疏矩阵的行数、列数和矩阵中非零元素的个数，随后tu行输入稀疏矩阵的非零元素所在的行、列值和非零元素的值，同一行数据之间用空格间隔。(矩阵以行序为主序)

Output

输出转置后的稀疏矩阵的三元组顺序表表示。

Sample Input

3 5 5
1 2 14
1 5 -5
2 2 -7
3 1 36
3 4 28

Sample Output

1 3 36
2 1 14
2 2 -7
4 3 28
5 1 -5

#include <stdio.h>
#include <stdlib.h>
#include <string.h>
#define MAXSIZE 100000
typedef struct
{
    int h;
    int l;//非零元素的行标和列标
    int x;
}Triple;//三元组
typedef struct
{
    Triple data[MAXSIZE+1];//非零元素的三元组表。data[0]未用
    int ha,lie,num;//矩阵的行数，列数的非零元素的个数
}TSMatrix;
int a[MAXSIZE],b[MAXSIZE];
void FastTransposeMatrix(TSMatrix *S,TSMatrix *T)
{
    int i,k;
    T->ha=S->lie;T->lie=S->ha;T->num=S->num;
    if(T->num!=0)
    {
        memset(a,0,MAXSIZE);
        for(i=1;i<=S->num;i++) a[S->data[i].l]++;//统计S中每一列所含非零元素的个数
        b[1]=1;//第一个元素的起始位置为1
        for(i=2;i<=S->lie;i++) b[i]=b[i-1]+a[i-1];
        //求T中每一行的起始元素的位置
        for(i=1;i<=S->num;i++)
        {
            int j=S->data[i].l;//记录元素的列标
            k=b[j];//记录该行的第一个元素的位置
            T->data[k].h=S->data[i].l;
            T->data[k].l=S->data[i].h;
            T->data[k].x=S->data[i].x;//交换行列号和元素数据
            b[j]++;//该行元素位置后移一位，加1
        }//矩阵的转置
    }
    for(i=1;i<=T->num;i++) printf("%d %d %d\n",T->data[i].h,T->data[i].l,T->data[i].x);
    //输出三元组
}
int main()
{
    TSMatrix M,T;
    int i,x,y,z;
    while(scanf("%d%d%d",&x,&y,&z)!=EOF)
    {
        M.ha=x;M.lie=y;M.num=z;
        for(i=1;i<=M.num;i++)
        {
            scanf("%d%d%d",&M.data[i].h,&M.data[i].l,&M.data[i].x);
        }
        FastTransposeMatrix(&M,&T);
    }
    return 0;
}