2

我想从相对较大的人群中抽取 n 个样本而不进行替换。所以我绘制随机数并跟踪我之前的选择，这样每当我两次绘制一个数字时我都可以重新采样：

boost::mt19937 generator;
boost::uniform_int<> distribution(0, 1669 - 1);
boost::variate_generator<boost::mt19937, boost::uniform_int<> >
        gen(generator, distribution);
int n = 100;

std::vector<int> idxs;
while(static_cast<int>(idxs.size()) < n)
{
    // get random samples
    std::generate_n(std::back_inserter(idxs), n - idxs.size(),
                    gen);
    // remove duplicates
    // keep everything that's not duplicates to save time
    std::sort(idxs.begin(), idxs.end());
    std::vector<int>::iterator it = std::unique(idxs.begin(), idxs.end());
    idxs.resize(std::distance(idxs.begin(), it));
}

不幸的是，我遇到了上面使用的常量的无限循环。

我添加了一些输出（表明它一直选择相同的数字）并在 10 次尝试后停止以显示问题：

boost::mt19937 generator;
boost::uniform_int<> distribution(0, 1669 - 1);
boost::variate_generator<boost::mt19937, boost::uniform_int<> >
        gen(generator, distribution);
int n = 100;

int repeat = 0;
std::vector<int> idxs;
while(static_cast<int>(idxs.size()) < n)
{
    if(repeat++ > 10) break;
    cout << "repeat " << repeat <<
            ", " << idxs.size() << " elements" << endl;
    std::generate_n(std::back_inserter(idxs), n - idxs.size(),
                    gen);
    cout << "last " << idxs.back() << endl;
    std::sort(idxs.begin(), idxs.end());
    std::vector<int>::iterator it = std::unique(idxs.begin(), idxs.end());
    idxs.resize(std::distance(idxs.begin(), it));
}

代码打印

repeat 1, 0 elements
last 1347
repeat 2, 99 elements
last 1359
repeat 3, 99 elements
last 1359

等等，如果我不杀死程序，这似乎永远循环。这不应该发生，对吧？我只是不走运吗？还是我做错了什么？

简短的解决方案 感谢@jxh！使用参考有助于：

boost::variate_generator<boost::mt19937&, boost::uniform_int<> >
        gen(generator, distribution);
4

1 回答 1

1

问题是创建了您创建generate_n的生成器的副本。gen因此，在对 的调用结束时generate_n， 的状态gen没有改变。因此，每次重新循环时，都会再次生成相同的序列。

解决此问题的一种方法是在您的variate_generator: *中使用对随机数生成器对象的引用

boost::variate_generator<boost::mt19937&, boost::uniform_int<> >
        gen(generator, distribution);

* 由于我对 Boost 的经验有限，我最初的建议相当笨拙。我在这个答案中采用了提问者实施的解决方案。

于 2013-08-19T22:11:44.147 回答