It really is (n/(n+1))^n for n sided dice, approximating to 1/e as n increases. Confirmed it by simulating in R program. This is the nice way to see the e value with probabilistic simulation. I remember that the probability of derangement approaches 1/e. But This one is really super cool!