稀疏矩阵的普通转置与快速转置算法
一般来说,对于系数矩阵,我们使用三元组来存储。即就是将矩阵的所有非零元素的三元组存放在一个顺序表中,如图所示:
注意一个转置的前提:该顺序表是排好序的,即行优先,列其次。
一、普通转置
这种算法比较简单,也很容易想到:
算法思想:
对
M.data
从头至尾扫描:
«
第一次扫描时,将
M.data
中列号为
1
的三元组赋值到
T.data
中
«
第二次扫描时,将
M.data
中列号为
2
的三元组赋值到
T.data
中
«
依此类推,直至将
M.data
所有三元组赋值到
T.data
中
这种算法往往理解起来较为困难。
代码
//转置运算算法
void transpose2(TSMatrix M, TSMatrix T)
{//采用三元组表存储表示,求稀疏矩阵M的转置矩阵T
T.rows=M.rows; T.size=M.size;;
q=1; // q为当前三元组在T.data[ ]存储位置(下标)
for (col=1; col<=M.cols; ++col)
for (p=1; p<=M.size; ++p)
if (M.data.get(p).j= =col){
T.data.get(q).i =M.data.get(p).j;
T.data.get(q).j=M.data.get(p).i;
T.data.get(q).value=M.data.get(p).value;
++q;
}
}
二、快速转置
这种算法往往理解起来较为困难。
算法思想:
«
按M中三元组次序转置,转置结果放入T中恰当位置。
«
关键
v
预先确定M中每一列第一个非零元在T中位置
v
为确定这些位置,转置前应先求得M的每一列中非零元个数
«
实现:设两个数组
num[ ]
、
cpos[ ]
v
num[col]
:存储
M
中第
col
列非零元个数
v
cpos[col]
:存储
M
中第
col
列第一个非零元三元组在
T.data
中的位置
«
cpos[col]
的计算方法:
v
cpos[1]=1
v
cpos[col]=cpos[col-1]+num[col-1]
(2<=col <=M.nu)
代码:
//稀疏矩阵快速转置算法(C语言描述)
void Transpose (TSMatrix M, TSMatrix &T)
{ //采用三元组顺序表存储表示,求稀疏矩阵M的转置矩阵T
T.m=M.n; T.n=M.m; T.tu=M.tu; //将M的行数赋值给T的列数,将M的列数赋值给T的行数,非零元素总数也赋值过去
for (col=1; col<=M.nu; ++col) num[col]=0; //初始化假设原三元组每一列都是0
for (t=1; t<=M.tu; ++t) ++num[M.data[t].j]; //枚举每一个非0元素, M.data[t].j 为其相应的列
cpot[1]=1; //求第 col列中第一个非零元在T.data中的序号
for(col=2; col<=M.t; ++col) //<span style="color:black;"><span style="color:black;"><span style="color:black;">存储</span><span style="color:black;">M</span><span style="color:black;">中第</span><span style="color:black;">col</span><span style="color:black;">列第一个非零元三元组在</span><span style="color:black;">T.data</span><span style="color:black;">中的位置</span></span></span>
cpot[col]=cpot[col-1]+num[col-1];
for(p=1; p<M.tu; ++p) { //枚举所有的三元组
col=M.data[p].j; //将第i个三元组的列数赋值给col
q=cpot[col];
/*取出cpot数组中cpot[col]的值,注意,这里很重要:因为在顺序表中,所有的行元素的大小已经是依次排好的,
所以我们遍历到的这个三元组时,当其与其后面的拥有相同列元素的三元组进行比较的时候,它一定是最小的,
所以应该放在前面时,所以我们取完这个值之后,将cpot[col]的值+1,即可在下面的遍历中搜索到与之邻近的
值,正好下标已经+1,直接添加上去即可*/
T.data[q].i=M.data[p].j;
T.data[q].j=M.data[p].i;
T.data[q].value=M.data[p].value;//赋值操作
++cpot[col]; //下标记录+1
}
}