参考(こちら) 正方行列を行と列の両方についてstochastic(足して１)になるように補正する処理 こちら これを使うと数独も解けるという(こちら) # sinkhorn sinkhorn<-function(A,epsilon=0.00001){ B<-A row_sum<-apply(A,1,sum) while(max(abs(t(row_sum)-1…

昨日までの記事で「すべての順列」について合算したかった 行列のパーマネントの計算であることがわかった(こちら) Wikiには小さい行列用の計算アルゴリズムRyserの方法が出ている Nはサイズ、P.Setは{1,2,...,N}のべき集合 Rで実装しよう # Ryser Permanent…