スポンサーリンク
※サイト運営にサーバーは必須です※
～このサイトもエックスサーバー
を使用しています～
目次
はじめに
※目次用の記事：ギャンブルの賭け方の種類をまとめてみた
賭け方にはいくつか種類がある。
そのうちの一つである逆ダランベール法というのがある。
※ダランベール法に関してはダランベール法（ピラミッド法）の検証（c++/c言語）で紹介した
この記事では、逆ダランベール法とはなにかを説明して、プログラムによるシミュレーション結果を示す。
逆ダランベール法とは
最初の賭け金額を自分で設定する。
負けた場合：次の勝負の賭け金を、1減らす
勝った場合：次の勝負の賭け金を、1増やす
※ちなみにダランベール法は、逆ダランベール法での、負けた場合と勝った場合の操作を反対にする
＜逆ダランベール法の特徴＞
- 勝った時は、多く賭けて、負けている時は少なく賭ける、手法といえる。
- 賭けるお金の変動が±1なので、賭け金の変動が少ない。
- 勝率が50％の勝負に対して適用するのが一般的
※勝率が50％を切る勝負だと、負ける回数の方がおおくなるので、賭ける金額が小さくなる。一方で、勝率が50％を超える勝負だと、賭け金がどんどん膨らむ。勝率が50％に近い勝負に使われる手法といえる。
＜実践上の問題＞
ダランベール法（ピラミッド法）の検証（c++/c言語）でも指摘したが、シミュレーションする際、いくつかのパターンに分けられる。
勝負の下限に関する問題（賭ける金額は0になった場合）
- 勝負をいったんやめる
- 賭ける金額を1にして、勝負を続行する
勝負の上限に関する問題
- 上限を設けない
※上限を設けた場合、上限に達して、勝負をやめるかやめないかも問題になる - 上限を設ける
※例えば、（最初にスタートする金額）の2倍を最大値とする
賭け金の従う式に関する問題
- 勝敗差が最初にスタートする金額を超えない限り、(賭ける金額)=（（最初にスタートする金額）＋（これまでの勝敗差）の式に従うか、そうでないか
※勝敗差が20以上になった時、賭ける金額が増えて
(賭ける金額)=（最初にスタートする金額）+（これまでの勝敗差）+α
になる場合がある
逆ダランベール法の検証
以下では2種類のプログラムを紹介する。
この記事では、上限を設けない逆ダランベール法を主に扱う。
上限を設けた場合のプログラムは、この記事の最後の部分で、参考までに乗せておく
言語はC++で作ったが、C言語としてコンパイルしても動くと思う
※このプログラムを組んだ人間は、大学時代に少しプログラムをかじった程度の戦闘能力しかない
|
1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60 61 62 63 64 65 66 67 68 69 70 71 72 73 74 75 76 77 78 79 80 81 82 83 84 85 86 87 88 89 90 91 92 93 94 95 96 |
#include <stdio.h> #include <stdlib.h> #include <math.h> #include <time.h> int main() { int i=0; int i_max=100000;//勝負回数 double w=2.0; //賭け金の変換率(2倍以上）(win) double p=50; //勝率（50％以下）(probability) double m_min=1;//最小の賭け金額(money) double m_ini=1000;//最初の手持ちの資金(money_initial) int c;//勝ち負け差のカウント（count) c=0;//初期化 //乱数（時間による変化） srand((unsigned)time(NULL)); //ファイルの出力準備 FILE *sf; sf = fopen("simulation.dat","w");//出力するファイルの名前を指定 //エラー if(sf==NULL){ printf("ファイルオープンエラー\n"); return -1; } //初期条件の表示 printf("#倍率は%lf\t勝率は%lf\t最小の賭け金は%lf\t手持ちのお金は%lf\n",w,p,m_min,m_ini); fprintf(sf,"#倍率は%lf\t勝率は%lf\t最小の賭け金は%lf\t手持ちのお金は%lf\n",w,p,m_min,m_ini); double m_now;//現在の資金 double m_bet;//賭け金 //初期化 m_now=m_ini;//最初は最初の手持ちの資金でスタート m_bet=20;//賭け金は資金の20でスタート for(i=1;i<=i_max;i++) { double r;//乱数（rand) r=100.0*rand()/(RAND_MAX+1.0);//0から100の乱数を出す printf("%d回目\t賭け金は%lf\t",i,m_bet); fprintf(sf,"%d回目\t賭け金は%lf\t",i,m_bet); if (p<=r){//負けの場合 m_now-=m_bet; m_bet-=1; c-=1; printf("lose\t"); fprintf(sf,"lose\t"); } else{//勝ちの場合 m_now +=m_bet*(w-1.0); m_bet+=1; c+=1; printf("win\t"); fprintf(sf,"win\t"); } printf("勝負差%d\t手持ちのお金は%lf\n",c,m_now); fprintf(sf,"勝負差%d\t手持ちのお金は%lf\n",c,m_now); if(m_bet==0){//賭け金が0になった場合、+1する m_bet+=1; } if(m_now<m_bet){//手持ちの金より賭けるお金が高い場合 printf("%d回目で手持ちのお金が足りません\n",i); fprintf(sf,"%d回目で手持ちのお金が足りません\n",i); break; } }//for(i)文の終わり //sfクローズ fclose(sf); return 0; }//プログラムの終わり |
プログラムで使用している変数の説明
※マーチンゲール法（２倍賭け）の破綻までのシミュレーションで紹介したプログラムと同じ変数の置き方をしている
i_max：勝負の回数
w：賭けたお金の戻る倍率を指定。ここでは2倍を指定
p：勝率を指定。ここでは1/2の確率なので50。
m_min：かけ金を指定
m_ini：最初の資金を指定
※ここではm_min=1でm_ini=1000と、最小の賭け金の1000倍を所持していると考えている。
例えば、最小の賭け金が1000円（千円）だとするなら、手元の軍資金は、1000000円（百万円）。割と現実的な設定だと思われる。
r:乱数（0~100の乱数）
※乱数rが勝率pより下の数で収まるなら勝ちの判定がでる。
そして、勝負を続けていくうちに、負け続けることもあるだろう。
そしてついに、賭けしようにも手持ちのお金が足りなくなるかもしれない。
この場合、これ以上勝負ができなくなり、プログラムは終了する。（破綻判定）
つまり、借金はNG
※マーチンゲール法になかった変数
c：トータルの勝ち負けの差がどのくらいあるか調べる用 。
10回勝負して6勝4負けなら+2 ←（+6−4）
10回勝負して3勝7負けなら−4 ←（+3−7）
実際にプログラムを走らせると以下のような結果となる
勝負回数i_maxは100000（10万回）までとする
※上限を設けないパターンの結果。最初の賭け金額は20でスタートする。
#倍率は2.000000 勝率は50.000000 最小の賭け金は1.000000 手持ちのお金は1000.000000
1回目 賭け金は20.000000 win 勝負差1 手持ちのお金は1020.000000
2回目 賭け金は21.000000 win 勝負差2 手持ちのお金は1041.000000
3回目 賭け金は22.000000 lose 勝負差1 手持ちのお金は1019.000000
4回目 賭け金は21.000000 lose 勝負差0 手持ちのお金は998.000000
（省略）
2222回目 賭け金は29.000000 lose 勝負差-18 手持ちのお金は68.000000
2223回目 賭け金は28.000000 lose 勝負差-19 手持ちのお金は40.000000
2224回目 賭け金は27.000000 lose 勝負差-20 手持ちのお金は13.000000
2224回目で手持ちのお金が足りません
この結果からわかるのは、上限を設けない場合だと、勝負を重ねていくと、(賭ける金額)=（最初にスタートする金額）+（これまでの勝敗差）という式から外れていく。
※srand((unsigned)time(NULL));の部分で時間を参照した上で乱数を発生させている。そのため、実行するタイミングで結果が変化する。
逆ダランベール法（上限なし）の期待値シミュレーション
上のプログラムを何度も走らせて期待値や破綻する率をもとめたい
そのためのプログラムは以下のようになる
|
1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60 61 62 63 64 65 66 67 68 69 70 71 72 73 74 75 76 77 78 79 80 81 82 83 84 85 86 87 88 89 90 91 92 93 94 95 96 97 98 99 100 101 102 103 104 105 106 107 108 109 110 111 112 113 114 115 116 117 118 119 120 121 122 123 124 125 126 127 128 129 130 131 132 133 134 135 136 137 138 139 140 141 142 143 144 145 146 147 148 149 150 151 152 153 154 155 156 |
#include <stdio.h> #include <stdlib.h> #include <math.h> #include <time.h> int main() { int k; int k_max=100000;//試行回数 int i=0; int i_max=100000;//勝負回数 double w=2.0; //賭け金の変換率(2倍以上）(win) double p=50; //勝率（50％以下）(probability) double m_min=1;//最小の賭け金額(money) double m_ini=1000;//最初の手持ちの資金(money_initial) int b;//破たん回数のカウント(Bankruptcy) b=0;//初期化 int i_sum;//勝負回数の合計 double m_sum;//最終的なお金の合計 i_sum=0;//初期化 m_sum=0.0; int c_sum;//勝負差の合計 c_sum=0;//初期化 //乱数（時間による変化） srand((unsigned)time(NULL)); //ファイルの出力準備 FILE *sf; sf = fopen("simulation.dat","w");//出力するファイルの名前を指定 //エラー if(sf==NULL){ printf("ファイルオープンエラー\n"); return -1; } //初期条件の表示 printf("#倍率は%lf\t勝率は%lf\t最小の賭け金は%lf\t手持ちのお金は%lf\n",w,p,m_min,m_ini); fprintf(sf,"#倍率は%lf\t勝率は%lf\t最小の賭け金は%lf\t手持ちのお金は%lf\n",w,p,m_min,m_ini); for(k=1;k<=k_max;k++) { int i_count;//勝負回数のカウント i_count=0; double m_now;//現在の資金 double m_bet;//賭け金 //初期化 m_now=m_ini;//最初は最初の手持ちの資金でスタート m_bet=20;//賭け金は資金の20でスタート int c;//勝ち負け差のカウント（count) c=0;//初期化 for(i=1;i<=i_max;i++) { double r;//乱数（rand) r=100.0*rand()/(RAND_MAX+1.0);//0から100の乱数を出す if (p<=r){//負けの場合 m_now-=m_bet; m_bet-=1; c-=1; } else{//勝ちの場合 m_now +=m_bet*(w-1.0); m_bet+=1; c+=1; } i_count=i;//i_countを定義しないと破たんしない（セーフ）場合、iと値が1ずれてしまうため if(m_bet==0){//賭け金が0になった場合、+1する m_bet+=1; } if(m_now<m_bet){//手持ちの金より賭けるお金が高い場合 break; } }//for(i)文の終わり if(i_count<i_max) { printf("%d番目%d回目で破たん\t手持ちのお金は%lf\n",k,i_count,m_now); fprintf(sf,"%d番目%d回目で破たん\t手持ちのお金は%lf\n",k,i_count,m_now); b=b+1; } else{ printf("%d番目\t%d回目までセーフ\t手持ちのお金は%lf\n",k,i_count,m_now); fprintf(sf,"%d番目\t%d回目までセーフ\t手持ちのお金は%lf\n",k,i_count,m_now); } i_sum=i_sum+i_count; m_sum=m_sum+m_now; c_sum=c_sum+c; }//for(k)文の終わり //破たん率 printf("勝負回数は%d\n",i_max); printf("試行回数は%d\n",k_max); printf("破たん回数は%d\n",b); fprintf(sf,"勝負回数は%d\n",i_max); fprintf(sf,"試行回数は%d\n",k_max); fprintf(sf,"破たん回数は%d\n",b); double b_pro= (double)b/(double)k_max*100.0;//(double)を抜くとb/k_maxがゼロと評価されてしまう printf("破たん率は%lf％\n",b_pro); fprintf(sf,"破たん率は%lf％\n",b_pro); double i_exp=(double)i_sum/(double)k_max; double m_exp=(double)m_sum/(double)k_max; printf("破たんするまでに行える勝負回数の期待値（i_exp）は%lf\n",i_exp); printf("破たんする直前で持っているお金の期待値（m_exp）は%lf\n",m_exp); fprintf(sf,"破たんするまでに行える勝負回数の期待値（i_exp）は%lf\n",i_exp); fprintf(sf,"破たんする直前で持っているお金の期待値（m_exp）は%lf\n",m_exp); double c_exp=(double)c_sum/(double)k_max; printf("破たんした時の勝負差の期待値は%lf\n",c_exp); fprintf(sf,"破たんした時の勝負差の期待値は%lf\n",c_exp); //sfクローズ fclose(sf); return 0; }//プログラムの終わり |
プログラムで使用している変数の説明
b_pro：破綻率
※ここでいう破綻率とは手持ちのお金が、賭けに耐えられない状況を意味する。借金してお金を用意することはできないとする。
i_exp: 破綻するまで何回勝負ができるかの期待値
m_exp: 勝負が終わった段階で持っているお金の期待値
k_max:試行回数。
※k_maxを大きくすればするほど、正確な期待値が求まる。
※マーチンゲール法になかった要素
c_exp:勝敗差の期待値
以下では、k_max=100000、i_max=100000で指定。
試行回数10万回、勝負回数10万回。
結果は以下のようになる
勝負回数は100000
試行回数は100000
破たん回数は99419
破たん率は99.419000％
破たんするまでに行える勝負回数の期待値（i_exp）は4500.586720
破たんする直前で持っているお金の期待値（m_exp）は994.091290
破たんした時の勝負差の期待値は-0.041980
逆ダランベールの破綻率は99.419％。ダランベール法の破綻率も97.722％ほど。破たんするまでに行える勝負回数の期待値も逆ダランベール法の方が少ない。そのため、上限を設けない場合は、逆ダランベール法の方が破綻しやすい。
お金の期待値は994となっており、初期の1000より減っている(誤差の範囲と思われる)。
10万回勝負を行うと、ほとんどの場合がもとの1000を下回り、一部10万回勝負が続行できた場合は大きな資金を手にするといった結果となった。
ちなみに、賭け金額の上限を設けた場合は以下のような結果となる
※上限は、最初の賭け金額の2倍の40とした。
※プログラムは下の方のページで紹介している
勝負回数は100000
試行回数は100000
破たん回数は89389
破たん率は89.389000％
破たんするまでに行える勝負回数の期待値（i_exp）は20066.632950
破たんする直前で持っているお金の期待値（m_exp）は999.959200
破たんした時の勝負差の期待値は-0.031230
上限を設けた分、破綻する率は99.419％から89.389％に少し落ちた。
ちなみに、ダランベール法（上限あり）の破綻率は89.76％だった。逆ダランベール法の方がわずかに破綻率が低かったことになる。
逆ダランベール法は、負けが続いている時に、賭ける金額を減らしていくので、ダランベール法と比べて、破綻率が低く出そうな気がしないでもない。
ただ、上限なしのバージョンの結果において、ダランベール法と逆ダランベール法で大きな差がなかったことから、0.22％で有意な差があると断定することができず、単に誤差の範囲という可能性がある。
書籍の紹介
『ギャンブルの必勝法が本当に儲かるかプログラミングで検証してみた
●勝率が50%の場合、利益を生み出す必勝法は存在するか？
●勝率が60%の場合、どのように賭けるのが最適か？
必勝法と思われている手法を15種類紹介します。必勝法には、例えば、マーチンゲール法(負けた時に2倍賭ける手法)などがあります。これらの手法が本当に儲かるかプログラミングを使用して検証します。また、検証するために必要なプログラミングの知識(C#)も紹介しています。
ギャンブルの必勝法が本当に儲かるかプログラミングで検証してみた
関連記事
※目次用の記事：ギャンブルの賭け方の種類をまとめてみた
参考
最後に、上限を設けたバージョンのプログラムを乗せておく。
|
1 2 3 |
if(m_bet>=40){//賭け金が40を超えたら、40にする（上限設定） m_bet=40; } |
ぶちゃけ、上のプログラムにこの記述を加えただけですが……
※賭ける金額の上限を40としている
逆ダランベール法（上限あり）の個別シミュレーション
|
1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60 61 62 63 64 65 66 67 68 69 70 71 72 73 74 75 76 77 78 79 80 81 82 83 84 85 86 87 88 89 90 91 92 93 94 95 96 97 98 99 100 |
#include <stdio.h> #include <stdlib.h> #include <math.h> #include <time.h> int main() { int i=0; int i_max=100000;//勝負回数 double w=2.0; //賭け金の変換率(2倍以上）(win) double p=50; //勝率（50％以下）(probability) double m_min=1;//最小の賭け金額(money) double m_ini=1000;//最初の手持ちの資金(money_initial) int c;//勝ち負け差のカウント（count) c=0;//初期化 //乱数（時間による変化） srand((unsigned)time(NULL)); //ファイルの出力準備 FILE *sf; sf = fopen("simulation.dat","w");//出力するファイルの名前を指定 //エラー if(sf==NULL){ printf("ファイルオープンエラー\n"); return -1; } //初期条件の表示 printf("#倍率は%lf\t勝率は%lf\t最小の賭け金は%lf\t手持ちのお金は%lf\n",w,p,m_min,m_ini); fprintf(sf,"#倍率は%lf\t勝率は%lf\t最小の賭け金は%lf\t手持ちのお金は%lf\n",w,p,m_min,m_ini); double m_now;//現在の資金 double m_bet;//賭け金 //初期化 m_now=m_ini;//最初は最初の手持ちの資金でスタート m_bet=20;//賭け金は資金の20でスタート for(i=1;i<=i_max;i++) { double r;//乱数（rand) r=100.0*rand()/(RAND_MAX+1.0);//0から100の乱数を出す printf("%d回目\t賭け金は%lf\t",i,m_bet); fprintf(sf,"%d回目\t賭け金は%lf\t",i,m_bet); if (p<=r){//負けの場合 m_now-=m_bet; m_bet-=1; c-=1; printf("lose\t"); fprintf(sf,"lose\t"); } else{//勝ちの場合 m_now +=m_bet*(w-1.0); m_bet+=1; c+=1; printf("win\t"); fprintf(sf,"win\t"); } printf("勝負差%d\t手持ちのお金は%lf\n",c,m_now); fprintf(sf,"勝負差%d\t手持ちのお金は%lf\n",c,m_now); if(m_bet==0){//賭け金が0になった場合、+1する m_bet+=1; } if(m_bet>=40){//賭け金が40を超えたら、40にする（上限設定） m_bet=40; } if(m_now<m_bet){//手持ちの金より賭けるお金が高い場合 printf("%d回目で手持ちのお金が足りません\n",i); fprintf(sf,"%d回目で手持ちのお金が足りません\n",i); break; } }//for(i)文の終わり //sfクローズ fclose(sf); return 0; }//プログラムの終わり |
逆ダランベール法（上限あり）の期待値シミュレーション
|
1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60 61 62 63 64 65 66 67 68 69 70 71 72 73 74 75 76 77 78 79 80 81 82 83 84 85 86 87 88 89 90 91 92 93 94 95 96 97 98 99 100 101 102 103 104 105 106 107 108 109 110 111 112 113 114 115 116 117 118 119 120 121 122 123 124 125 126 127 128 129 130 131 132 133 134 135 136 137 138 139 140 141 142 143 144 145 146 147 148 149 150 151 152 153 154 155 156 157 158 159 160 |
#include <stdio.h> #include <stdlib.h> #include <math.h> #include <time.h> int main() { int k; int k_max=100000;//試行回数 int i=0; int i_max=100000;//勝負回数 double w=2.0; //賭け金の変換率(2倍以上）(win) double p=50; //勝率（50％以下）(probability) double m_min=1;//最小の賭け金額(money) double m_ini=1000;//最初の手持ちの資金(money_initial) int b;//破たん回数のカウント(Bankruptcy) b=0;//初期化 int i_sum;//勝負回数の合計 double m_sum;//最終的なお金の合計 i_sum=0;//初期化 m_sum=0.0; int c_sum;//勝負差の合計 c_sum=0;//初期化 //乱数（時間による変化） srand((unsigned)time(NULL)); //ファイルの出力準備 FILE *sf; sf = fopen("simulation.dat","w");//出力するファイルの名前を指定 //エラー if(sf==NULL){ printf("ファイルオープンエラー\n"); return -1; } //初期条件の表示 printf("#倍率は%lf\t勝率は%lf\t最小の賭け金は%lf\t手持ちのお金は%lf\n",w,p,m_min,m_ini); fprintf(sf,"#倍率は%lf\t勝率は%lf\t最小の賭け金は%lf\t手持ちのお金は%lf\n",w,p,m_min,m_ini); for(k=1;k<=k_max;k++) { int i_count;//勝負回数のカウント i_count=0; double m_now;//現在の資金 double m_bet;//賭け金 //初期化 m_now=m_ini;//最初は最初の手持ちの資金でスタート m_bet=20;//賭け金は資金の20でスタート int c;//勝ち負け差のカウント（count) c=0;//初期化 for(i=1;i<=i_max;i++) { double r;//乱数（rand) r=100.0*rand()/(RAND_MAX+1.0);//0から100の乱数を出す if (p<=r){//負けの場合 m_now-=m_bet; m_bet-=1; c-=1; } else{//勝ちの場合 m_now +=m_bet*(w-1.0); m_bet+=1; c+=1; } i_count=i;//i_countを定義しないと破たんしない（セーフ）場合、iと値が1ずれてしまうため if(m_bet==0){//賭け金が0になった場合、+1する m_bet+=1; } if(m_bet>=40){//賭け金が40を超えたら、40にする（上限設定） m_bet=40; } if(m_now<m_bet){//手持ちの金より賭けるお金が高い場合 break; } }//for(i)文の終わり if(i_count<i_max) { printf("%d番目%d回目で破たん\t手持ちのお金は%lf\n",k,i_count,m_now); fprintf(sf,"%d番目%d回目で破たん\t手持ちのお金は%lf\n",k,i_count,m_now); b=b+1; } else{ printf("%d番目\t%d回目までセーフ\t手持ちのお金は%lf\n",k,i_count,m_now); fprintf(sf,"%d番目\t%d回目までセーフ\t手持ちのお金は%lf\n",k,i_count,m_now); } i_sum=i_sum+i_count; m_sum=m_sum+m_now; c_sum=c_sum+c; }//for(k)文の終わり //破たん率 printf("勝負回数は%d\n",i_max); printf("試行回数は%d\n",k_max); printf("破たん回数は%d\n",b); fprintf(sf,"勝負回数は%d\n",i_max); fprintf(sf,"試行回数は%d\n",k_max); fprintf(sf,"破たん回数は%d\n",b); double b_pro= (double)b/(double)k_max*100.0;//(double)を抜くとb/k_maxがゼロと評価されてしまう printf("破たん率は%lf％\n",b_pro); fprintf(sf,"破たん率は%lf％\n",b_pro); double i_exp=(double)i_sum/(double)k_max; double m_exp=(double)m_sum/(double)k_max; printf("破たんするまでに行える勝負回数の期待値（i_exp）は%lf\n",i_exp); printf("破たんする直前で持っているお金の期待値（m_exp）は%lf\n",m_exp); fprintf(sf,"破たんするまでに行える勝負回数の期待値（i_exp）は%lf\n",i_exp); fprintf(sf,"破たんする直前で持っているお金の期待値（m_exp）は%lf\n",m_exp); double c_exp=(double)c_sum/(double)k_max; printf("破たんした時の勝負差の期待値は%lf\n",c_exp); fprintf(sf,"破たんした時の勝負差の期待値は%lf\n",c_exp); //sfクローズ fclose(sf); return 0; }//プログラムの終わり |
～プログラミングを勉強してみませんか？～
TechAcademy [テックアカデミー] 
は無料の体験講座が用意されているので、気軽に体験できます。
※私(サイト主)も無料体験講座を実際に受けてみました(→感想)
コメント