稀疏矩阵一次定位快速转置算法

//基于矩阵的三元组表示，将矩阵A转置为矩阵B

FastTransposeTSMatrix(TSMatrix A,TSMatrix *B)

{

  int col,t,p,q;

  int num[MAXSIZE],position[MAXSIZE];

  B->len=A.len;B-N=A.m;B->m=A.n;

  if(B->len)

  {

   for(col=1;col<=A.n;col++)

     num[col]=0;

   for(t=1;t<=A.len;t++)

     num[A.data[t].col]++; //采用数组下表计数法计算每一列的非零元素个数    

    position[1]=1;

   for(col=2;col<=A.n;col++)   //求col列中第一个非零元素在B.data[ ]中的正确位置

     position[col]=position[col-1]+num[col-1];

   for(p=1;p<=A.len;p++)        //将表A从头至尾扫描一次，实现矩阵转置

    {

     col=A.data[p].col;q=position[col];

     B->data[q].row=A.data[p].col;  //1

     B->data[q].col=A.data[p].row;  //2

     B->data[q].e=A.data[p].e;   //3      1、2、3为行列互换，元素赋值

     position[col]++;      //指向下一个列号为col的非零元素在表B中的存放位置

    }

  }

}