Re: [問題] 矩陣!!限定只能用兩個矩陣...

[LPH66 ((short)(-15074))]

剛剛想到了XD

把 B 做 LUP decomposition

也就是求得一個下三角矩陣 L 一個上三角矩陣 U 及一個排列矩陣 P

使得 B = LUP

(講起來好像很複雜  其實做一次高斯消去法把它「紀錄」下來就是了 XD

 看維基: http://zh-tw.wikipedia.org/wiki/LU%E5%88%86%E8%A7%A3

 做法修改一下還可以在 B 裡面 in-place 完成!)

所以 AB 就成了 ALUP

而乘一個三角矩陣是可以 in-place 來做的

(也就是把 乘B 分成了三步 乘L 乘U 和 乘P

 其中 乘L 和 乘U 都是乘三角矩陣  可以 in-place

 乘P 只代表我們要把 A 的 column 間換來換去 XD)

--
実琴：「河野！你真的就這樣被物質慾望給吸引過去了嗎?!」
亨：「只要穿著女裝擺出親切的樣子，所有必要花費就能全免，似乎一點都不壞啊。」
実琴：「難道你沒有男人的尊嚴了嗎?!」
亨：(斷然道)「沒有。在節衣縮食且生活吃緊的學生面前，沒有那種東西。」
                                            －－プリンセス・プリンセス  第二話

--
※ 發信站: 批踢踢實業坊(ptt.cc)
◆ From: 140.112.28.92
※ 編輯: LPH66           來自: 140.112.28.92        (11/11 17:40)
推 stilltin:好像真的可以耶 ~ 推
推 HuangTzHuan:可是要把 P 存在哪裡？
→ LPH66:我如果沒搞錯的話  P的動作是可以直接把A和B做對應的對換
推 HuangTzHuan:我想問：L和U可以塞在B裡，但P是不是要另外的空間放？
推 chz:應該不用..寫成AB=APLU的話，AP可以直接做。
推 pizza0117:L 跟 U 怎麼塞進B? 對角線的空間不是重複用到嗎?
推 HudsonE:對角空間都是 1 啊
推 pizza0117:了解 3Q