乱数を実行ごとではなく, 関数の呼び出しごとに異なる値を発生させたい.

Question

前提・実現したいこと

実行ごとではなく, その都度に異なる乱数を代入したい.

配列を10000個用意して14の乱数を発生させています.
14の値が発生する確率は, 25%ずつだと考えられますので,
1は2500個, 2は2500個, 3は2500個, 4は2500個
という風になるようにしたいのですが, 下のようになってしまいます.

発生している問題

$ ./XCSC
1は0個, 2は10000個, 3は0個, 4は0個
$ ./XCSC
1は0個, 2は10000個, 3は0個, 4は0個
$ ./XCSC
1は0個, 2は0個, 3は10000個, 4は0個
$ 

該当のソースコード

#include<stdio.h>
#include<stdlib.h>
#include<math.h>
#include<time.h>
#include"MT.h"


/*[0,1]の一様実乱数発生*/
double MT(void)
{
    init_genrand((unsigned)time(NULL));
    return genrand_real1();
}

/*乱数分類(乱数に基づいて，1～n(引数)の値を返す)*/
int cla(int n)
{
    init_genrand((unsigned)time(NULL));
    return genrand_int32()%(n)+1;
}

int main() {
    int a[10000];
    int i, j, k, l, m, n;
    k = 0;
    l = 0;
    m = 0;
    n = 0;
    for(i=0; i < 10000; i++) {
        a[i] = cla(4);
        if(a[i] == 1)
            k++;
        if(a[i] == 2)
            l++;
        if(a[i] == 3)
            m++;
        if(a[i] == 4)
            n++;
    }
    printf("1は%d個, 2は%d個, 3は%d個, 4は%d個\n" , k, l, m, n);
    return 0;
}

Answer 1

init_genrandをmain関数の先頭で一度だけ呼べば良いです。
逆に言うと複数回呼んではいけません。

Comments

[rosbergf1]

ありがとうございました。おかげで欲しい乱数が使えてよかったです。

Answer 2

init_genrand((unsigned)time(NULL));

現在時刻（秒）をシードに乱数を初期化する。
同じシードで初期化して先頭だけ使うと毎回同じ値になっちゃいますね。

ルイスさんの言うとおり一回だけ呼ぶようにすれば直ります。
スコープの見通しをきにするならstatic変数使うと綺麗に書けると思います。

Comments

[rosbergf1]

ありがとうございました。static変数の勉強にもなりました。

Answer 3

参考情報

• C言語による乱数生成

http://www.sat.t.u-tokyo.ac.jp/~omi/random_variables_generation.html#Prepare_rand

> ...
> #include <stdio.h>
> #include "MT.h"
>  
> int main(void){
>   int i;
>   init_genrand(10);
>   for(i=0;i<100;i++){  printf("%ld\n",genrand_int32());}
> }
> ...

上で init_genrand(...)を for ループの中で呼んだらどうなるかをためしてみると良いです。