稀疏矩阵的两个算法：

矩阵的转置：

 用了两个数组 num[ ],cpot[ ]

num[i]表示第i 列有多少元素；cpot[i]表示第i 列的第一个元素转置后的位置。

cpot[1]=1;

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

status FastTrans(TSMatrix M,TSMatrix&T)

{

   T.mu=M.nu;T.nu=M.mu;T.tu=M.tu;

   if(T.tu=0) return OK;

   for(col=1;col<=M.nu;col++)

   num[col]=0;

   for(t=1;t<=M.tu;t++)

   num[M.data[t].j]++;

   cpot[1]=0;

   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].value=M.data[p].value;

  ++cpos[col];

  }

return OK;

}