稀疏矩阵快速转置


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

}