%I #31 Aug 05 2023 22:01:14
%S 7,9,11,13,14,17,17,19,20,25,23,29,26,27,29,37,31,40,34,35,38,46,39,
%T 41,44,43,44,54,47,58,49,51,56,53,54,67,62,59,59,70,62,73,64,65,74,78,
%U 69,71,71,75,74,86,76,77,79,83,92,93,83,103
%N a(n) > n is the smallest integer such that there exist integers n < c <= d < a(n) satisfying n^2 + a(n)^2 = c^2 + d^2.
%C n^2 + a(n)^2 belongs to A007692.
%C The identity n^2 + (2*n + 5)^2 = (n+4)^2 + (2*n + 3)^2 shows that a(n) <= 2*n + 5. The last case when the equality holds is n = 16.
%C a(n) = a(n+1) has infinitely many solutions. This holds, in particular, when n = (u*v + u + v - 1) * (u*v - 2)/2 - 1 for positive integers u, v satisfying v+2 <= u <= 6*v - 3.
%C a(n-1) = a(n) = a(n+1) holds for n = (3*v^2 + 5*v + 1) * (6*v^2 + 3*v - 2), v >= 3.
%H Giedrius Alkauskas, <a href="https://web.vu.lt/mif/g.alkauskas/math/squares.pdf">On the function which is defined by the minimal solution to n^2 + f(n)^2 = a^2 + b^2, a, b, f(n) > n</a>.
%e a(10) = 25, since 10^2 + 25^2 = 14^2 + 23^2, and no integers b, c, d exist satisfying 10 < c <= d < b < 25 and 10^2 + b^2 = c^2 + d^2.
%p a :=proc(n::integer) local found::boolean; local N, SQ, i;
%p found:=false; N:=n+1; SQ:={};
%p while not found do SQ:=SQ union {N^2}; N:=N+1;
%p for i from n+1 to N-1 do
%p if evalb(N^2+n^2-i^2 in SQ) then found:=true; end if;
%p end do; end do; N end proc;
%Y Cf. A360619, A000404, A007692.
%K nonn
%O 1,1
%A _Giedrius Alkauskas_, Feb 21 2023
|