回答編集履歴
1
追記
test
CHANGED
@@ -65,3 +65,29 @@
|
|
65
65
|
https://wandbox.org/permlink/8NOijeLGjtDs6ang
|
66
66
|
|
67
67
|
こうなります
|
68
|
+
|
69
|
+
|
70
|
+
|
71
|
+
---
|
72
|
+
|
73
|
+
|
74
|
+
|
75
|
+
**追記分への回答**
|
76
|
+
|
77
|
+
|
78
|
+
|
79
|
+
d = a[i][k] / a[k][k]なので
|
80
|
+
|
81
|
+
a[k][k] * d = a[k][k] * a[i][k] / a[k][k] = a[i][k] となり
|
82
|
+
|
83
|
+
a[i][k] -= a[k][k] * dは常に0になります。(というか、0にしようとしているわけです。)
|
84
|
+
|
85
|
+
よって j = k の場合はサボって計算せずに
|
86
|
+
|
87
|
+
`for(j = k + 1; j <= N; j++)`でループさせています。
|
88
|
+
|
89
|
+
|
90
|
+
|
91
|
+
やってる事を一口で言うと
|
92
|
+
|
93
|
+
`i`行目の式の`k`項の係数を0にした係数行列を生成しています。
|