comment on

I saw this same problem a while back myself. If I recall correctly, the answer was 50 before the odds would be better than 50-50 of 2 people in the room having the same birthday. The 1st person as 1/365th of a chance of being on a particular day, the odds the 2nd person having one other than the first is 364/365ths, the 3rd having a different birthday 363/365ths, and so forth.

As to modeling this situation, you would have to keep track of the day (potentially an array 0..364), and add people until you hit the same day for two of them. Then, you would need to simulate it multiple times (such as your $TRIALS=1024) to get an estimate.

Enjoy. (The code below might give you a starting place, but may require you to run considerably longer to get a better estimate.)

#!/usr/bin/perl -w

use strict;
use warnings;

srand();

my $maxday = 365;
my $TRIALS = 1024000;
my $average = 0;
my $sum = 0;
my $tenpercent = $TRIALS / 10;

for (my $count = 0; $count < $TRIALS; $count++) {
    $sum += &test_bdays($maxday);
    if (($count) and (!($count % $tenpercent)))    {
        $average = ($sum * 1.0) / $count;
        print("After $count trials, ", 
"the average number of people in the room was $average\n");
    }
}
$average = ($sum * 1.0) / $TRIALS;
print("After $TRIALS trials, ", "
the average number of people in the room was $average\n");

sub test_bdays {
    my ($mday) = @_;
    my @days = ();
    my $d = $mday;
    my $people = 0;
    while ($d) {
        $days[$d] = 0;
        $d--;
    }
    my $flag = 0;
    while (!($flag)) {
        $people++;
        my $index = int(rand() * $mday);
        $days[$index]++;
        $flag = ($days[$index] > 1);
    }
    return($people);
}
[download]

Update:

thraxil is correct-I bow both to his math skill and to his remembering the correct answer. An article at http://www.people.virginia.edu/~rjh9u/birthday.html discusses it more. It had made me wonder when the simulations I ran with the script above seemed to fall within the range between 20 and 25....

In reply to Re: Birthday Chances by atcroft
in thread Birthday Chances by YuckFoo

Are you posting in the right place? Check out Where do I post X? to know for sure.
Posts may use any of the Perl Monks Approved HTML tags. Currently these include the following:
<code> <a> <b> <big> <blockquote> <br /> <dd> <dl> <dt> <em> <font> <h1> <h2> <h3> <h4> <h5> <h6> <hr /> <i> <li> <nbsp> <ol> <p> <small> <strike> <strong> <sub> <sup> <table> <td> <th> <tr> <tt> <u> <ul>
Snippets of code should be wrapped in <code> tags not <pre> tags. In fact, <pre> tags should generally be avoided. If they must be used, extreme care should be taken to ensure that their contents do not have long lines (<70 chars), in order to prevent horizontal scrolling (and possible janitor intervention).
Want more info? How to link or How to display code and escape characters are good places to start.


go ahead... be a heretic
	PerlMonks