数据结构_稀疏矩阵转置（trans_sparse_matrix）

最新推荐文章于 2023-05-08 11:24:21 发布

原创最新推荐文章于 2023-05-08 11:24:21 发布 · 1.8k 阅读

3 ·

CC 4.0 BY-SA版权

本文介绍了一种求解稀疏矩阵转置的有效方法，并提供了两种实现方式：一种是按列优先顺序进行存储的方式；另一种是快速转置算法。通过这两种方法，可以有效地将稀疏矩阵转换为其转置形式。

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

经过几个小时的奋战，AC～～~ ^_^

数据结构_稀疏矩阵转置_方法一_用M的列优先表示（trans_sparse_matrix） - Essence - ACM_Of_Essence

下面是用M的列优先表示，

Description
/************************************/
/*求转置矩阵*/
/*问题描述：已知一个稀疏矩阵的三元组表，*/
/*求该矩阵转置矩阵的三元组表*/
/************************************/

Sample Output

code：

//author:essence_of_acmer;
//time:2011.8.7
//theme:trans_sparse_matrix
//用M的列优先表示

#include <iostream>
using namespace std;
#define M 20

typedef struct node
{
int i,j;
int v;
} JD ;

int trans_Sparse_matrix(JD ma[],JD mb[])  //way 1：按照M的列序存储
{
    if(ma[0].v==0)
       return (0);
       int n=ma[0].j;
       int t=ma[0].v;
       mb[0].i=n;
       mb[0].j=ma[0].i;
       mb[0].v=t;
       int k=1;
       for(int col=1;col<=n;col++)
           for(int p=1;p<=t;p++)
               if (ma[p].j==col)
               {
                   mb[k].i=ma[p].j;
                   mb[k].j=ma[p].i;
                   mb[k].v=ma[p].v;
                   k++;
               }
return (1);
}

int fast_transpos(JD ma[],JD mb[])  //way 2:快速转置
{  int n,col,p,k,t;
   int num[M],cpot[M];
   n=ma[0].j;
   t=ma[0].v;
   mb[0].i=n;  mb[0].j=ma[0].i;  mb[0].v=t;
   if(t<=0)
      return(0);
   for(col=0;col<=n;col++)
       num[col]=0;   //num[col]：表示矩阵M中第col列中非零元个数

for(p=1;p<=t;p++)//优秀
   {   k=ma[p].j;
       num[k]++;// 统计矩阵M中第col列中非零元个数
   }

   cpot[0]=0; cpot[1]=1;
   for(col=2;col<=n;col++)
  cpot[col]=cpot[col-1]+num[col-1];//统计M中第col列第一个非零元在mb中位置

   for(p=1;p<=t;p++)
   {  col=ma[p].j;  //ma[]中第col列
      k=cpot[col];  //第col列第一个非零元在mb中位置
      mb[k].i=ma[p].j;
      mb[k].j=ma[p].i;
      mb[k].v=ma[p].v;
      cpot[col]++;  //优秀     //经典代码段
   }
   return(1);
}

int main()
{
int a[6][7]={0,12,9,0,0,0,0,
                 0,0,0,0,0,0,0,
                -3,0,0,0,0,14,0,
                 0,0,24,0,0,0,0,
                 0,18,0,0,0,0,0,
                 15,0,0,-7,0,0,0};
JD ma[M],mb[M];
int k=1;

for(int i=0;i<6;i++)  //初始化
  for(int j=0;j<7;j++)
    if(a[i][j]) {ma[k].i=i+1;ma[k].j=j+1;ma[k].v=a[i][j];k++;}
    ma[0].i=6;ma[0].j=7;ma[0].v=k-1;

//trans_Sparse_matrix(ma,mb);  //way 1;
fast_transpos(ma,mb);  //way 2;

for(int i=0;i<k;i++)  //打印原始
  cout<<ma[i].i<<"  ,"<<ma[i].j<<"  ,"<<ma[i].v<<endl;

cout<<endl<<endl;

for(int i=0;i<k;i++)  //打印改变后的
  cout<<mb[i].i<<"  ,"<<mb[i].j<<"  ,"<<mb[i].v<<endl;
return (0);
}