DiscretePlot[StirlingS2[n, 365]*365!/365^n, {n, 1000, 5000}]