%I #87 Sep 17 2023 21:40:41
%S 2,2,3,3,3,4,4,4,4,5,5,5,5,5,5,5,6,6,6,6,6,6,6,7,7,7,7,7,7,7,7,7,8,8,
%T 8,8,8,8,8,8,8,8,9,9,9,9,9,9,9,9,9,9,9,9,10,10,10,10,10,10,10,10,10,
%U 10,10,10,10,10,11,11,11,11,11,11,11,11,11,11,11,11,11,11,12,12,12,12
%N Number of people needed so that probability of at least two sharing a birthday out of n possible days is at least 50%.
%C a(365) = 23 is the solution to the Birthday Problem.
%H T. D. Noe and Charles R Greathouse IV, <a href="/A033810/b033810.txt">Table of n, a(n) for n = 1..10000</a> (366 terms from T. D. Noe)
%H S. E. Ahmed and R. J. McIntosh, <a href="https://cms.math.ca/crux/v26/n3/page151-155.pdf">An Asymptotic Approximation for the Birthday Problem</a>, Crux Mathematicorum 26(3) 151-5 2000 CMS.
%H D. Brink, <a href="http://dx.doi.org/10.1007/s11139-011-9343-9">A (probably) exact solution to the Birthday Problem</a>, Ramanujan Journal, June 2012, Volume 28, Issue 2, pp. 223-238. - From _N. J. A. Sloane_, Oct 08 2012
%H Mathforum, <a href="http://mathforum.org/dr.math/faq/faq.birthdayprob.html">The Birthday Problem</a>.
%H Eric Weisstein's World of Mathematics, <a href="http://mathworld.wolfram.com/BirthdayProblem.html">Birthday Problem</a>.
%F a(n) = ceiling(sqrt(2*n*log(2)) + (3 - 2*log(2))/6 + (9 - 4*log(2)^2) / (72*sqrt(2*n*log(2))) - 2*log(2)^2/(135*n)) for all n up to 10^18. It is conjectured that this formula holds for all n.
%t lst = {}; s = 1; Do[Do[If[Product[(n - i + 1)/n, {i, j}] <= 1/2, If[j > s, s = j]; AppendTo[lst, j]; Break[]], {j, s, s + 1}], {n, 86}]; lst (* _Arkadiusz Wesolowski_, Apr 29 2012 *)
%t A033810[n_] := Catch@Do[If[1/2 >= n!/(n - m)!/n^m, Throw[m]], {m, 2, Infinity}]; Array[A033810, 86] (* _JungHwan Min_, Mar 27 2017 *)
%Y Essentially the same as A088141. See also A182008, A182009, A182010.
%Y Cf. A014088 (n people on 365 days), A225852 (3 on n days), A225871 (4 people on n days).
%K nonn,easy,nice
%O 1,1
%A _Eric W. Weisstein_