Numerical probability and Monte Carlo solution for any k shared birthdays for year with any d days and any group size n