三对角矩阵的压缩存储

a_1,1   a_1,2          

a_2,1  a_2,2  a_2,3  　　　　0

　　a_3,2  a_3,3  a_3,4

　　  ...    ...     ...

　0　   　 a_n-1,n-2  a_n-1,n-1 a_n-1,n

　　　　　　　　a_n,n-1    a_n,n

三对角矩阵指n阶方阵的非零元素a_i,j聚集在主对角线及其两边的两条线上，即|i-j|≤1,其余位置均为0，如果使用二维数组进行存储，则会浪费大量空间，对此可以使用一维数组将其压缩存储。

其中元素总数为：2+3*(n-2)+2=3n-2，则可以定义一个一维数组B[3n-2]，则a_i,j在B中的位置为k（注意k从0开始）

则在a_i,j之前的元素个数为

第1行：2

第2行：3

第3行：3

...

第i-1行：3

第i行：j-i+1

则k=2+3*(i-2)+j-i+1=2i+j-3

若已知k，则a_i,j为

i=(k+1)/3+1，(k+1:向右偏移一个单位；3：每行3个元素；+1：i从1开始)

j=k-2i+3（由k的计算公式可得）

 

转载于:https://www.cnblogs.com/6wenhong6/p/9833107.html