Да, читаю "Вероятность на алгебраических структурах (Гренадер)" и там характеристические функции - основной аппарат. Вопрос, как представить задачу в алгебраических, групповых или терминах колец?..
Пока в голову приходит лишь то, что P(n) будут стоять на диагонали в матрице, которая получается как произведение матриц, одна из которых: последовательно выписанные значения

со сдвигом вправо, на матрицу: последовательно выписанные значения

после сдвига вправо (остальные строки сдвигаются последовательно вправо). Причем, в остальных элементах матрицы будут выражения вида

, где i и j только одновременно < и > чем n (если одновременно равны, то будет

).