新吉宗 俵確率からの設定看破(2)

俵確率は
設定1で1/6.97
設定6で1/6.58
たった0.39の確率差。されど0.39の差。こいつはどの程度露呈するのか。
プログラムしてみてみようと思います。
今日の実験は設定1の台を打っていて設定6以上の確率で俵が成立してしまう確率です。
実験の結果は以下の通り

1000G 25.60%
2000G 16.17%
3000G 10.51%
4000G 6.88%
5000G 4.76%
6000G 3.41%
7000G 2.40%
8000G 1.66%

子役カウントは偉大なり。設定1or6という状況ならかなりいい精度なのでは。
しかし、1以外の設定になってくるととたんに弱くなってくる可能性もありますね。

今回のプログラムソース
#include <stdio.h>
#include <stdlib.h>
#include <time.h>

int main () {
int i,loop,sinrai=0;;
double count=0.0,x,kekka;

srand((int)time(NULL));

for(loop=0;loop<100000;loop++){//100000回挑戦
for(i=0;i<8000;i++){//1000G試行
x = Random(1,65536);
if(x>=1 && x<=9402){//設定1
//printf("%d\n",count);
count=count+1.0;
}
}
//printf("%lf\n",count);
if(8000.0/count<=6.58){//設定6の値
sinrai+=1;
}
// printf("1000回の内%lf回がベルでした\n",count);
count=1.0;//reset
}
//printf("%d",sinrai);
kekka=(sinrai/100000.0)*100.0;
printf("確率は%lf パーセントです。\n",kekka);
return(0);
}

int Random(int min,int max)
{
return min + (int)(rand()*(max-min+1.0)/(1.0+RAND_MAX));
}