矩阵的快速转置

最新推荐文章于 2023-11-17 19:50:43 发布

原创最新推荐文章于 2023-11-17 19:50:43 发布 · 801 阅读

0 ·

CC 4.0 BY-SA版权

数据结构专栏收录该内容

12 篇文章

订阅专栏

本文介绍了一种使用三元组表示稀疏矩阵的方法，并实现了稀疏矩阵的快速转置算法。通过定义结构体来存储非零元素的位置和值，有效地减少了存储空间的需求。此外，还提供了一个示例程序，演示了如何输入一个稀疏矩阵并进行快速转置。

摘要生成于 C知道，由 DeepSeek-R1 满血版支持，前往体验 >

#include<stdio.h>
#include<malloc.h>
#define MAXSIZE 100
typedef struct
{
int i,j;
int e;
}Triple;
typedef struct{
Triple data[MAXSIZE+1];
int mu,nu,tu;
} Tsmatrix;
int a,b; /*定义全局变量数组的行数a和列数b*/ /*用数组创建三元组*/
Tsmatrix * creatarray(Tsmatrix *M) {
int m,n,p=1;
int c;
printf("please input the array A:\n");
for(m=1;m<=a;m++)
  for(n=1;n<=b;n++)
  {
   scanf("%d",&c);
   if(c!=0)
   {
    M->data[p].e=c;
    M->data[p].i=m;
    M->data[p].j=n;
    p++;
   }
  }
  M->tu=p;
  M->mu=a;
  M->nu=b;
  printf("原矩阵三元组表示为 :\n\n");
  for(m=1;m<M->tu;m++)
   printf("%3d%3d%3d\t\n",M->data[m].i,M->data[m].j,M->data[m].e); printf("\n"); return M;
}
/*三元组快速转置*/
Tsmatrix * fasttrans(Tsmatrix *M,Tsmatrix *T)
{
int p,col,q,t,m; int num[100];
int cpot[100];
T->mu=M->nu;
T->nu=M->mu;
T->tu=M->tu;
if(T->tu!=0)
{
  for(col=1;col<=M->nu;++col)
   num[col]=0;
  for(t=1;t<M->tu;++t)
   num[M->data[t].j]++; cpot[1]=1;
  for(col=2;col<=M->nu;++col)
   cpot[col]=cpot[col-1]+num[col-1];
  for(p=1;p<M->tu;++p)
  {
   col=M->data[p].j;
   q=cpot[col];
   T->data[q].i=M->data[p].j;
   T->data[q].j=M->data[p].i;
   T->data[q].e=M->data[p].e;
   ++cpot[col];
  }
}
printf("\n转置之后的三元组表示 :\n");
for(m=1;m<T->tu;m++)
  printf("%3d%3d%3d\t\n",T->data[m].i,T->data[m].j,T->data[m].e);
printf("\n");
return T;
}
/*输出三元组函数*/
void print(Tsmatrix *T,int x,int y)
{
int m,n,p=1;
int d;
for(m=1;m<=x;m++)
{
  printf("\n");
  for(n=1;n<=y;n++)
  {
   if(T->data[p].i==m&&T->data[p].j==n)
   { d=T->data[p].e;
   p++;
   }
   else d=0;
   printf("%6d",d);
  }
}
printf("\n");
}
void main()
{
Tsmatrix *M,*T;
M=(Tsmatrix *)malloc(sizeof(Tsmatrix));
T=(Tsmatrix *)malloc(sizeof(Tsmatrix));
printf("请输入行数和列数:\n");
scanf("%d %d",&a,&b);
/*输入行列数*/
M=creatarray(M);
printf("创建的矩阵为:\n");
print(M,a,b);
T=fasttrans(M,T);
printf("转置后的矩阵为:\n");
print(T,b,a);

}