使用 std::uniform_real_distribution 初始化一个 N 大小的 std::vector<double>

Initialize an N-sized std::vector with std::uniform_real_distribution<double>

本文关键字:std vector gt lt double distribution real uniform 初始化 使用 一个      更新时间:2023-10-16

我有一个std::vector<double> x,在编译时我不知道它的大小。假设向量的大小是N的,我想给它分配N均匀分布的随机值。我目前在循环中执行此操作

std::default_random_engine generator;
std::uniform_real_distribution<double> distribution_pos(0.0,1.0);
for (auto it = x.begin(); it != x.end(); it++)
{
*it = distribution_pos(generator);
}

这对我来说似乎不是很优雅,所以我想知道是否有更聪明的方法可以做到这一点?

我会使用std::generate

std::vector<double> x(10);
std::default_random_engine gen{std::random_device{}()};
std::uniform_real_distribution<double> dist(0.0, 1.0);
std::generate(std::begin(x), std::end(x), [&]{ return dist(gen); });

注意:

您需要为随机数生成器设定种子,否则每次都会得到相同的数字序列。

我不建议std::default_random_engine,因为该标准不保证跨实现或编译器版本的质量或可重现性。

您可以使用如下所示std::generate

// First create an instance of an engine.
random_device rnd_device;
// Specify the engine and distribution.
mt19937 mersenne_engine {rnd_device()};  // Generates random doubles
std::uniform_real_distribution<double> dist (0, 1.0);
auto gen = [&dist, &mersenne_engine](){
return dist(mersenne_engine);
};
std::vector<double> x(N);
std::generate(begin(x), end(x), gen);

demo

您可以使用std::generatestd::generate_n作为:

#include <algorithm>
std::generate(std::begin(x), std::end(x), [&] {
return distribution_pos(generator);
});
std::generate_n(std::begin(x), N, [&] {
return distribution_pos(generator);
});

如果您不想事先在向量中保留内存,那么这将起作用。

#include <vector>
#include <random>
#include <algorithm>
#include <iterator>
int main() {
std::vector<double> x;
const int n = 10;
std::default_random_engine gen(std::random_device{}());
std::uniform_real_distribution<double> dist(0.0, 1.0);
std::generate_n(std::back_inserter(x), n, [&]{ return dist(gen); });
}