cp's OEIS Frontend

This is a front-end for the Online Encyclopedia of Integer Sequences, made by Christian Perfect. The idea is to provide OEIS entries in non-ancient HTML, and then to think about how they're presented visually. The source code is on GitHub.

A014088 Minimal number of people to give a 50% probability of having at least n coincident birthdays in one year.

Original entry on oeis.org

1, 23, 88, 187, 313, 460, 623, 798, 985, 1181, 1385, 1596, 1813, 2035, 2263, 2494, 2730, 2970, 3213, 3459, 3707, 3959, 4213, 4470, 4728, 4989, 5252, 5516, 5783, 6051, 6320, 6592, 6864, 7138, 7413, 7690, 7968, 8247, 8527, 8808, 9090, 9373, 9657, 9942, 10228
Offset: 1

Views

Author

Keywords

Crossrefs

Cf. A033810 (2 people on n days), A225852 (3 people on n days), A225871 (4 people on n days).

Programs

  • Mathematica
    q[1][n_, d_] := q[1][n, d] = d!/((d-n)!*d^n) // N; q[k_][n_, d_] := q[k][n, d] = Sum[ n!*d!/(d^(i* k)*i!*(k!)^i*(n-i*k)!*(d-i)!)*Sum[ q[j][n-i*k, d-i]*(d-i)^(n-i* k)/d^(n-i*k), {j, 1, k-1}], {i, 1, Floor[n/k]}] // N; p[k_][n_, d_] := 1 - Sum[q[i][n, d], {i, 1, k-1}]; a[1] = 1; a[k_] := a[k] = For[n = a[k-1], True, n++, If[p[k][n, 365] >= 1/2, Return[n]]]; Table[ Print["a(", k, ") = ", a[k]]; a[k], {k, 1, 15}] (* Jean-François Alcover, Jun 12 2013, after Eric W. Weisstein *)

Extensions

Broken links corrected by Steven Finch, Jan 27 2009
a(16)-a(45) from Hiroaki Yamanouchi, Mar 19 2015