Skip to content

Latest commit

 

History

History
114 lines (96 loc) · 6.91 KB

random.md

File metadata and controls

114 lines (96 loc) · 6.91 KB

#random <random>ヘッダ は擬似乱数を取り扱う為、乱数生成器 (Random Number Generator) や分布生成器 (Distribution) 及び関連する一連のクラス等を定義する C++11 で仕様化された乱数ライブラリである。

##乱数生成器

※多くのユーザーにとって以下に挙げる生の乱数生成器テンプレートクラスそのものよりも、この次に挙げるパラメータ定義済みの乱数生成器クラスを使用すれば十分である。

linear_congruential_engine 線形合同法エンジン
mersenne_twister_engine メルセンヌツイスターエンジン
subtract_with_carry_engine キャリー付き減算エンジン
discard_block_engine ディスカードブロックエンジン
independent_bits_engine インデペンデントビッツエンジン
shuffle_order_engine シャッフルオーダーエンジン

##パラメータ定義済み乱数生成器

※先に挙げた乱数生成器に対し一般的なパラメータを定義し、使い易い様に用意された乱数生成器クラス。

minstd_rand0
minstd_rand
mt19937
mt19937_64
ranlux24_base
ranlux48_base
ranlux24
ranlux48
knuth_b
default_random_engine

##分布生成器及びその他のクラス

random_device 予測不能な真の乱数値を得る為のクラス
seed_seq 疑似乱数エンジンの為のシード数列クラス
generate_canonical 実数区間[0.0-1.0)に展開(事実上正規化)された一様分布乱数を得る為のテンプレート関数

| | uniform_int_distribution | 一様分布/一様整数分布生成器 | | uniform_real_distribution | 一様分布/一様実数分布生成器 | | bernoulli_distribution | ベルヌーイ分布/ベルヌーイ分布生成器 | | binomial_distribution | ベルヌーイ分布/二項分布生成器 | | geometric_distribution | ベルヌーイ分布/幾何分布生成器 | | negative_binomial_distribution | ベルヌーイ分布/負の二項分布生成器 | | poisson_distribution | ポワソン分布/ポワソン分布生成器 | | exponential_distribution | ポワソン分布/指数分布生成器 | | gamma_distribution | ポワソン分布/ガンマ分布生成器 | | weibull_distribution | ポワソン分布/ワイブル分布生成器 | | extreme_value_distribution | ポワソン分布/極値分布生成器 | | normal_distribution | 正規分布/正規分布生成器 | | lognormal_distribution | 正規分布/対数正規分布生成器 | | chi_squared_distribution | 正規分布/χ二乗分布生成器 | | cauchy_distribution | 正規分布/コーシー分布生成器 | | fisher_f_distribution | 正規分布/フィッシャーの基準分布生成器 | | student_t_distribution | 正規分布/ステューデントのt分布生成器 | | discrete_distribution | 標本分布/整数のインデックスごとに離散した確率分布生成器 | | piecewise_constant_distribution | 標本分布/区間ごとの重み付けを定数値とした分布生成器 | | piecewise_linear_distribution | 標本分布/区間ごとの重み付けを線形に接続した分布生成器 |

##例

以下に示す例では、C++11の標準の乱数ライブラリを用いてランダムデバイスから初期化数列を取得してメルセンヌツイスタエンジンを初期化、単精度浮動少数型で[-1.0f - 1.0f]の一様分布及び1.0fを中心として標準偏差0.5fの正規分布に基づく擬似乱数を1メガ個生成し"random.tsv"にタブ区切り形式のファイルとして結果を保存する。

#include <iostream>
#include <exception>
#include <random>
#include <algorithm>
#include <functional>
#include <fstream>

int main()try{
  static const size_t seed_size = 8;
  typedef std::random_device device_type;
  typedef std::mt19937_64 engine_type;
  typedef std::uniform_real_distribution<float> distribution_type_1;
  typedef std::normal_distribution<float> distribution_type_2;
  auto s = [seed_size](){
    device_type r;
    std::vector<device_type::result_type> i(seed_size);
    std::generate(i.begin(), i.end(), std::ref(r));
    return std::seed_seq(i.begin(), i.end());
  }();
  engine_type e(s);
  distribution_type_1 d1(-1.0f, 1.0f);
  distribution_type_2 d2(1.0f, 0.5f);
  std::ofstream o("random.tsv");
  for(size_t n = 1000*1000; n; --n)
    o << d1(e) << "\t" << d2(e) << "\n";
  o.close();
}catch(const std::exception& e){
  std::cerr << e.what();
}
  • random[color ff0000]
  • std::random_device[color ff0000]
  • std::mt19937_64[color ff0000]
  • std::uniform_real_distribution[color ff0000]
  • std::normal_distributio<[color ff0000]
  • std::seed_seq[color ff0000]

この例である時得られた random.tsv (ファイルサイズが大きいので添付する上では random.tsv.xz に圧縮) を元に、得られたデータの密度を図示すると、以下のような図が得られた。

緑のラインは d1 (一様分布; min=-1.0f, max=1.0f) 、赤のラインは d2 (正規分布; mean=1.0f, stdev=0.5f) 、横軸は値、縦軸は密度(値の件数を区間ごとに数えたヒストグラムを全体に占める割合で表したもの)である。

参照: さまざまな確率分布 (probability distributions)