使用 boost::random 作为 std::random_shuffle 的 RNG

我有一个程序，使用来自 boost::random 的 mt19937 随机数生成器。我需要进行 random_shuffle 并希望为此生成的随机数来自此共享状态，以便它们可以相对于梅森扭曲器先前生成的数字具有确定性。

我尝试过这样的事情：

void foo(std::vector<unsigned> &vec, boost::mt19937 &state)
{
    struct bar {
        boost::mt19937 &_state;
        unsigned operator()(unsigned i) {
            boost::uniform_int<> rng(0, i - 1);
            return rng(_state);
        }
        bar(boost::mt19937 &state) : _state(state) {}
    } rand(state);

    std::random_shuffle(vec.begin(), vec.end(), rand);
}

但我收到一个模板错误，用 rand 调用 random_shuffle 。然而这有效：

unsigned bar(unsigned i)
{
    boost::mt19937 no_state;
    boost::uniform_int<> rng(0, i - 1);
    return rng(no_state);
}
void foo(std::vector<unsigned> &vec, boost::mt19937 &state)
{
    std::random_shuffle(vec.begin(), vec.end(), bar);
}

可能是因为它是一个实际的函数调用。但显然这并不能阻止原始梅森扭曲器的状态。是什么赋予了？有没有办法在没有全局变量的情况下完成我想做的事情？

在评论中，罗伯特·古尔德要求为后代提供一个工作版本：

#include <algorithm>
#include <functional>
#include <vector>
#include <boost/random.hpp>

struct bar : std::unary_function<unsigned, unsigned> {
    boost::mt19937 &_state;
    unsigned operator()(unsigned i) {
        boost::uniform_int<> rng(0, i - 1);
        return rng(_state);
    }
    bar(boost::mt19937 &state) : _state(state) {}
};

void foo(std::vector<unsigned> &vec, boost::mt19937 &state)
{
    bar rand(state);
    std::random_shuffle(vec.begin(), vec.end(), rand);
}