麻雀において、対子が暗刻になるまでの平均順目2

質問

以前、数学カテゴリーでした質問（対子が暗刻になるまでの平均順目http://okwave.jp/qa/q9209998.html）をもとに、質問します。手計算に限界を感じたため、一念発起してC++の学習を始めて、プログラムを書いてみたのですがどうでしょうか？計算結果は13.492（順目）だったのですが、予想していた約8順に対して少し大きい値が得られてしまいました。
もっとこういう書き方をしたほうがいいよ、とか、立式が間違っているよ、とかなんでもいいので意見をいただけると嬉しいです。
あと、rand関数の精度？がいまいちだという話も聞いているので、もっといい乱数の発生させかたに関する意見も欲しいです。しかし、約5順のずれは乱数の問題だけではないような気がしますが・・

#include <iostream>
#include <stdio.h>

using namespace std;

inline double chusen()
{
return rand() % 100 ;//0～100までの値を生成
}

int main() {

//牌姿は例えば、(1)(1)(9)(9)123456789のような形とする
//ここに(1)(1)(9)(9)123456789+(1)or(9)のように、アンコができる場合を考える

int nhai[18] = { 123, 119, 115, 111, 107, 103, 99, 95, 91, 87, 83, 79, 75, 71, 67, 63, 59, 55 };//各順目における残りの牌
double tumochusen[18];//各順目にツモってこれるかどうかを判定するために用いる。rand関数によって、ランダムな値を得る
int SIMULATE = 10000;//シミュレーションの回数
int junme[10000];//アンコができた順目。終局までにできなければ18とする
double sum=0;//順目の合計値

for (int j = 0; j < SIMULATE; j++)
{
junme[j] = 18;
for (int i = 0; i <= 17; i++)
{
tumochusen[i] = chusen();
cout << "tumochusen[" << i << "]の値は" << tumochusen[i] << "です。\n";

if (tumochusen[i] < (400 / nhai[i])) {
cout << i + 1 << "順目でアンコになりました。\n";
junme[j] = i + 1;
goto end;//アンコができたら、その局は終了する
}

}

end:
cout << "シミュレーション" << j + 1 << "を終了します。\n";
}

for (int k = 0; k < SIMULATE; k++)
{
cout<<"junme["<<k<<"]の値は"<<junme[k] << "です。\n";
}

for (int l = 0; l < SIMULATE; l++) {
sum += junme[l];
}

cout << "よって、求める平均順目は" << sum / SIMULATE << "です。\n";

return 0;
}

回答No.1

f272 · ベストアンサー

if (tumochusen[i] < (400 / nhai[i])) {
というように比較していますが，これって意図通りなのだろうか？
まず400 / nhai[i]はint/intだから答えはintにしかならない。(切り捨てです。)
tumochusen[i]は，rand() % 100を計算してdoubleにしていますが，実質的にはintです。
比較は，doubleとdoubleの比較にすべきなんじゃなかろうか？

他にはrand()を使うのならsrand()は使わないの？とか，そもそもrand()を使うべきなの？とかありますが，大勢には影響ありません。

もう一点だけ。
平均巡目はsum / SIMULATEでいいのだろうか？これについてはよくわからないので明確な言い方は避けますが，全員が無作為に打牌するわけではないという事実を完全に無視していますよね。それで求めた計算にどれほどの意味があるのだろうか？
少なくとも，自分からは見えていない牌の集合を考えて，そこから各プレイヤの手牌を推測した上で，ツモる牌の確率分布を考えるといったことまでしないと，計算がかなりずれそうです。実際にどういう計算をすればよいかは聞かないでください。私も答えは持っていません。

ちなみに，あなたが書いたプログラムでも18巡目にならないものだけを取り出して，その平均を取れば9巡目くらいにはすぐなりますよ。

お礼

回答ありがとうございます。とても参考になりました。今後もよろしくお願いします。他の方からの意見も聞けるかもしれないので、回答はもうしばらくの間受け付けます

指摘1に関して：doubleとdoubleの比較に変更します。そのほうがいいと思います
指摘2に関して：srand()って、乱数列の並び？を初期化する関数でしたっけ？まだ、乱数の扱いが不慣れですので、どうしたらいいかよくわからなかったのですが、大きくは影響しないとのことなので、このままrand()を使います
指摘3に関して：平均順目はsum / SIMULATEでも、まあいいんじゃないかな、と思います。今回のシュミレーションでは他家は和了放棄のbotとして考えます。他家の手牌の推測等は今後の課題とします
指摘4に関して：暗刻ができる場合だけを取り出して、平均を取ろうとしたことがなかったので、意外でした。著者が言っている平均8順とは、これのことをいっているのかもしれませんね