xiangze's sparse blog

機械学習、ベイズ統計、コンピュータビジョンと関連する数学について

2011-01-01から1年間の記事一覧

固有AKB

レコード大賞には間にあいませんが、紅白のまえにごり押しました。 こちらの画像を使いました。 http://akb48matomemo.com/wp-content/uploads/a14738fa6ac358c1f3e90347f2e32e91.jpg 手順 まず単純に48等分します。 split.cpp #include <iostream> #include <string> #include <stdio.h></stdio.h></string></iostream>…

ラマヌジャン数

アンサイクロペディアのPythonの項に感銘を受けたので書いてみました。 ramanujan1 = head [y|(y,i,j)<-pre, (x,k,l)<-pre, y==x,i/=k ,j/=l,i/=l] where pre=[(i^3+j^3,i,j)|i<-ll, j<-ll] where ll=[1..] しかし無限リストにすると数分待っても計算が終わ…

EigenでMDS

C++の行列演算ライブラリのEigenを使うと比較的簡単にアルゴリズムを実装することができます。 MDS(多次元尺度法)を実装してみました。 #include <eigen3/Eigen/Dense> #ifdef _MDS_WITH_REDSVD #include <redsvd/redsvd.hpp> #else #include <eigen3/Eigen/SVD> #endif using namespace Eigen; //d: distance matrix o</eigen3/eigen/svd></redsvd/redsvd.hpp></eigen3/eigen/dense>…