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