一、问题背景
在对稀疏矩阵的三元组表使用一次定位快速转置法时,我们通常会借用两个辅助向量,即num向量和positiong向量。
其中,num向量用于统计三元组表A中第col列的非零元素的个数,并将其赋值给num[col]。
position[col]用于记录第col列中,第一个非0元素在第几行
但是用过一次定位快速转置法的同学都知道,真正在三元组表转置中起作用的,只有position向量,而num向量只是起到辅助作用,通过代码
poiston[col] = position[col-1] + num[col-1]
来计算position[col]的值
二、思路分析
让我们深入思考一下,既然最终只会用到position向量,那我们是否可以自己定义变量sum,用来记录position[col - 1]的值,用中间变量temp来记录num[col-1]的值,这样做我们就不需要num向量来全程存放所有的每一列非零元素的个数。
因此,我们可以在成功计算出num向量后,执行操作
//变量i表示从第i列开始,稀疏矩阵出现非零元素
//变量n表示稀疏矩阵的列数
for(; i<n ; i++){
temp = num[i]; //先把该列非零元素暂存temp中
num[i] = sum; //将原num[i]替换为原position[i]
sum = sum + temp; //计算原position[i+1]
}
如此操作后,我们整个过程就只使用到了num向量。
原来的算法是完整求出num向量,再完整求出position向量;在计算完后,两个向量都被保存下来。
但是在该算法中,我们计算position向量时,会让其一步一步代替num向量的元素。计算完后,num向量的值将不复存在,而position向量将以num向量的身份存活下去。