【C言語】偏りのない乱数を生成する方法を解説

2023-09-132024-07-24

この記事では、C言語を使って乱数を生成する方法について詳しく解説します。

乱数は、ゲームやシミュレーションなどでよく使われる重要な要素です。

初心者の方でも理解しやすい内容になっていますので、ぜひ参考にしてください。

目次から探す

C言語における乱数生成の基本

C言語では、乱数を生成するために標準ライブラリを利用します。

乱数は、ゲームやシミュレーション、統計的な処理など、さまざまな場面で必要とされる重要な要素です。

ここでは、C言語における乱数生成の基本について解説します。

標準ライブラリの利用

C言語で乱数を生成するためには、stdlib.hという標準ライブラリを使用します。

このライブラリには、乱数生成に関する関数が含まれており、プログラム内で簡単に乱数を扱うことができます。

stdlib.hの役割

stdlib.hは、C言語の標準ライブラリの一部で、メモリ管理や乱数生成、変換などの機能を提供します。

乱数生成に関しては、特にrand()とsrand()という2つの関数が重要です。

これらの関数を使うことで、簡単に乱数を生成し、必要に応じてその挙動を制御することができます。

乱数生成関数の紹介

C言語では、主に以下の2つの関数を使用して乱数を生成します。

rand()関数

rand()関数は、0からRAND_MAX(通常は32767)までの整数の乱数を生成します。

この関数を呼び出すたびに、異なる乱数が返されますが、初期化を行わない場合、同じシーケンスの乱数が生成されることがあります。

以下は、rand()関数を使った簡単な例です。

#include <stdio.h>
#include <stdlib.h>
int main() {
    // 乱数を生成する
    int random_number = rand();
    printf("生成された乱数: %d\n", random_number);
    return 0;
}

このプログラムを実行すると、毎回異なる乱数が表示されますが、初期化を行わない場合、同じ乱数が生成されることがあります。

srand()関数

srand()関数は、乱数生成の初期化を行うための関数です。

この関数にシード値を渡すことで、rand()関数が生成する乱数のシーケンスを制御できます。

シード値が同じであれば、rand()関数は同じ乱数のシーケンスを生成します。

以下は、srand()関数を使って乱数を初期化する例です。

#include <stdio.h>
#include <stdlib.h>
#include <time.h>
int main() {
    // 現在の時刻をシード値として使用
    srand(time(NULL));
    // 乱数を生成する
    int random_number = rand();
    printf("生成された乱数: %d\n", random_number);
    return 0;
}

このプログラムでは、time(NULL)を使って現在の時刻をシード値として設定しています。

これにより、プログラムを実行するたびに異なる乱数のシーケンスが生成されるようになります。

以上が、C言語における乱数生成の基本的な部分です。

次のセクションでは、偏りのない乱数を生成するための方法について詳しく解説します。

偏りのない乱数を生成するための方法

乱数を生成する際に、偏りのない乱数を得るためには、いくつかの重要なポイントがあります。

ここでは、乱数の初期化や範囲指定の方法について詳しく解説します。

乱数の初期化

乱数を生成する前に、初期化を行うことが重要です。

初期化を行うことで、毎回異なる乱数列を生成することができます。

シード値の重要性

乱数生成器は、内部的にシード値を使用して乱数を生成します。

シード値が同じであれば、生成される乱数列も同じになります。

したがって、シード値を毎回異なる値に設定することで、異なる乱数列を得ることができます。

一般的に、シード値には現在の時刻を使用することが多いです。

これにより、プログラムを実行するたびに異なるシード値が設定され、偏りのない乱数を生成することが可能になります。

srand()の使い方

rand()関数を使用する前に、srand()関数を使ってシード値を設定します。

以下は、srand()の使い方の例です。

#include <stdio.h>
#include <stdlib.h>
#include <time.h>
int main() {
    // 現在の時刻をシード値として設定
    srand(time(NULL));
    // 乱数を生成
    for (int i = 0; i < 5; i++) {
        int random_number = rand();
        printf("%d\n", random_number);
    }
    return 0;
}

このプログラムでは、time(NULL)を使って現在の時刻を取得し、それをシード値としてrand()関数を初期化しています。

これにより、毎回異なる乱数が生成されます。

乱数の範囲を指定する

生成される乱数の範囲を指定することも重要です。

rand()関数は、0からRAND_MAXまでの整数を生成しますが、特定の範囲の乱数が必要な場合は、スケーリングを行う必要があります。

乱数のスケーリング

特定の範囲の乱数を生成するためには、次のようにスケーリングを行います。

int lower_bound = 1; // 下限
int upper_bound = 10; // 上限
int random_number = (rand() % (upper_bound - lower_bound + 1)) + lower_bound;

この式では、rand()で生成された乱数をupper_bound - lower_bound + 1で割り、その余りを使って指定した範囲に収めています。

範囲指定の実装例

以下は、1から10の範囲で偏りのない乱数を生成するプログラムの例です。

#include <stdio.h>
#include <stdlib.h>
#include <time.h>
int main() {
    srand(time(NULL)); // シード値の設定
    // 1から10の範囲で乱数を生成
    for (int i = 0; i < 5; i++) {
        int random_number = (rand() % 10) + 1; // 1から10の乱数
        printf("%d\n", random_number);
    }
    return 0;
}

このプログラムでは、1から10の範囲で5つの乱数を生成し、出力しています。

rand() % 10で0から9の乱数を生成し、1を加えることで1から10の範囲に調整しています。

以上の方法を用いることで、C言語において偏りのない乱数を生成することができます。

シード値の設定や範囲指定を適切に行うことで、より信頼性の高い乱数を得ることができるでしょう。

より高品質な乱数生成

乱数生成にはさまざまなアルゴリズムが存在し、それぞれ特性や用途が異なります。

ここでは、より高品質な乱数を生成するためのアルゴリズムとライブラリについて解説します。

乱数生成アルゴリズムの選択

乱数生成アルゴリズムは、生成される乱数の品質に大きな影響を与えます。

以下に代表的なアルゴリズムを紹介します。

線形合同法

線形合同法(Linear Congruential Generator, LCG)は、最も古典的な乱数生成アルゴリズムの一つです。

このアルゴリズムは、次の式を用いて新しい乱数を生成します。

ここで、X_nは前回生成した乱数、a、c、mは定数です。

この方法は実装が簡単で、計算が速いという利点がありますが、周期が短く、偏りが生じやすいという欠点もあります。

メルセンヌ・ツイスタ

メルセンヌ・ツイスタ(Mersenne Twister)は、1997年に発表された高品質な乱数生成アルゴリズムです。

このアルゴリズムは、非常に長い周期(2^19937−1)を持ち、均一性や独立性が高い乱数を生成します。

メルセンヌ・ツイスタは、特にシミュレーションやモンテカルロ法など、精度が求められる場面で広く使用されています。

C言語では、メルセンヌ・ツイスタを実装したライブラリも存在しますので、必要に応じて利用することができます。

高品質な乱数生成ライブラリの紹介

高品質な乱数を生成するためには、専用のライブラリを使用するのが効果的です。

ここでは、いくつかの代表的なライブラリを紹介します。

random.hライブラリ

random.hは、C言語で高品質な乱数を生成するためのライブラリです。

このライブラリは、メルセンヌ・ツイスタを基にした乱数生成器を提供しており、簡単に高品質な乱数を生成することができます。

random.h 入手先(Github

Github – contiki-os / random.h

Github – contiki-os / random.c

以下は、random.hを使用したサンプルコードです。

#include <stdio.h>
#include "random.h" // random.hをインクルード
int main() {
    // 乱数生成器の初期化
    random_init();
    // 10個の乱数を生成して表示
    for (int i = 0; i < 10; i++) {
        printf("%d\n", random_int(1, 100)); // 1から100の範囲の乱数
    }
    return 0;
}

このコードでは、random_init()で乱数生成器を初期化し、random_int(1, 100)で1から100の範囲の乱数を生成しています。