ちょいと検証したけど、LU分解でいけるみたい。
あとは、k<Nの場合の処理を入れれば終了やね。