Smoker has two matchboxes

А smoker has two matchboxes in his pocket, each containing $n$ mathes. When he wants to smoke, he selects randomly one of the boxes, and takes a match from it. At some point he finds out, that one of the boxes is empty. What is the probability that in the other box there are exactly $r$ matches left (for $r\in\{0,1,\ldots,n\}$). What about the expected number of the matches in the other box.
 
If $X$ is the random variable "number of matches in the other box", then one directly finds that $$\text{P}(X=r)={2n-r \choose n-r}\frac{1}{2^{2n-r}}$$
 Now we will evaluate $$\text{E}X=\sum_{r=0}^{n}r{2n-r \choose n-r}\frac{1}{2^{2n-r}}$$
 The sum of all probabilities being $1$ means here that $$\sum_{r=0}^{n}{2n-r \choose n-r}\frac{1}{2^{2n-r}}=1 $$
Change the variables $r=n-d$ in the last and and we obtain $$\sum_{d=0}^{n}{n+d \choose d}\frac{1}{2^{n+d}}=1 \ \ \ \ \ \ \  (*)$$ Make the same change in the sum for the expectation $$\text{E}X=\sum_{d=0}^{n}(n-d){n+d \choose d}\frac{1}{2^{n+d}}=n-\sum_{d=0}^{n}d{n+d \choose d}\frac{1}{2^{n+d}}$$
It remains to calculate the last sum which could be rewritten in the form $$\frac{1}{2^n n!}\sum_{d=0}^{n}d\frac{(n+d)!}{2^{d}d!}$$ First of all rewrite $(*)$ in the form
$$\sum_{d=0}^{n}\frac{(n+d)!}{2^d d!}=2^n n!$$
Write the same for $n+1$
$$\sum_{d=0}^{n+1}\frac{(n+d+1)!}{2^d d!}=2^{n+1} (n+1)!$$ and manipulating it we obtain
$$\sum_{d=0}^{n}(n+d+1)\frac{(n+d)!}{2^d d!}+\frac{(2n+2)!}{2^{n+1}(n+1)!}=\sum_{d=0}^{n}d\frac{(n+d)!}{2^d d!}+(n+1)\sum_{d=0}^{n}\frac{(n+d)!}{2^d d!}+\frac{(2n+2)!}{2^{n+1}(n+1)!}=\\= \sum_{d=0}^{n}d\frac{(n+d)!}{2^d d!}+(n+1)2^n n!+\frac{(2n+2)!}{2^{n+1}(n+1)!}=2^{n+1} (n+1)!$$
From the last we find $$ \sum_{d=0}^{n}d{n+d \choose d}\frac{1}{2^{n+d}}=\frac{1}{2^n n!}\left(2^{n+1} (n+1)!-2^n (n+1)!-\frac{(2n+2)!}{2^{n+1}(n+1)!}\right)$$
hence $$\text{E}X=n-\frac{1}{2^n n!}\left(2^n (n+1)!-\frac{(2n+2)!}{2^{n+1}(n+1)!}\right)=n-(n+1)+\frac{(2n+2)!}{2^{2n+1}(n+1)!n!}=\\=-1+\frac{n+1}{4^n}{2n+1\choose n}$$
 Using the Stirling approximation one finds that $$\text{E}X \sim \frac{ 2\sqrt{n}}{\sqrt{\pi}}$$

Comments. One can similarly find the variance of $X$ - just write out $(*)$ for $n+2$ so that squares of $d$ appear.

In essence the problem breaks down to evaluating $$ \sum_{d=0}^{n}d{n+d \choose d}\frac{1}{2^{n+d}}$$ Observe that this is quite easier if the sum was infinite - then this could be obtained from differentiating the series of $(1-s)^{-(n+1)}$ (up to some translation).
Another interesting piece of the problem is the identity $(*)$
$$\sum_{d=0}^{n}{n+d \choose d}\frac{1}{2^{n+d}}=1$$
From probabilistic point of view this is trivial. Is there any other way to evaluate it? This could lead to a different tackling of the previous sum (not using $(*)$). One more thing to note - if the sum in $(*)$ was infinite (it is easier to find as well) it is 2, and up to $n$ it is 1.
Something that could be related is Negative binomial distribution.

Addendum (11/17/2018) Probabilistic expanation of the identity $$ \sum_{d=0}^{n}{n+d \choose d}\frac{1}{2^{n+d}}=\frac{1}{2}\sum_{d=0}^{\infty}{n+d \choose d}\frac{1}{2^{n+d}}.$$ (Based on ideas of Emilian Rogachev)
Consider we perform a walk on the integer points of the Cartesian plane, being allowed only moves up and right. If we count the probability of going through the line $x=n+1$ before we go through $y=n+1$ it is $$ \sum_{d=0}^{n}{n+d \choose d}\frac{1}{2^{n+d+1}}.$$ The probability that we will ever go through $x=n+1$ is $1$ and could be written as $$ \sum_{d=0}^{\infty}{n+d \choose d}\frac{1}{2^{n+d+1}}.$$ However, since everything is symmetric, we see that going through $x=n+1$ before $y=n+1$ is equiprobable with going through $y=n+1$ before $x=n+1$ and since both sum up to $1$, each of those probabilities is $1/2$. So $$\sum_{d=0}^{n}{n+d \choose d}\frac{1}{2^{n+d+1}}=\frac{1}{2}\sum_{d=0}^{\infty}{n+d \choose d}\frac{1}{2^{n+d+1}}.$$

Comments

Popular posts from this blog

Basel problem type sum, proposed by prof. Skordev

Modification on a sequence from VJIMC

A problem proposed by prof. Babev