阅读本文,默认你已经看了三元组及其基本存储结构,这里只有快速转置方法
初始三元组为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]++;//为下一个本裂解点得到制作准备;
}
}