std::seed_seq::generate
template< class RandomIt > void generate( RandomIt begin, RandomIt end ); |
(C++11以上) | |
この seed_seq のコンストラクタに元々提供されたデータを基にして、範囲 [begin, end) を符号なし整数値 i (0 ≤ i < 232
) で埋めます。 生成される値は、初期値が強くバイアスされていたとしても、32ビットの範囲全体に均等に分布されます。
以下のアルゴリズム (松本眞と西村拓士によるメルセンヌ・ツイスタの初期化シーケンスをベースに、2007年の斎藤睦夫によって行われた改善を取り込んだもの) が使用されます。
begin == endの場合は、何もしません。 そうでなければ、
- まず、出力範囲の各要素を値
0x8b8b8b8bに設定します。
- 以下のアルゴリズムに従って、出力範囲の要素を変換します。
For k = 0,..., m-1
where m=max(s+1, n)
and n=end-begin
and s=v.size()
and v is the private container holding the values originally provided by the constructor of this seed_seq object,
begin[k+p] += r1begin[k+q] += r2begin[k] = r2,
where p=(n-t)/2
and q=p+t
and t=(n >= 623) ? 11 : (n >= 68) ? 7 : (n >= 39) ? 5 : (n >= 7) ? 3 : (n - 1) / 2
and r1=1664525 * T(begin[k]^begin[k+p]^begin[k−1])
and T(x) = x ^ (x >> 27)
and r2=r1+s if k==0, r2=r1 + k%n + v[k-1] if 0<k<=s, r2=r1 + k%n if k>s.
For k = m,..., m+n-1,
begin[k+p] ^= r3begin[k+q] ^= r4begin[k]=r4
where r3 = 1566083941 * T(begin[k]+begin[k+p]+begin[k-1])
and r4=r3 - k%n
ただし、すべての計算は modulo 232
で行われ、出力範囲のインデックスアクセス (begin[x]) は modulo n が取られます。
引数
| begin, end | - | std::iterator_traits<>::value_type が32ビット値を格納するのに適した符号なし整数型である可変ランダムアクセスイテレータ
|
| 型の要件 | ||
-RandomIt は LegacyRandomAccessIterator の要件を満たさなければなりません。
| ||
戻り値
なし。 結果は [begin, end) の範囲に書き込まれます。
例外
begin または end の操作が例外を投げる場合に限り投げます。
例
#include <random>
#include <iostream>
int main()
{
std::seed_seq seq({1,2,3,4,5});
std::vector<std::uint32_t> seeds(10);
seq.generate(seeds.begin(), seeds.end());
// Step 1: fill with 0x8b8b8b8b
// seeds = {2341178251, 2341178251, 2341178251, 2341178251, 2341178251,
// 2341178251, 2341178251, 2341178251, 2341178251, 2341178251 }
//
// Step 2:
// n = 10, s = 5, t = 3, p = 3, q = 6, m = 10
//
// 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 }
//
// 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 n : seeds)
std::cout << n << '\n';
}
出力:
4204997637
4246533866
1856049002
1129615051
690460811
1075771511
46783058
3904109078
1534123438
1495905678