三元组快速转置算法

阅读本文，默认你已经看了三元组及其基本存储结构，这里只有快速转置方法

初始三元组为M，并且是以行为标准排列，
转置重点考虑列的情况。
用T来存储转制后的三元组

#define MAX 3
typedef struct {
    int i,j;
    int num;
}triple;//存放三元组中每行三个数
typedef struct{   //用triple创建三元组
    triple data[MAX+1];//data[0]保存原始数组的行列数及不为零的元素个数。
}tsmatrix;
int main()//书中要求三元组存储元素为小标。
{
    tsmatrix M,T;//此处直接使用，并未赋值，读者假装已赋值
     //快速转置
     //先统计每列元素个数，并计算得出每列第一个元素，再转置后三元组中的位置；
     //First
     int amount[M.data[0].j+1];
     for(int k=0;k<=M.data[0].j;k++)
         amount[k]=0;
     for(int k=1;k<MAX+1;k++)
         amount[M.data[k].j]++;//找到每列有多少个元素

     // second
     int Place[M.data[0].j+1];
     for(int k=0;k<=M.data[0].j;k++)
         Place[k]=0;
     Place[1]=1;
     for(int k=2;k<=M.data[0].j;k++)
         Place[k]+=(Place[k-1]+amount[k-1]);//找到原矩阵每列第一个元素所在位置。
     //Third
     for(int k=1;k<M.data[0].num;k++)
     {
         //开始转置数组

         //直接硬写代码
         T.data[Place[M.data[k].j]].i=M.data[k].j;
         T.data[Place[M.data[k].j]].j=M.data[k].i;
         T.data[Place[M.data[k].j]].num=M.data[k].num;
         Place[M.data[k].j]++;

         //或者找一个中介，这样可以减少访问次数。
         int lie=M.data[k].j;
         int pla=Place[lie];
         T.data[pla].i=M.data[k].j;
         T.data[pla].j=M.data[k].i;
         T.data[pla].num=M.data[k].num;
         Place[lie]++;//为下一个本列结点得到制作准备；
     }

}