完全２部グラフから、条件を満たす辺のリストを作る ２部グラフの一塊はx+y+z+w=n-1のベクトルの集まり、もう一塊がx+y+z+w=nのベクトルの集まり それぞれ、相当のベクトル数ではあるけれど、全ペアの距離を調べてリストアップする ただし、出来上がるエッジ…

頂点数を数える 頂点数はただの漸化式。C++化する必要もないだろう。してもよいけど… n <- 500 a <- rep(0,n) a[1] <- 1 for(i in 2:n){ a[i] <- (i+2)*a[i-1]/(i-1) } a Rcppを使ってみる コメント行は"//"で始まる RcppArmadilloを使うので"//Rcpp::depend…

こちらでRcppを覚えてみようか、そのための地ならし情報収集をしてみよう、という記事を書いた。 その実践編。『自分のためにRcppを使ってみる 〜RstudioでC++を覚える』 課題を以下のように設定する。 d=4次元立方格子の非負整数格子を考える 格子はハミン…