稀疏矩阵的转置

//稀疏矩阵的转置
//使用了两个辅助数组
//num  : num[j] : 矩阵A中第j列非零元素的个数
//cpot : cpot[j]: 矩阵A中第j列的第一个非零元素在其转置矩阵B的三元组顺序表的位置。
//推出2个规则：
//1 cpot[1] = 1；
//2 cpot[j] = cpot[j-1]+num[j-1];
//算法思路:
//根据1,2规则计算num数组和cpot数组
//根据cpot数组的指示求转置,将A的三元组顺序表中每个元素直接放入B的三元组顺序表中的适当位置.
#define max 100
struct element
{
 int row,col; //行下标,列下标
 int value;   //值
};

struct triple
{
 int m; //行数
 int n; //列数
 int t; //非0元素个数
 element list[max+1]; //存放的位置
};

void transpos(triple &A,triple &B)
{//求A的转置数组B
 int num[max+1]; //A每列非0元素的个数 
 int cpot[max+1];//A的某列非0元素在B的三元组中的位置

 //基本信息的复制
 B.m = A.n;
 B.n = A.m;
 B.t = A.t;
 if(B.t>0)
 {//如果存在非0的元素,这个时候才需要转置
  //初始化num数组 num[j] : j列中非0元素的个数
  //初始化 
  for(j=1;j<=A.n;j++)
  {
   num[j] = 0;
  }
  //求A中每一列中非0元素的个数
  for(j=1;j<=A.t;j++)
  {
   num[A.list[j].col]++; //
  }
  //num已经初始化好了
  //初始化cpot;
  cpot[1] = 1;
  for(j=2;j<=A.n;j++)
  {
   cpot[j] = cpot[j-1] +num[j-1];
  }
  //cpot也初始化好了

  //求转置
  for(p=1;p<=A.t;p++)
  {
   j = A.list[p].col; //j是某元素的列数
   i = cpot[j];       //i是A的j列应该存储在B的三元组的首位
   B.list[i].row = A.list[p].col;
   B.list[i].col = A.list[p].row;
   B.list[i].data =A.list[p].value;
   cpot[j]++;
  }

 }
}