Почти решил. Полного доказательства пока нет.
Обозначим
![$$m+k=a,\,\,\,m-k=b$$ $$m+k=a,\,\,\,m-k=b$$](http://fx.ifz.ru/tex2.php?d=120&i=%24%24m%2Bk%3Da%2C%5C%2C%5C%2C%5C%2Cm-k%3Db%24%24)
.
Тогда решение
![$$(x,y)$$ $$(x,y)$$](http://fx.ifz.ru/tex2.php?d=120&i=%24%24%28x%2Cy%29%24%24)
выражается через
![$$(a,b),\,\,a,b>0$$ $$(a,b),\,\,a,b>0$$](http://fx.ifz.ru/tex2.php?d=120&i=%24%24%28a%2Cb%29%2C%5C%2C%5C%2Ca%2Cb%3E0%24%24)
таким образом:
![$$x=\frac{ab-1 \pm \sqrt{(a^2+1)(b^2+1)}}{2}$$ $$x=\frac{ab-1 \pm \sqrt{(a^2+1)(b^2+1)}}{2}$$](http://fx.ifz.ru/tex2.php?d=120&i=%24%24x%3D%5Cfrac%7Bab-1%20%5Cpm%20%5Csqrt%7B%28a%5E2%2B1%29%28b%5E2%2B1%29%7D%7D%7B2%7D%24%24)
![$$y=\frac{-ab-1 \pm \sqrt{(a^2+1)(b^2+1)}}{2}$$ $$y=\frac{-ab-1 \pm \sqrt{(a^2+1)(b^2+1)}}{2}$$](http://fx.ifz.ru/tex2.php?d=120&i=%24%24y%3D%5Cfrac%7B-ab-1%20%5Cpm%20%5Csqrt%7B%28a%5E2%2B1%29%28b%5E2%2B1%29%7D%7D%7B2%7D%24%24)
Тогда любая пара
![$$(a,b)$$ $$(a,b)$$](http://fx.ifz.ru/tex2.php?d=120&i=%24%24%28a%2Cb%29%24%24)
, для которых
![$$(a^2+1)(b^2+1)$$ $$(a^2+1)(b^2+1)$$](http://fx.ifz.ru/tex2.php?d=120&i=%24%24%28a%5E2%2B1%29%28b%5E2%2B1%29%24%24)
- полный квадрат, определяет целочисленное решение (полуцелого быть не может, достаточно рассмотреть остатки дискриминанта по модулю 4). Назовем пары
![$$(a,b)$$ $$(a,b)$$](http://fx.ifz.ru/tex2.php?d=120&i=%24%24%28a%2Cb%29%24%24)
, которым соответствуют целочисленные решения, хорошими, т.e. для них
![$$(a^2+1)(b^2+1)$$ $$(a^2+1)(b^2+1)$$](http://fx.ifz.ru/tex2.php?d=120&i=%24%24%28a%5E2%2B1%29%28b%5E2%2B1%29%24%24)
- полный квадрат.
Верно утверждение, что если
![$$(a,b)$$ $$(a,b)$$](http://fx.ifz.ru/tex2.php?d=120&i=%24%24%28a%2Cb%29%24%24)
и
![$$(a,c)$$ $$(a,c)$$](http://fx.ifz.ru/tex2.php?d=120&i=%24%24%28a%2Cc%29%24%24)
- хорошие пары, то
![$$(b,c)$$ $$(b,c)$$](http://fx.ifz.ru/tex2.php?d=120&i=%24%24%28b%2Cc%29%24%24)
- тоже хорошая пара (это несложно проверить). Далее, для любого
![$$a$$ $$a$$](http://fx.ifz.ru/tex2.php?d=120&i=%24%24a%24%24)
, существует как минимум одно число, составляющее c данным хорошую пару, например
![$$4a^3+3a$$ $$4a^3+3a$$](http://fx.ifz.ru/tex2.php?d=120&i=%24%244a%5E3%2B3a%24%24)
.
Получаем, все множество натуральных чисел разбивается на подмножества, в каждом из которых все числа друг c другом образуют хорошую пару. Написал программку и выяснил, что структура этих подмножеств такова: каждое является рекурсивной последовательностью, c первыми двумя элементами
![$$z_0=a,\,\, z_1=a(4a^2+3)$$ $$z_0=a,\,\, z_1=a(4a^2+3)$$](http://fx.ifz.ru/tex2.php?d=120&i=%24%24z_0%3Da%2C%5C%2C%5C%2C%20z_1%3Da%284a%5E2%2B3%29%24%24)
и зависимостью
![$$z_{n+1}=(4a^2+2)z_n-z_{n-1}$$ $$z_{n+1}=(4a^2+2)z_n-z_{n-1}$$](http://fx.ifz.ru/tex2.php?d=120&i=%24%24z_%7Bn%2B1%7D%3D%284a%5E2%2B2%29z_n-z_%7Bn-1%7D%24%24)
. Теперь эту запись надо бы поупрощать, т.e вывести сразу формулу
![$$z_n=f(a)$$ $$z_n=f(a)$$](http://fx.ifz.ru/tex2.php?d=120&i=%24%24z_n%3Df%28a%29%24%24)
, после этого скорее всего доказательство, что все числа в этой последовательности образуют c a хорошую пару, не составит труда.
Теперь для каждого натурального
![$$a$$ $$a$$](http://fx.ifz.ru/tex2.php?d=120&i=%24%24a%24%24)
можно найти свою последовательность, и им будет соответствовать свой набор решений. Правда, если некоторое
![$$d$$ $$d$$](http://fx.ifz.ru/tex2.php?d=120&i=%24%24d%24%24)
находится в последовательности, порождаемой
![$$a$$ $$a$$](http://fx.ifz.ru/tex2.php?d=120&i=%24%24a%24%24)
, то последовательность, порождаемая
![$$d$$ $$d$$](http://fx.ifz.ru/tex2.php?d=120&i=%24%24d%24%24)
, будет полностью содержаться в последовательности, порождаемой
![$$a$$ $$a$$](http://fx.ifz.ru/tex2.php?d=120&i=%24%24a%24%24)
, то есть решения системы будут повторяться, но это ничего, я думаю.
B общем, структура решений есть, осталось только полное и строгое доказательство.