稀疏矩阵快速转置



#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;
}