C
CoreyWhite
/*
It is possible to use martingale probability theory to beat some games
of chance. In a fair game of coin toss, where the odds reach an
equilibrium of 50/50 chain reactions do occur. This can be explained
using martingale probability theory, but in simpler terms it only shows
an example of how order emerges out of chaos.
Example: One player has 3 pennies, and another player has only 1 penny.
A fair coin is tossed every round to determine if a penny is won or
lost for either player. The odds are 3/4 that player A (Who begins with
3 pennies) will win the game. This is entirely different than the
martingale betting strategy, because only 1 penny is bet for each round
of the game.
Because there are 3 ways player A may win, and only one way player B
can win, player A has a concrete advantage. Player B, only wins in the
event that the coin is tossed in his favor 3 times in a row, while
player A can win on the first throw. Or he can win after losing the
first coin toss, or he can win after losing the second coin toss. So
the odds are 75% that he (or she) will win in this game.
Upon further analysis it is possible to calculate the average number of
coin flips before player A is likely to win. The equation k(n-k) works
for perfectly fair games according to martingale probability theory to
solve this problem. In this case 3(4-3) solves the problem, so on
average it takes 3 coin flips for player A to win.
To show that chain reactions occur you only have to move from the
probability of winning the first game, and multiply it by the
probabilities of winning the following games. For example, if 3 pennies
are used to play this game in an attempt to win one penny, the odds are
3/4. And once that penny is collected there is now a 4/5th chance of
winning another penny.
So statistics tells us that there is a (3/4) * (4/5) * (5/6) * (7/8) *
(8/9) * (9/10) = 30% chance of the 3 pennies growing into a pile of 10.
But in repeatable tests you will find that on average there is not a
net win or loss in this game. If there is a 75% chance of winning 1
penny, and a 25% chance of losing 3. The two odds cancel each other
out, to create an equilibrium in 50/50 games.
And at the same time we can see that despite the fact that the initial
value of coins reaches an equilibrium when the pattern is extended to
any length, we can show a concrete advantage to beginning with 3
pennies, instead of beginning with one.
In the last example player A had a 30% chance of winning 7 pennies, and
totaling 10 in all. If we started with only one penny then player A
would just have to total 8 pennies in order to earn 7. So lets look at
the math:
(1/2) * (2/3) * (3/4) * (4/5) * (5/6) * (6/7) * (7/8) = 12.5%
So we can cleary see that even though winning 7 pennies has the same
expected value as losing 1 penny. Outside of repeatable tests the odds
of earning 7 pennies is clearly higher if you begin with 3.
I can also explain the laws of nature with these prinicples. If we
look at the equation for gravity on earth, which accelerates at 9.8 m/s
we can derive an acceptable answer from the earlier equations. The
gravity equation I am using is sqrt(2*n/9.8).
In this example we are dropping a ball from 4.9 meters, and you can see
it takes one second to land.
t = sqrt( ( 2(4.9 m) ) / ( 9.8 m/s^2 ) ) = 1 s
So here is my gravity theory. We are using the quadratic formula to
solve: 2*n/9.8 = k(n-k) , for k. (The formula k(n-k) finds the average
number of coin flips).
k=(1/14) (7n +- sqrt(49 n^2 - 40 n)).
So now an example...
We are dropping a ball from 10 meters above the ground. So we plug 10
meters into n to solve for k.
k=(1/14) (7n +- sqrt(49 n^2 - 40 n))
k=9.791574237
My question to calculate the average number of coin flips in my game is
k(n-k), so we plug in k & n:
k*(10-k) = 2.040816327 = average number of coin flips
Now we take the square root of the average number of flips to get the
actual time it takes to land:
sqrt(avg flips) = 1.428571429 = number of seconds to land.
Now finally to factor in a problem with my equation we say that if k is
9.791574327, that means our large gravity pile is that many pennies.
And our small gravity pile is exactly 0.208425673 pennies!
Now for the source code. You can actually prove everything I have
written by running a few simple test cases. In the program when you
set the initial beans to 5, and set 1, 2, 3, 4, or 5 beans as your
goal, the output should look like this:
5:1 = 83.5%
5:2 = 71.6%
5:3 = 62.6%
5:4 = 56.1%
5:5 = 50%
But if you only play with 5 beans every time and only go after 1 bean
with those 5 each game, then your output will look like this;
5:1 = 83.5%
5:1 X 5:1 = 69%
5:1 X 5:1 X 5:1 = 58%
5:1 X 5:1 X 5:1 X 5:1 = 48%
5:1 X 5:1 X 5:1 X 5:1 X 5:1 = 40%
So there is all the proof you need. Which experiment would you rather
play?
Another experiment you can try:
Modify the program to run 10,000,000. Starting with 3 beans each time
with a target of 5000. You will win 5000 beans 273 times, for winnings
of 1,365,000 beans in total. And you would lose 3 beans 443323 times
for a loss of 1,329,969 beans. So you ended up 35,031 beans ahead.
Try changing the seed and you will still be ahead in the long run.
Even if you play 100 million games you will still be ahead.
*/
.. #include <stdio.h>
.. #include <stdlib.h>
..
..
.. main ()
.. {
.. double r;
.. long int M;
.. double x;
.. int y;
.. int z;
.. int count;
..
..
.. int seed = 10000;
.. printf("Enter seend for RNG: ");
.. scanf("%d", &seed);
.. srand (seed);
.. M = 2;
..
..
.. int score = 0;
..
.. //Score keeps track of the number of beans won every game
..
..
.. int games = 0;
..
.. // games keeps track of the number of games we have played before
.. //losing all of the beans, which is equal to score.
..
..
.. int beans1 = 0;
..
.. // Initial value set to zero and defined within the loop
..
.. int wins = 0;
.. int lost = 0;
.. int quit = 0;
.. int init = 0;
.. int rounds = 0;
.. int live = 0;
..
.. printf ("Initial Beans: ");
.. scanf ("%d", &init);
.. printf ("Stop after winning X number of beans: ");
.. scanf ("%d", &quit);
.. printf ("Number of rounds: ");
.. scanf("%d", &rounds);
.. printf("Show live output (1 or 0): ");
.. scanf("%d", &live);
..
.. for (int cnt = 0; cnt < rounds; cnt++)
.. {
.. // We play up to (int) rounds
..
..
.. int count = 0;
.. beans1 = init + score;
..
.. // Beans gets defined here, as starting with 3 beans
.. // and having a 0 bonus score (It changes as you
.. // win more beans per round)
..
..
.. int beans2 = 1;
..
.. // The program attempts to win just one
.. // bean for every game.
..
..
.. while (beans1 != 0 && beans2 != 0)
..
.. // The battle begins
..
..
.. {
.. r = ((double) rand () / ((double) (RAND_MAX) + (double) (1)));
..
..
.. x = (r * M);
.. y = (int) x;
..
.. z = y + 1;
..
.. // A coin is flipped and is either 1 or 2 in value
..
.. if (z == 1)
.. {
.. // Heads wins.
..
.. beans1++;
..
.. // Beans1 gains one bean from Beans2
..
.. beans2--;
.. }
.. if (z == 2)
.. {
.. // Tails loses
..
.. beans1--;
..
.. // Beans2 gains one bean from Beans1
..
.. beans2++;
.. }
..
.. count++;
..
.. // We keep track of the number of rounds in the battle
..
.. }
..
..
.. if (beans1 > score + init)
.. {
.. // If beans1 is greater than the initial value
.. // of beans plus the total number of beans
.. // that have been won so far in this game, then
.. // the score goes up, and we go on to the next
.. // game. We check this at the end of every game.
..
.. score++;
..
.. games++;
.. }
..
.. if (beans1 <= 0)
.. {
.. //If beans1 has lost the game and doesn't
.. //have anymore beans then we know the
.. //game is over, so we reset score, and reset
.. //games.
..
..
.. if(live==1){
.. printf ("Lost at: %d beans , %d games.\n", score + init, games);
.. }
..
.. // And we print out the total number of
.. // games played on this trial and show the
.. // total score plus the initial value of beans.
..
.. lost++;
.. score = 0;
.. games = 0;
..
.. }
..
.. if (score >= quit)
.. {
.. wins++;.
..
.. if(live==1){
.. printf ("Won at: %d beans , %d games.\n", score + init, games);
.. }
..
.. beans1 == 0;
.. score = 0;
.. games = 0;
..
.. }
..
.. }
..
.. printf ("Total Won: %d/%d\n", wins, wins + lost);
.. printf ("Net win: %d beans.\n",(wins*quit)-( (wins+lost)*init ) );
..
.. }
It is possible to use martingale probability theory to beat some games
of chance. In a fair game of coin toss, where the odds reach an
equilibrium of 50/50 chain reactions do occur. This can be explained
using martingale probability theory, but in simpler terms it only shows
an example of how order emerges out of chaos.
Example: One player has 3 pennies, and another player has only 1 penny.
A fair coin is tossed every round to determine if a penny is won or
lost for either player. The odds are 3/4 that player A (Who begins with
3 pennies) will win the game. This is entirely different than the
martingale betting strategy, because only 1 penny is bet for each round
of the game.
Because there are 3 ways player A may win, and only one way player B
can win, player A has a concrete advantage. Player B, only wins in the
event that the coin is tossed in his favor 3 times in a row, while
player A can win on the first throw. Or he can win after losing the
first coin toss, or he can win after losing the second coin toss. So
the odds are 75% that he (or she) will win in this game.
Upon further analysis it is possible to calculate the average number of
coin flips before player A is likely to win. The equation k(n-k) works
for perfectly fair games according to martingale probability theory to
solve this problem. In this case 3(4-3) solves the problem, so on
average it takes 3 coin flips for player A to win.
To show that chain reactions occur you only have to move from the
probability of winning the first game, and multiply it by the
probabilities of winning the following games. For example, if 3 pennies
are used to play this game in an attempt to win one penny, the odds are
3/4. And once that penny is collected there is now a 4/5th chance of
winning another penny.
So statistics tells us that there is a (3/4) * (4/5) * (5/6) * (7/8) *
(8/9) * (9/10) = 30% chance of the 3 pennies growing into a pile of 10.
But in repeatable tests you will find that on average there is not a
net win or loss in this game. If there is a 75% chance of winning 1
penny, and a 25% chance of losing 3. The two odds cancel each other
out, to create an equilibrium in 50/50 games.
And at the same time we can see that despite the fact that the initial
value of coins reaches an equilibrium when the pattern is extended to
any length, we can show a concrete advantage to beginning with 3
pennies, instead of beginning with one.
In the last example player A had a 30% chance of winning 7 pennies, and
totaling 10 in all. If we started with only one penny then player A
would just have to total 8 pennies in order to earn 7. So lets look at
the math:
(1/2) * (2/3) * (3/4) * (4/5) * (5/6) * (6/7) * (7/8) = 12.5%
So we can cleary see that even though winning 7 pennies has the same
expected value as losing 1 penny. Outside of repeatable tests the odds
of earning 7 pennies is clearly higher if you begin with 3.
I can also explain the laws of nature with these prinicples. If we
look at the equation for gravity on earth, which accelerates at 9.8 m/s
we can derive an acceptable answer from the earlier equations. The
gravity equation I am using is sqrt(2*n/9.8).
In this example we are dropping a ball from 4.9 meters, and you can see
it takes one second to land.
t = sqrt( ( 2(4.9 m) ) / ( 9.8 m/s^2 ) ) = 1 s
So here is my gravity theory. We are using the quadratic formula to
solve: 2*n/9.8 = k(n-k) , for k. (The formula k(n-k) finds the average
number of coin flips).
k=(1/14) (7n +- sqrt(49 n^2 - 40 n)).
So now an example...
We are dropping a ball from 10 meters above the ground. So we plug 10
meters into n to solve for k.
k=(1/14) (7n +- sqrt(49 n^2 - 40 n))
k=9.791574237
My question to calculate the average number of coin flips in my game is
k(n-k), so we plug in k & n:
k*(10-k) = 2.040816327 = average number of coin flips
Now we take the square root of the average number of flips to get the
actual time it takes to land:
sqrt(avg flips) = 1.428571429 = number of seconds to land.
Now finally to factor in a problem with my equation we say that if k is
9.791574327, that means our large gravity pile is that many pennies.
And our small gravity pile is exactly 0.208425673 pennies!
Now for the source code. You can actually prove everything I have
written by running a few simple test cases. In the program when you
set the initial beans to 5, and set 1, 2, 3, 4, or 5 beans as your
goal, the output should look like this:
5:1 = 83.5%
5:2 = 71.6%
5:3 = 62.6%
5:4 = 56.1%
5:5 = 50%
But if you only play with 5 beans every time and only go after 1 bean
with those 5 each game, then your output will look like this;
5:1 = 83.5%
5:1 X 5:1 = 69%
5:1 X 5:1 X 5:1 = 58%
5:1 X 5:1 X 5:1 X 5:1 = 48%
5:1 X 5:1 X 5:1 X 5:1 X 5:1 = 40%
So there is all the proof you need. Which experiment would you rather
play?
Another experiment you can try:
Modify the program to run 10,000,000. Starting with 3 beans each time
with a target of 5000. You will win 5000 beans 273 times, for winnings
of 1,365,000 beans in total. And you would lose 3 beans 443323 times
for a loss of 1,329,969 beans. So you ended up 35,031 beans ahead.
Try changing the seed and you will still be ahead in the long run.
Even if you play 100 million games you will still be ahead.
*/
.. #include <stdio.h>
.. #include <stdlib.h>
..
..
.. main ()
.. {
.. double r;
.. long int M;
.. double x;
.. int y;
.. int z;
.. int count;
..
..
.. int seed = 10000;
.. printf("Enter seend for RNG: ");
.. scanf("%d", &seed);
.. srand (seed);
.. M = 2;
..
..
.. int score = 0;
..
.. //Score keeps track of the number of beans won every game
..
..
.. int games = 0;
..
.. // games keeps track of the number of games we have played before
.. //losing all of the beans, which is equal to score.
..
..
.. int beans1 = 0;
..
.. // Initial value set to zero and defined within the loop
..
.. int wins = 0;
.. int lost = 0;
.. int quit = 0;
.. int init = 0;
.. int rounds = 0;
.. int live = 0;
..
.. printf ("Initial Beans: ");
.. scanf ("%d", &init);
.. printf ("Stop after winning X number of beans: ");
.. scanf ("%d", &quit);
.. printf ("Number of rounds: ");
.. scanf("%d", &rounds);
.. printf("Show live output (1 or 0): ");
.. scanf("%d", &live);
..
.. for (int cnt = 0; cnt < rounds; cnt++)
.. {
.. // We play up to (int) rounds
..
..
.. int count = 0;
.. beans1 = init + score;
..
.. // Beans gets defined here, as starting with 3 beans
.. // and having a 0 bonus score (It changes as you
.. // win more beans per round)
..
..
.. int beans2 = 1;
..
.. // The program attempts to win just one
.. // bean for every game.
..
..
.. while (beans1 != 0 && beans2 != 0)
..
.. // The battle begins
..
..
.. {
.. r = ((double) rand () / ((double) (RAND_MAX) + (double) (1)));
..
..
.. x = (r * M);
.. y = (int) x;
..
.. z = y + 1;
..
.. // A coin is flipped and is either 1 or 2 in value
..
.. if (z == 1)
.. {
.. // Heads wins.
..
.. beans1++;
..
.. // Beans1 gains one bean from Beans2
..
.. beans2--;
.. }
.. if (z == 2)
.. {
.. // Tails loses
..
.. beans1--;
..
.. // Beans2 gains one bean from Beans1
..
.. beans2++;
.. }
..
.. count++;
..
.. // We keep track of the number of rounds in the battle
..
.. }
..
..
.. if (beans1 > score + init)
.. {
.. // If beans1 is greater than the initial value
.. // of beans plus the total number of beans
.. // that have been won so far in this game, then
.. // the score goes up, and we go on to the next
.. // game. We check this at the end of every game.
..
.. score++;
..
.. games++;
.. }
..
.. if (beans1 <= 0)
.. {
.. //If beans1 has lost the game and doesn't
.. //have anymore beans then we know the
.. //game is over, so we reset score, and reset
.. //games.
..
..
.. if(live==1){
.. printf ("Lost at: %d beans , %d games.\n", score + init, games);
.. }
..
.. // And we print out the total number of
.. // games played on this trial and show the
.. // total score plus the initial value of beans.
..
.. lost++;
.. score = 0;
.. games = 0;
..
.. }
..
.. if (score >= quit)
.. {
.. wins++;.
..
.. if(live==1){
.. printf ("Won at: %d beans , %d games.\n", score + init, games);
.. }
..
.. beans1 == 0;
.. score = 0;
.. games = 0;
..
.. }
..
.. }
..
.. printf ("Total Won: %d/%d\n", wins, wins + lost);
.. printf ("Net win: %d beans.\n",(wins*quit)-( (wins+lost)*init ) );
..
.. }
What I