- ベストアンサー
高精度乱数関数
Fookyの回答
- Fooky
- ベストアンサー率71% (59/82)
どうもしつこいようですが…。 nakashiさんのあげられてるようなページを見て、 一般的なアルゴリズムで乱数発生器をつくるのが、 最もお勧めだと私も思います。 > long doubleが10進で35桁、32768が5桁だから、 > 35/5=7と出してしまったのですが、内部では2進なので、 > 実際には何次くらいまで計算するべきかを考えると、 > 方法がわかりません。 についてですが、long doubleの最大値＋１を 10進数Ｍとすると、可能な繰り返し回数は、 log_2(M) / log_2(32768) で表されます。 Ｍは１０進数で３５桁なので、 34 <= log_10(M) < 35 ですが、今回はlog_2(M) / log_2(32768)が取りうる 最大の整数を求めたいので、上限値だけ考えればよく、 log_10(M) < 35 log_p(q) / log_p(r) = log_r(q)より、 両辺をlog_10(2)で割ると、2進数の桁数に直り、 log_2(M) < 35 / log_10(2) log_2(M) < 116.3 一方、log_2(32768)は１５。 よって、両辺を１５で割ると、 log_2(M) / log_2(32768) < 7.75 となりますので、最大の整数を取ると７。 よって、long doubleが本当に35桁であるなら、 ７回で正解だと思います。
関連するQ&A
- 標準正規分布の乱数
RAND()関数は ((double)rand() / (1.0 + RAND_MAX))と定義します。 中心極限定理により、一様乱数を足し合わせると正規分布に近づくことから、 x = 分散 * (Σ[1～12]RAND() - 6) + 平均 で正規乱数が作れる。標準正規分布は分散1、平均0なのでその乱数は x = Σ[1～12]RAND() - 6 ですよね。この乱数を例えば100個羅列するにはどうしたらいいのでしょうか？ もし間違ってたら指摘してください。 参考文献「Cによるシミュレーションプログラム 石川宏」 #include <stdio.h> #include <stdlib.h> #define RAND() ((double)rand() / (1.0 + RAND_MAX)) #define NUMBER 10000 /* 発生させる乱数の数 */ main(void) { int j; double u, x; srand(5); for (j = 0; j <= 11; j++) { u = u + RAND(); } x = u - 6.0; }
- ベストアンサー
- C・C++・C#
- 大きな数の乱数を作るには
C 初心者です。 表題のように、unsigned longのスケールの乱数をつくりたいんですが、以下のように記述すると値がいつも同じになります。この理由と、正しく動作するにはどう直したらいいのか教えてください。 unsigned long ul; ul = 4294967295UL * rand() / (RAND_MAX + 1); 値は常に131071でした。
- ベストアンサー
- C・C++・C#
- 乱数について
乱数の分布を見るために以下のようなプログラムを書きました。 #include <stdio.h> #include <stdlib.h> #include <math.h> int main() { int i,imax, S[RAND_MAX], r; double x,y; FILE *output1; output1=fopen("random2.data","w"); imax=100000; for(i=0;i<=imax;i++){ r = rand(); S[r] += 1; } for(i=0;i<=RAND_MAX;i++){ fprintf(output1,"%d %d \n",i,S[i]); } return 0; } するとコンパイルできて実行もできるのですが、なぜか乱数が30000を 超えるくらいのところでおかしな値になりました。 原因がわからないのでどなたか教えてください。
- ベストアンサー
- C・C++・C#
- 乱数について(Ｃ言語)
C言語において，乱数の範囲を 0 ≦ r < 1 とする場合には double r=(double)rand()/(RAND_MAX+1); とするのは知っているのですが0 < r ≦ 1にする場合の方法がわからず困っています． アドバイスいただきたいです．
- ベストアンサー
- その他（インターネット・Webサービス）
- 乱数のdouble型について
JSPから下記のクラスファイルを呼び出し、戻り値を返すように作りたいのですがうまくいきません。 1.Math.floor(Math.random()*100)-50 上記で実行しても小数点以下がでてしまいます。 出ないようにするのは無理なのでしょうか？ 2.乱数はdouble型以外だめなのでしょうか？ コンパイルするとdouble型なので間接参照できません。というようなエラーが出てしまいます。 これは結果を文字型に変更して戻り値として返したいのですができません。 どうすれば理想どおりにできるようになりますか？ public String getR(){ double rnd = Math.floor(Math.random()*100)-50; if(rnd.length = 4){ rnd = rnd.substring(0,2) ; } else{ rnd = rnd.substring(0,3) ; } _R = rmd; return _R; }
- ベストアンサー
- Java
- 正規分布の乱数生成
C言語で正規分布の乱数を発生させたいのですがどうすればいいのでしょうか？ 自分なりにネットで検索して調べたのですが void gaussrand() { static double V1, V2, S; static int phase = 0; double X; if(phase == 0) { do { double U1 = (double)rand() / RAND_MAX; double U2 = (double)rand() / RAND_MAX; V1 = 2 * U1 - 1; V2 = 2 * U2 - 1; S = V1 * V1 + V2 * V2; } while(S >= 1 || S == 0); X = V1 * sqrt(-2 * log(S) / S); } else X = V2 * sqrt(-2 * log(S) / S); phase = 1 - phase; } こうありました。 例えば平均50の分散9の正規分布の乱数を1000個発生させて、配列seiki[1000]に代入したいときは、このプログラミングをどのようにすればいいのでしょうか？ もちろん、このソースではなく、他のもので説明していただけても全然構いません。 また、もしよろしければ、正規分布の他に、二項分布など他の分布でのデータの生成方法もお教えいただけたら幸いです。 よろしくお願いいたします。
- 締切済み
- C・C++・C#
- 乱数を発生させるプログラムを教えてください。
タイトルのままなのですが、１から１００までの乱数を発生させるプログラムを知りたいです。 乱数をｘとおくと、ｘの値は、０＜ｘ＜１の範囲内でお願いします。 C言語で、rand関数を用いて、どうかお願いします。
- ベストアンサー
- C・C++・C#
- 精度を上げたいのですが…
#include <stdio.h> #include <time.h> #include <stdlib.h> #define MAX 1000 main(void) { int i; float x1, x2, en, sum=0.0, s; srand( (unsigned)time( NULL ) ); for(i=0;i<MAX;i++) { x1=((float)rand()/(float)RAND_MAX); x2=((float)rand()/(float)RAND_MAX); if(en=(x1-0.5)*(x1-0.5)+(x2-0.5)*(x2-0.5)<=(0.5)*(0.5)) { sum++; } } s=sum/MAX; printf("円の面積:%15.6e\n",s); } この方法で円の面積を求めたんですが、もう少し精度を上げたいと思います。ただそのプログラムをどうやって書けばいいのかさえわからずとまどっています。円全体でなくその一部を考え、またその部分を簡単に面積が求められるようにわけるプログラムを組みたいのですがどのようにすればいいのか教えてください。
- ベストアンサー
- C・C++・C#
お礼
そうか! 対数をとればよかったのか。 2進と10進の変換については、あまり頭が回らなくて。 どうもありがとうございました。