std::seed_seq::generate
来自 cppreference.com
| template< class RandomIt > void generate( RandomIt begin, RandomIt end ); |
(自 C++11) | |
使用在该 seed_seq 的构造函数中最初提供的 i 的无符号整数值填充范围 [begin, end)，0 ≤ i < 232
。即使初始值存在严重偏差，生成的数值也将在整个 32 位范围内分布。
使用以下算法（改编自 松本真和西村拓志 的梅森旋转生成器的初始化序列，结合了 斎藤睦男在 2007 年的改进）
- 如果 begin == end，则不执行任何操作。否则，
- 首先，将输出范围的每个元素都设置为值 0x8b8b8b8b。
- 根据以下算法转换输出范围的元素
对于 k = 0,..., m - 1
其中 m = max(s + 1, n)
以及 n = end - begin
以及 s = v.size()
以及 v 是包含该 seed_seq 对象构造函数最初提供的值的私有容器，
- begin[k + p] += r1
- begin[k + q] += r2
- begin[k] = r2,
其中 p = (n - t) / 2
以及 q = p + t
以及 t = (n >= 623) ? 11 : (n >= 68) ? 7 : (n >= 39) ? 5 : (n >= 7) ? 3 : (n - 1) / 2
以及 r1 = 1664525 * T(begin[k] ^ begin[k + p] ^ begin[k − 1])
以及 T(x) = x ^ (x >> 27)
以及 r2 = r1 + s 如果 k == 0，r2 = r1 + k % n + v[k - 1] 如果 0 < k <=s，r2 = r1 + k % n 如果 k > s。
对于 k = m,..., m + n - 1,
- begin[k + p] ^= r3
- begin[k + q] ^= r4
- begin[k] = r4
其中 r3 = 1566083941 * T(begin[k] + begin[k + p] + begin[k - 1])
以及 r4 = r3 - k % n
其中所有计算均为 模 232
，并且输出范围的索引（begin[x]）为 模 n。
内容 |
[编辑] 参数
| begin, end | - | 可变随机访问迭代器，其 std::iterator_traits<>::value_type 是适用于存储 32 位值的无符号整型 |
| 类型要求 | ||
-RandomIt 必须满足 LegacyRandomAccessIterator 的要求。 | ||
[编辑] 返回值
无，结果将写入到 [begin, end) 范围。
[编辑] 异常
仅当对 begin 和 end 的操作抛出异常时才会抛出异常。
[编辑] 示例
运行此代码
#include <algorithm> #include <cassert> #include <cstdint> #include <iostream> #include <random> // Prototyping the main part of std::seed_seq... struct seed_seq { std::vector<std::uint32_t> v; seed_seq(std::initializer_list<std::uint32_t> const il) : v{il} {} template <typename RandomIt> void generate(RandomIt first, RandomIt last) { if (first == last) return; // // Assuming v = {1,2,3,4,5} and distance(first, last) == 10. // // Step 1: fill with 0x8b8b8b8b // seeds = {2341178251, 2341178251, 2341178251, 2341178251, 2341178251, // 2341178251, 2341178251, 2341178251, 2341178251, 2341178251 } // std::fill(first, last, 0x8b8b8b8b); // // Step 2: // n = 10, s = 5, t = 3, p = 3, q = 6, m = 10 // const std::uint32_t n = last - first; const std::uint32_t s = v.size(); const std::uint32_t t = (n < 7) ? (n - 1) / 2 : (n < 39) ? 3 : (n < 68) ? 5 : (n < 623) ? 7 : 11; const std::uint32_t p = (n - t) / 2; const std::uint32_t q = p + t; const std::uint32_t m = std::max(s + 1, n); // // First iteration, k = 0; r1 = 1371501266, r2 = 1371501271 // // seeds = {1371501271, 2341178251, 2341178251, 3712679517, 2341178251, // 2341178251, 3712679522, 2341178251, 2341178251, 2341178251 } // // Iterations from k = 1 to k = 5 (r2 = r1 + k % n + v[k - 1]) // // r1 = 2786190137, 3204727651, 4173325571, 1979226628, 401983366 // r2 = 2786190139, 3204727655, 4173325577, 1979226636, 401983376 // // seeds = {3350727907, 3188173515, 3204727655, 4173325577, 1979226636, // 401983376, 3591037797, 2811627722, 1652921976, 2219536532 } // // Iterations from k = 6 to k = 9 (r2 = r1 + k % n) // // r1 = 2718637909, 1378394210, 2297813071, 1608643617 // r2 = 2718637915, 1378394217, 2297813079, 1608643626 // // seeds = { 434154821, 1191019290, 3237041891, 1256752498, 4277039715, // 2010627002, 2718637915, 1378394217, 2297813079, 1608643626 } // auto begin_mod = [first, n](std::uint32_t u) -> decltype(*first)& { return first[u % n]; // i.e. begin[x] is taken modulo n }; auto T = [](std::uint32_t x) { return x ^ (x >> 27); }; for (std::uint32_t k = 0, r1, r2; k < m; ++k) { r1 = 1664525 * T(begin_mod(k) ^ begin_mod(k + p) ^ begin_mod(k - 1)); r2 = (k == 0) ? r1 + s : (k <= s) ? r1 + k % n + v[k - 1] : r1 + k % n; begin_mod(k + p) += r1; begin_mod(k + q) += r2; begin_mod(k) = r2; } // // Step 3 // iterations from k = 10 to k = 19, using ^= to modify the output // // r1 = 1615303485, 3210438310, 893477041, 2884072672, 1918321961, // r2 = 1615303485, 3210438309, 893477039, 2884072669, 1918321957 // // seeds = { 303093272, 3210438309, 893477039, 2884072669, 1918321957, // 1117182731, 1772877958, 2669970405, 3182737656, 4094066935 } // // r1 = 423054846, 46783064, 3904109085, 1534123446, 1495905687 // r2 = 423054841, 46783058, 3904109078, 1534123438, 1495905678 // // seeds = { 4204997637, 4246533866, 1856049002, 1129615051, 690460811, // 1075771511, 46783058, 3904109078, 1534123438, 1495905678 } // for (std::uint32_t k = m, r3, r4; k < m + n; ++k) { r3 = 1566083941 * T(begin_mod(k) + begin_mod(k + p) + begin_mod(k - 1)); r4 = r3 - k % n; begin_mod(k+p) ^= r3; begin_mod(k+q) ^= r4; begin_mod(k) = r4; } } }; int main() { const auto input = std::initializer_list<std::uint32_t>{1,2,3,4,5}; const auto output_size = 10; // using std version of seed_seq std::seed_seq seq(input); std::vector<std::uint32_t> seeds(output_size); seq.generate(seeds.begin(), seeds.end()); for (const std::uint32_t n : seeds) std::cout << n << '\n'; // using custom version of seed_seq seed_seq seq2(input); std::vector<std::uint32_t> seeds2(output_size); seq2.generate(seeds2.begin(), seeds2.end()); assert(seeds == seeds2); }
输出
4204997637 4246533866 1856049002 1129615051 690460811 1075771511 46783058 3904109078 1534123438 1495905678
