Ты думал в правильном направлении, и я ему поддался..
2013х+1=а²
2013х=(а-1)(а+1)
частные решения х очевидны:
х=2015, х=2011
должны подойти.
попробую сформулировать общий вид:
(a+1)z=x
(2013z+2)z=x
аналогично с обратным
получили:
х=(2013*к-2)*к, где к∈N
x=(2013*n+2)*n, где n∈N
решения.. то есть у - соответствующее тоже будет натуральным