当矩阵的非零个数和位置在操作过程中变化大时，就不宜采用顺序存储结构来表示三元组的线性表。例如，在作“将矩阵B加到矩阵A上”的操作时，由于非零元的插入或删除将会引起A.data中元素的移动。为此，对这种类型的矩阵，采用链式存储结构表示三元组的线性表更为恰当。
在链表中，每个非陵园可用一个含5个域的结点表示，其中i，j和e这3个域分别表示该非零元所在的行，列和非零元的值，向右域right用以链接同一行中下一个非零元，向下域down用以链接同一列中下一个非零元。同一行的非零元通过right域链接成一个线性表，同一列中的非零元通常down域链接成一个线性链表，每一个非零元既是某个行链表中的一个结点，又是某个列链表中的一个结点，整个矩阵构成了一个十字交叉的链表。
可用两个分别存储行链表的头指针和列链表的头指针的一维数组来表示。
代码：
1 /**
2 * 十字链表
3 *
4 * 当矩阵的非零个数和位置在操作过程中变化大时，就不宜采用顺序存储结构来表示三元组的线性表。例如，在作“将矩阵B加到矩阵A上”的操作时，由于非零元的插入或删除将会引起A.data中元素的移动。为此，对这种类型的矩阵，采用链式存储结构表示三元组的线性表更为恰当。
5 *
6 * 在链表中，每个非陵园可用一个含5个域的结点表示，其中i，j和e这3个域分别表示该非零元所在的行，列和非零元的值，向右域right用以链接同一行中下一个非零元，向下域down用以链接同一列中下一个非零元。同一行的非零元通过right域链接成一个线性表，同一列中的非零元通常down域链接成一个线性链表，每一个非零元既是某个行链表中的一个结点，又是某个列链表中的一个结点，整个矩阵构成了一个十字交叉的链表。
7 *
8 * 可用两个分别存储行链表的头指针和列链表的头指针的一维数组来表示。
9 */
10
11 // 稀疏矩阵的十字链表存储表示
12
13 function OLNode(i, j, e) {
14 // 该非零元的行和列下标
15 this.i = i || 0;
16 this.j = j || 0;
17 this.e = e;
18 // 该非零元所在行表和列表的后继链域
19 this.right = null; // type: OLNode
20 this.down = null; // type: OLNode
21 }
22
23 function CrossList() {
24 // 行和列链表头指针向量基址由CreateSMatrix分配
25 this.rhead = [];
26 this.chead = [];
27 // 稀疏矩阵的行数，列数
28 this.mu = 0;
29 this.nu = 0;
30 this.tu = 0;
31 }
32 /**
33 * 矩阵初始化
34 * @param m
35 * @param n
36 * @param t
37 * @param {Array} list 二维数组，每行的元素分别是[i, j, e]
38 */
39 CrossList.prototype.createSMatrix = function (m, n, t, list) {
40 this.mu = m;
41 this.nu = n;
42 this.tu = t;
43
44 for (var row = 0; row < list.length; row++) {
45 var p = {};
46 OLNode.apply(p, list[row]);
47 var i = list[row][0];
48 var j = list[row][1];
49 var q;
50
51 if (this.rhead[i] == null || this.rhead[i].j > j) {
52 p.right = this.rhead[i];
53 this.rhead[i] = p;
54 } else {
55 // 查询在行表中的插入位置
56 for (q = this.rhead[i]; q.right && q.right.j < j; q = q.right);
57 p.right = q.right;
58 q.right = p;
59 }
60
61 if (this.chead[j] == null || this.chead[j].i > i) {
62 p.down = this.chead[j];
63 this.chead[j] = p;
64 } else {
65 for (q = this.chead[j]; q.down && q.down.i < i; q = q.down);
66 p.down = q.down;
67 q.down = p;
68 }
69 }
70 };
71
72 // 矩阵相加
73 CrossList.prototype.addMatrix = function (crossList) {
74 var hl = [];
75 //hl初始化
76 for (var j = 0; j <= this.nu; j++)
77 hl[j] = this.chead[j];
78
79 for (var i = 0; i <= this.mu; i++) {
80 //pa和pb指向每一行的第一个非0元结点，直至最后一行
81 var pa = this.rhead[i];
82 var pb = crossList.rhead[i];
83 var pre = null;
84
85 //处理B的一行，直至本行中无非0元素的结点
86 while (pb) {
87 var p, q;
88 // 新插入一个结点到pa的左侧
89 if (!pa || pa.j > pb.j) {
90 p = new OLNode(pb.i, pb.j, pb.e);
91
92 //行表的指针变化
93 if (!pre) this.rhead[p.i] = p;
94 else pre.right = p;
95
96 p.right = pa;
97 pre = p;
98
99 //列表的指针变化
100 if (hl[p.j]) {
101 // 从hl[p.j]开始找到新结点在同一列中的前驱结点，并让hl[p.j]指向它
102 for (q = hl[p.j]; q && q.i < p.i;q = q.down)
103 hl[p.j] = q;
104 }
105
106 //在列表中插入新结点，根据行数判断插入前面还是后面
107 if (!this.chead[p.j] || this.chead[p.j].i > p.i) {
108 p.down = this.chead[p.j];
109 this.chead[p.j] = p;
110 } else {
111 p.down = hl[p.j].down;
112 hl[p.j].down = p;
113 }
114
115 hl[p.j] = p;
116 pb = pb.right;
117 } else if (pa.j < pb.j) {
118 pre = pa;
119 pa = pa.right;
120 } else {
121 //当pa.j === pb.j时，将B中当前结点的值加到A中当前结点上
122 pa.e += pb.e;
123
124 //当pa.e === 0时，删除该结点
125 if (pa.e === 0) {
126 // 若无前驱结点，将第一个非0元结点置为当前结点的后继结点，
127 // 否则前驱结点的后继结点为当前结点的后继结点
128 if (!pre) this.rhead[pa.i] = pa.right;
129 else pre.right = pa.right;
130
131 p = pa;
132 pa = pa.right;
133
134 //列表的指针变化
135 if (hl[p.j]) {
136 //从hl[p.j]开始找到新结点在同一列中的前驱结点，并让hl[p.j]指向它
137 for (q = hl[p.j]; q && q.i < p.i; q = q.down)
138 hl[p.j] = q;
139 }
140
141 if (this.chead[p.j] == p)
142 this.chead[p.j] = hl[p.j] = p.down;
143 else
144 hl[p.j].down = p.down;
145 }
146
147 pb = pb.right;
148 }
149 }
150 }
151 };
152
153 var lists = [
154 [1, 4, 5],
155 [2, 2, -1],
156 [1, 1, 3],
157 [3, 1, 2]
158 ];
159 var a = new CrossList();
160 a.createSMatrix(4, 4, 4, lists);
161 console.log(a);
162
163 var lists2 = [
164 [1, 4, -5],
165 [2, 3, 1],
166 [1, 1, 3],
167 [3, 2, 2]
168 ];
169 var b = new CrossList();
170 b.createSMatrix(4, 4, 4, lists2);
171 console.log(b);
172
173 a.addMatrix(b);
174 console.log(a);
时间： 2024-10-23 02:50:40