Sorry, I made a calculator typo yesterday. Yes, what you did works. It gets .1848 and 1-that=.8152
Now I have to figure out why it works...
Quote:
(1 - 4x4/48)*(1 - 4x4/47)*(1 - 4x4/46) x (1-4x4/45)
4×4/48 is the overestimated chance one player has an ace, overcounting the chance of multiple aces.
1-that = underestimated chance they lack an ace
Then you did the same thing for the other 3 players, then multiplied the 4 underestimated chances of each player lacking aces. And for some reason you only decreased the denominators by 1, even though ostensibly 4 cards had been removed each time.
But somehow the combination of those wonky things equals a working formula, and it beats the hell out of me. I love when this happens! Such a routine problem and yet I'm about to learn something.
It looks like if we expand the terms, an inclusion-exclusion formula will emerge.
1 - 16(1/45 + 1/46 + 1/47 + 1/48) +
256[1/(46*45) + 1/(47*45) + 1/(48*45) + 1/(47*46) + 1/(48*46) + 1/(48*47)] -
4096[1/(47*46*45) + 1/(48*46*45) + 1/(48*47*45) + 1/(48P3)] +
2^16 / 48P4
Sanity check: there should be 2^4 terms in the alternating sum and there are.
Great but that's not helping me much right now. I'll come back to this later.
OP, I'm glad you posted this! If you have a reason that this works then I'm all ears!