Ну я особо не полировал и не универсализировал скрипт.
На PARI/GP
Сначала заполняем два массива 100x100 (для
и
)
На ходу это считать долго, поэтому так делал.
Потом в четырехкратно вложенном цикле складываем каждый элемент одного массива с каждым элементом второго и смотрим не степень ли это.
Если степень, то печатаем порядковый номер решения, целую часть десятичного логарифма суммы, и само решение.
(Код)
Код:
Fg=matrix(100,100,i,j,sum(k=1,j,i^k));
Fk=matrix(100,100,i,j,i^j);
c=0;t=0;n=0;s=0;
for(i=2,100,for(j=2,100,for(k=2,100,for(l=2,100,s=Fg[k,l]+Fk[i,j];t=ispower(s,,&n);if(t>0,c=c+1;print(c," ",floor(log(s)/log(10))," digits: ",i,"^",j,"+F(",k,",",l,")=",n,"^",t))))))
Первые 10 строк вывода:
1 1 digits: 2^2+F(3,2)=2^4
2 2 digits: 2^2+F(5,4)=28^2
3 2 digits: 2^3+F(3,4)=2^7
4 1 digits: 2^3+F(7,2)=2^6
5 2 digits: 2^3+F(31,2)=10^3
6 1 digits: 2^4+F(4,2)=6^2
7 2 digits: 2^4+F(4,3)=10^2
8 3 digits: 2^4+F(4,6)=74^2
9 3 digits: 2^4+F(9,4)=86^2
10 2 digits: 2^4+F(15,2)=2^8Всего вывелось 449 строк. Ну ясно, что есть повторы в решениях потому что например
-- 18.10.2018, 13:10 --я посмотрю ещё варианты типа
и
.
Есть еще случаи
, например
:)