Is there a way to really solve it?

If I win when I flip more consecutive heads than I've ever flipped tails, what's my probability of winning?

#Combinatorics

Note by Sumukh Bansal
3 years, 7 months ago

No vote yet
1 vote

  Easy Math Editor

This discussion board is a place to discuss our Daily Challenges and the math and science related to those challenges. Explanations are more than just a solution — they should explain the steps and thinking strategies that you used to obtain the solution. Comments should further the discussion of math and science.

When posting on Brilliant:

  • Use the emojis to react to an explanation, whether you're congratulating a job well done , or just really confused .
  • Ask specific questions about the challenge or the steps in somebody's explanation. Well-posed questions can add a lot to the discussion, but posting "I don't understand!" doesn't help anyone.
  • Try to contribute something new to the discussion, whether it is an extension, generalization or other idea related to the challenge.
  • Stay on topic — we're all here to learn more about math and science, not to hear about your favorite get-rich-quick scheme or current world events.

MarkdownAppears as
*italics* or _italics_ italics
**bold** or __bold__ bold

- bulleted
- list

  • bulleted
  • list

1. numbered
2. list

  1. numbered
  2. list
Note: you must add a full line of space before and after lists for them to show up correctly
paragraph 1

paragraph 2

paragraph 1

paragraph 2

[example link](https://brilliant.org)example link
> This is a quote
This is a quote
    # I indented these lines
    # 4 spaces, and now they show
    # up as a code block.

    print "hello world"
# I indented these lines
# 4 spaces, and now they show
# up as a code block.

print "hello world"
MathAppears as
Remember to wrap math in \( ... \) or \[ ... \] to ensure proper formatting.
2 \times 3 2×3 2 \times 3
2^{34} 234 2^{34}
a_{i-1} ai1 a_{i-1}
\frac{2}{3} 23 \frac{2}{3}
\sqrt{2} 2 \sqrt{2}
\sum_{i=1}^3 i=13 \sum_{i=1}^3
\sin \theta sinθ \sin \theta
\boxed{123} 123 \boxed{123}

Comments

This is a tricky one. I counted the number of allowable strings of length nn starting at n=1n = 1 and obtained the sequence 1,0,1,0,1,1,1,2,3,4,6,10,15,....1,0,1,0,1,1,1,2,3,4,6,10,15, ...., which matches this OEIS sequence. So the desired probability of winning is

P=12+123+125+126+127+228+329+4210+6211+10212+.....P = \dfrac{1}{2} + \dfrac{1}{2^{3}} + \dfrac{1}{2^{5}} + \dfrac{1}{2^{6}} + \dfrac{1}{2^{7}} + \dfrac{2}{2^{8}} + \dfrac{3}{2^{9}} + \dfrac{4}{2^{10}} + \dfrac{6}{2^{11}} + \dfrac{10}{2^{12}} + ......

Now the sequence of numerators is Fibonacci-like in that the ratio of successive terms approaches the golden ratio ϕ\phi. This allows us to estimate that 0.798<P<0.7990.798 \lt P \lt 0.799.

Edit: As Milo Koumouris has pointed out, my initial calculation was incorrect; the estimate should actually be P0.7112P \approx 0.7112.

Brian Charlesworth - 3 years, 7 months ago

Log in to reply

Thanks Sir

Sumukh Bansal - 3 years, 7 months ago

Wow - I didn't expect this! Nice... Is there away to prove the link?

Miles Koumouris - 3 years, 7 months ago

Log in to reply

While I've confirmed the match up to n=13n = 13 I haven't been able to prove it yet. It involves recurrence relations, but seems quite complicated. I'll keep working at it.

Brian Charlesworth - 3 years, 7 months ago

Log in to reply

@Brian Charlesworth Also, would there be a way to calculate the exact value of PP? It's pretty easy to do for most recurrences, but do you think it can be done for this one?

Miles Koumouris - 3 years, 7 months ago

Log in to reply

@Miles Koumouris Also, regarding your approximation, since (from observation) a(n+1)a(n)<ϕ        n12\dfrac{a(n+1)}{a(n)}<\phi \;\; \forall \;\; n\geq 12, shouldn't P<k=112a(k)2k+a(13)213(1ϕ2)=14392048+158192(15+14)0.7122P<\sum_{k=1}^{12}\dfrac{a(k)}{2^k}+\dfrac{a(13)}{2^{13}\left(1-\frac{\phi }{2}\right)}=\dfrac{1439}{2048}+\dfrac{15}{8192\left(1-\frac{\sqrt{5}+1}{4}\right)}\approx 0.7122? I feel like I missed some information about the sequence when I made the assumption - would it be possible for you to show how you got your approximation?

Miles Koumouris - 3 years, 7 months ago

Log in to reply

@Miles Koumouris Sorry for the slow reply. Yes, you're calculation is correct; I should have triple-checked mine. :/

Brian Charlesworth - 3 years, 7 months ago

I believe that

P=1(n=1(1(12)n))=0.7112119049133975787211002780707692P=1-\left(\prod_{n=1}^{\infty }\left(1-\left(\dfrac{1}{2}\right)^n\right)\right)=0.7112119049133975787211002780707692\ldots

Based on the assumption that the ratio between successive terms of a(n)a(n) is monotonically increasing and upper-bounded by ϕ\phi for all n13n\geq 13, we can design an approximation technique for PP:

k=1M(a(k)2k)+a(M+1)2M+1(1(a(M+2)2a(M+1)))<P<k=1N(a(k)2k)+a(N+1)2N+1(1(ϕ2))\sum_{k=1}^M\left(\dfrac{a(k)}{2^k}\right)+\dfrac{a(M+1)}{2^{M+1}\left(1-\left(\frac{a(M+2)}{2a(M+1)}\right)\right)}<P<\sum_{k=1}^N\left(\dfrac{a(k)}{2^k}\right)+\dfrac{a(N+1)}{2^{N+1}\left(1-\left(\frac{\phi }{2}\right)\right)}

On the OEIS link, there are only 6060 terms provided, so setting M=58M=58 and N=59N=59 yields

34677104745555667007539167054875776756717555040447889408<P<715765590735+8199713396797771591152921504606846976\dfrac{3467710474555566700753916705}{4875776756717555040447889408}<P<\dfrac{71576559073\sqrt{5}+819971339679777159}{1152921504606846976}

which implies

0.7112119048<P<0.71121190510.7112119048\ldots <P<0.7112119051\ldots

and hence we can confirm PP rounded to 77 decimal places as 0.71121190.7112119. Upon searching this value, I found this link to an almost identical problem, with an answer of 0.71121190.7112119\ldots , found by evaluating the product

P=1(n=1(1(12)n))P=1-\left(\prod_{n=1}^{\infty }\left(1-\left(\dfrac{1}{2}\right)^n\right)\right)

So thanks to the much smarter people who actually did the work, came up with this, and proved it, I believe it wouldn't be too hard to use their solutions and verify that this is indeed the answer.

Miles Koumouris - 3 years, 7 months ago

Log in to reply

Great work! With several approaches to the problem yielding the same value I think that it's safe to conclude that the answer is verified. The formula you have for PP is surprisingly elegant. :)

Brian Charlesworth - 3 years, 7 months ago

Log in to reply

Thanks! Although I didn't really do any work - the others on stack exchange solved it, and you were the one that initially found the recurrence relation. I am surprised too by how simple the formula is... do you think PP can be expressed in closed form? Something tells me it can't, but I'm not sure how one would go about proving it.

Miles Koumouris - 3 years, 7 months ago

Log in to reply

@Miles Koumouris I would doubt that there is a closed form. It appears to be related to the Pentagonal Number Theorem.

Brian Charlesworth - 3 years, 7 months ago

Log in to reply

@Brian Charlesworth Thanks! I'll have a read.

Miles Koumouris - 3 years, 7 months ago

Dunno if this is right but....

Define PnP_n as the probability of winning starting right after you flipped a total of nn tails.

Pn=(12)n+1+(1(12)n+1)Pn+1P_n=\left(\dfrac{1}{2}\right)^{n+1}+\left(1-\left(\dfrac{1}{2}\right)^{n+1}\right)P_{n+1}

(2n1)Pn=2nPn11\Leftrightarrow (2^n-1)P_n=2^nP_{n-1}-1

Pretty obvious that P()=0.P(\infty)=0.

Aaaand I dunno if this is solvable x'D

Boi (보이) - 3 years, 7 months ago

0.5

Vivek Khatri - 3 years, 6 months ago

Since there is no way to lose, and it is possible to win given any previous sequence of tosses, the reasoning to solve this problem changes. The only way to 'not win' is to toss an endless sequence, since the game can only terminate in a win. There are clearly infinitely many of these endless games, but since the probability of one of these endless games is infinitesimally small ((12)N(\dfrac{1}{2})^N as NN\rightarrow \infty), it would be sufficient to prove that the infinite number of these endless games is not an exponential divergence. I'm not really sure how to address this...

Miles Koumouris - 3 years, 7 months ago
×

Problem Loading...

Note Loading...

Set Loading...