Ок. Предлагаю тривиальные случаи

и

не рассматривать.
Т.е. условие тогда начинается так: "Пусть для простого

и натуральных

выполняется..."
Подсказка. Простое

образует две пары

=
(127,889),(127,5334).
Примеры (на pari/gp; в Вольфраме есть PowerMod; наверное, в любых CAS есть возведение в степень по модулю):
Код:
? j=0;a=2;x=a^a+j*889;Mod(x,127)^x
%1 = Mod(2, 127)
?
? j=777;a=2;x=a^a+j*889;Mod(x,127)^x
%2 = Mod(2, 127)
?
? j=787878;a=2;x=a^a+j*889;Mod(x,127)^x
%3 = Mod(2, 127)
?
?
? j=0;a=5;x=a^a+j*5334;Mod(x,127)^x
%4 = Mod(5, 127)
?
? j=888;a=5;x=a^a+j*5334;Mod(x,127)^x
%5 = Mod(5, 127)
?
? j=898989;a=5;x=a^a+j*5334;Mod(x,127)^x
%6 = Mod(5, 127)