c++ - 特征值:t 为 1s 的随机二元 vector

Question

我想计算 K*es，其中 K 是一个 Eigen 矩阵(维度 pxp)和 es 是一个 px1 随机二进制 vector ，值为 1。

例如，如果 p=5 和 t=2 一个可能的 es 是 [1,0,1,0, 0]' 或 [0,0,1,1,0]' 等等...

如何使用 Eigen 轻松生成 es？

Answer 1

我想出了一个更好的解决方案，它是 std::vector 的组合, Egien::Map和 std::shuffle .

std::vector<int> esv(p,0);
std::fill_n(esv.begin(),t,1);
Eigen::Map<Eigen::VectorXi> es (esv.data(), esv.size());
std::random_device rd;
std::mt19937 g(rd());
std::shuffle(std::begin(esv), std::end(esv), g);

这个解决方案是内存高效的(因为 Eigen::Map 不复制 esv )并且有很大的优势，如果我们想置换 es几次(就像在这种情况下)，然后我们只需要重复 std::shuffle(std::begin(esv), std::end(esv), g);

也许我错了，但这个解决方案似乎比以前的解决方案更优雅、更高效。