Процесс можно представить по-другому. Очередной неудачник не уходит домой, а становится в конец очереди.
Кажется, Вы всё-таки не поняли условия. Цепь Маркова неоднородная тут, конечно, есть изначально - количество оставшихся шляп после каждого разбора образует ЦМ, толку-то от неё. В исходном сообщении всё, что верно - это только

, а последний абзац вообще никак. Всем хорошо известно, что вероятность никому не выбрать свои шляпы есть

, где

- субфакториал, а никакие не степени из-за зависимости

.
Доказывать проще всего по индукции. Пусть

- число шагов, которые потребовалось для полного разбора всех

шляп. Пусть для

известно, что

. Докажем, что

.
Пусть

- число шляп, которые ушли за первый шаг, если их изначально было

. Пусть распределение

описывается набором вероятностей

,

. Они выписываются, конечно, через те же субфакториалы, но мы без них обойдёмся. Достаточно того, что

.
Тогда имеет место следующая рекуррентная формула. Один раз разобрали

шляп. Тогда

(уже осталось на 1 шаг меньше) совпадает по распределению с величиной

Последние два слагаемых нулевые, но это неважно. В этой сумме

- это количества шагов после того, как осталось

шляп. Они независимы с результатом предыдущего извлечения шляп

. Их надо было бы, конечно, новыми буковками обозначить, потому как

"старая" от

зависит, а "новая" уже нет, но мне стало лениво ещё индексы заводить.
(используем индукционное предположение)

(перегруппируем слагаемые)

Отсюда

,

.