長周期乱数生成器の徹底解説：代表アルゴリズム、統計検証、並列化、用途別の選び方
エバープレイ編集部はじめに — 「長周期乱数生成器」とは何か
長周期乱数生成器（long-period pseudorandom number generator）は、その名の通り「周期（周期長）」が非常に長い擬似乱数生成器を指します。周期とは内部状態が繰り返されるまでの出力長で、周期が短いと大規模なシミュレーションや長時間のランダム列利用で周期境界に達し、統計的偏りや相関が生じます。とはいえ「長い周期＝良い乱数」ではなく、周期長は重要な要素の一つに過ぎません。本コラムでは代表的な長周期生成器、その設計思想と利点・欠点、実務上の注意点、並列化・復現性、そして用途別の選び方まで詳しく掘り下げます。
代表的な長周期生成器とその特徴
Mersenne Twister（MT19937）
1997年に松本・西村らが提案した有名な生成器で、周期は 2^19937−1、内部状態は約19937ビット（32ビット出力版では624語）です。非常に長い周期と高次元の均等分布性（32ビット出力で623次元までの均等分布性）が特徴で、科学技術計算や統計シミュレーションで長く使われてきました。欠点としては、線形性（GF(2) 上の線形再帰）に起因する低次の線形関係が残ること、状態サイズが大きくジャンプ・アヘッドや軽量な並列ストリーム作成がやや扱いづらいことが挙げられます。SFMT（SIMD-oriented Fast Mersenne Twister）
Mersenne Twister の改良版で、SIMD 命令に最適化され高速化とメモリレイアウトの改善を図ったものです。高速で良好な統計特性を持ちますが、根本的な線形性は残ります。WELL 系列（WELL19937 等）
WELL（Well Equidistributed Long-period Linear）生成器は MT の欠点を改善する目的で設計され、均等分布性や「初期化後の回復性（状態が部分的に悪いときの出力品質回復の速さ）」で優れることが報告されています。線形再帰に基づくため、やはり暗号用途には不向きです。xoroshiro / xorshift 系（xorshift128+、xoshiro256** など）
ビット単位のXORとシフト、ローテーションを組み合わせた非常に高速な生成器群です。小さな状態サイズで高速に乱数を生成できますが、設計や定数選定により品質差が大きく、低位ビットや特定の統計試験で弱点を示すことがあります。後続の改良版（xoshiro 系など）は出力の混合を工夫して性能を向上させています。PCG（Permuted Congruential Generator）
Melissa O'Neill が提唱した設計で、線形合同法（LCG）系の高速さと出力の「写像（permutation）」によるビット混合を組み合わせ、低位ビットの弱点を改善しています。小さい状態で良好な統計特性を示し、実用上の選択肢として人気です。Counter-based / Stateless（Random123、Philox、Threefry）
カウンタ＋鍵（key）を入力にしてブロック暗号や簡易関数で出力を得る方式で、任意の位置に直接ジャンプできる、並列化に非常に適した設計です。HPCやGPU上での大規模並列シミュレーションでよく使われます。暗号的な安全性は目的の関数次第です（特に設計による）。
周期の意味とその限界
周期が長いことは長い列を繰り返さずに利用できる点で有利ですが、周期だけで質を評価してはいけません。典型的な誤解は「周期が2^19937なら何しても大丈夫」というものです。実際には以下の点が重要です。
- 統計的性質（分布の均一性、相関の無さ、低次元の投影での偏り）
- 線形性や代数的構造による欠陥（特にGF(2)上の線形再帰系）
- ビット位置ごとの品質（低位ビットが弱いLCGなど）
- シード依存性・初期化の問題（初期化方法によっては短周期や偏りを生むことがある）
よって、用途ごとに適切な統計試験（TestU01、Dieharder 等）で検証された生成器を選ぶことが重要です。
統計試験とファクトチェック
生成器の品質を評価する主要なテストスイートには、George Marsaglia の Diehard、Dieharder（改良版）、そして L'Ecuyer と Simard による TestU01（SmallCrush, Crush, BigCrush）があります。これらは乱数列に潜む相関や分布の偏りを検出する手段として広く使われます。実務では、候補生成器がこれらのテストに合格しているか、あるいは既知の欠点が許容範囲かを確認してください。
実務上の注意点 — シード、再現性、初期化
十分で適切なシードを用意する
単に現在時刻だけ（time()）をシードにするのは、複数インスタンスを同時起動する用途や並列実行時に衝突を招くことがあります。可能ならば OS の CSPRNG（/dev/urandom、getrandom）や高品質なエントロピーソースを用いて初期シードを生成してください。初期化アルゴリズムに注意
特に線形系では「全零状態」や一部の弱い状態に陥ると短周期や不良出力になることがあります。MT 等は専用の初期化手順を持ち、PCG 系では splitmix64 のような良質な seeder が推奨されます。再現性を確保する
サイエンス用途では同じランを再現できることが重要です。利用したシードと生成器のバージョン（アルゴリズムの実装や定数）を必ず保存してください。
並列化とストリーム分割
大規模並列シミュレーションでは、独立した乱数ストリームを多数用意する必要があります。代表的な方法は次の通りです。
- パラメータ化：異なる定数列（パラメータ）で生成器インスタンスを構築する（ただし相互独立性の保証が難しい）
- ブロックスプリッティング / Leapfrog：一定間隔ごとに出力を分配する方法（ジャンプ・アヘッドが可能なら有効）
- カウンタベース（Random123 等）：state を持たないため、任意のスレッドが異なるカウンタ範囲を与えられて安全に独立ストリームを生成できる
並列環境では counter-based な設計や、容易にジャンプ・アヘッドできる生成器が実用的です。MT のような大きな状態を持つ生成器でも行列冪乗などでジャンプ・アヘッドは可能ですが計算コストがかかります。
暗号用途と非暗号用途の区別
多くの長周期生成器は「擬似乱数」であり、暗号的に安全（CSPRNG）ではありません。攻撃者が内部状態の一部を知れば将来の出力を予測できる可能性があります。暗号・セキュリティ用途では、ChaCha20 ベースや AES-CTR、OS 提供の暗号用乱数（getrandom, /dev/urandom）など、暗号的証明に基づく生成器を使ってください。ゲーム、シミュレーション、統計学的用途では非暗号PRNG（MT, PCG, xoshiro など）が高速で適していますが、セキュリティは考慮しないでください。
用途別の選び方（実務ガイド）
- 学術的シミュレーション（Monte Carlo 等）：TestU01 等で検証済みの高品質生成器（MT19937、SFMT、WELL、PCG、xoshiro 系の適切なパラメータ）を選ぶ。再現性確保のためシード管理を厳密に。
- 高性能並列計算（HPC/GPU）：Random123（Philox/Threefry）や counter-based アプローチ、あるいはジャンプ・アヘッドが効率的にできる生成器。
- 組み込みや低リソース環境：状態サイズが小さく高速な xorshift 系や PCG 系。ただし統計特性は用途に合わせて検証する。
- 暗号・セキュリティ：ChaCha20、AES-CTR、OS の CSPRNG を使用。擬似乱数生成器（MT 等）は不可。
- ゲーム・プロシージャル生成：再現性と速度が重要。PCG や xoshiro 系がよく選ばれる。共有シードの扱いに注意。
実践的なチェックリスト
- 用途（速度・並列性・再現性・安全性）を明確にする
- 生成器の既知の欠点（低位ビット、線形性、初期化）を把握する
- TestU01 / Dieharder 等で必要なら追加検証を行う
- シード源はできるだけ高品質なエントロピーを使う。複数実行でぶつからない対策をする
- 長期実行や大量サンプルが必要なら「周期」だけでなく統計的性質を評価する
まとめ
長周期乱数生成器は大規模シミュレーションや長時間の利用に適した基盤を提供しますが、周期長だけで安心してはいけません。線形性、ビット単位の特性、初期化の弱点、並列化のしやすさ、そして暗号的安全性の有無といった複数の観点から評価する必要があります。用途に応じて MT19937 や SFMT、WELL のような伝統的な長周期生成器、あるいは PCG や xoshiro、Random123 といった近年の設計を使い分けるのが実務的です。必ず既存の統計試験結果や公開された評価・実装を参照し、適切にシード管理を行ってください。
参考文献
- Mersenne Twister — 松本 昌樹, 西村 浩（公式ページ）
- M. Matsumoto and T. Nishimura, "Mersenne Twister: A 623-dimensionally equidistributed uniform pseudorandom number generator"
- WELL generator（Panneton, L'Ecuyer, Matsumoto）情報
- PCG — Permuted Congruential Generators（Melissa O'Neill）
- xorshift / xoshiro 系（Vigna）
- TestU01 — L'Ecuyer & Simard
- Diehard (George Marsaglia)
- Random123（counter-based RNG、GitHub）
- RFC 7539 — ChaCha20 and Poly1305
- getrandom(2) — Linux のランダム API
投稿者プロフィール
