SUMMARY
Here's a problem from Harvard's Stats 110 course.
For a group of 7 people, find the probability that all 4 seasons (winter, spring, summer, fall) occur at least once each among their birthdays, assuming that all seasons are equally likely.
The problem and its solution can be found here.
The hard part for me comes when calculating the probability that no people out of seven were born in winter. I don't understand why "order matters" here. (I understand the inclusion-exclusion part of the main problem, just not the "order matters" part.)
WHAT I TRIED
As one part of this problem, I need to calculate the probability that none of the seven people has a birthday in the winter. I tried to calculate this by considering the people as indistinguishable and lumping them into 4 different categories: birthdays in winter, spring, fall, or summer. Thinking like this, I used the "stars and bars" formula to calculate the probability that none of the birthdays fell in the winter category. This can be calculated as
$$P(A) = \dbinom{7 + 3 - 1}{3} / \dbinom{7+4-1}{4} $$
WHAT HARVARD DID
Harvard says this is wrong and that since order matters, the calculation is (the much easier)
$$P(A) = (3/4)^7$$
It seems like it could go either way, and I can actually model either in R. So what am I doing wrong?