В данной задаче 4 неизвестных - это распределение входов

. На них есть одно дополнительное соотношение, так как сумма равна 1, плюс условие неотрицательности. Условные веротности переходов

заданы. Выходные вероятности

являются при этом линейными функциями входных, а совместные вероятности

- квадратичными функциями. Ваша задача при этом записывается в виде задачи на максимум при ограничениях типа равенств и неравенств. Для этих задач есть алгоритмы решения.
Но в данной задаче явно просматривается симметрия. Если я не ошибаюсь (нужно уточнить), то совместная информация есть выпуклый функционал по распределениям. А тогда из симметрии легко показать, что симметричные переменные должны быть равными. Подозреваю, что максимум достигнется на равномерном распределении.