По-видимому, альтернативный метод (вместо использования расширенного алгоритма Евклида) получения экспоненты для расшифровки - это сделать d = e ** (phi (phi (n)) - 1) mod (phi ( п)). Почему это работает?