|
%I #31 Jun 29 2024 10:57:17
%S 0,7,91,637,1729,31213,12103,405769,53599,157339,593047
%N a(n) is the smallest nonnegative integer k where there are exactly n solutions to x^2 + x*y + y^2 = k with 0 < x < y.
%C a(n) is the smallest nonnegative k such that A374092(k) = n.
%C a(11) > 10^8. - _Robert Israel_, Jun 28 2024
%F a(n) <= 13 * 7^(n-1).
%p N:= 10^6:
%p V:= Array(0..N):
%p for x from 1 to floor(sqrt(N/3)) do
%p for y from x+1 do
%p v:= x^2 + x*y + y^2;
%p if v > N then break fi;
%p V[v]:= V[v]+1;
%p od od:
%p W:= Array(0..10);
%p for i from 1 to N while count < 11 do
%p v:= V[i];
%p if W[v] = 0 then W[v]:= i; count:= count+1 fi
%p od:
%p 0, seq(W[i],i=1..10); # _Robert Israel_, Jun 28 2024
%o (Python)
%o from itertools import count
%o from sympy.abc import x,y
%o from sympy.solvers.diophantine.diophantine import diop_quadratic
%o def A374094(n): return next(m for m in count(0) if sum(1 for d in diop_quadratic(x*(x+y)+y**2-m) if 0<d[0]<d[1]) == n) # _Chai Wah Wu_, Jun 28 2024
%Y Cf. A093195, A374095.
%Y Cf. A374090, A374092.
%K nonn,more
%O 0,2
%A _Seiichi Manyama_, Jun 28 2024
| 