Rには、ベクタxのそれぞれに確率を割り振るsample()という関数があると思うのですが (sample(x, size, replace = FALSE, prob = NULL))、同じことをC++ですることは (ライブラリを使うなどすれば)可能でしょうか。
Eigenライブラリでは、すべての要素に均一の確率しか割り振れないようでしたが、工夫すればできるのでしょうか。