OFFSET
1,1
COMMENTS
All terms are even.
EXAMPLE
24 is a term: j*k*(j+k) = 24^2 for j=2, k=16.
MATHEMATICA
Select[2*Range@500,
Length@Select[Table[(Sqrt[b^2 + 4 #^2/b] - b)/2, {b, #}], IntegerQ] >
0 &]
Select[Union@
Flatten@Table[Sqrt[a*b (a + b)], {a, 1, 80}, {b, a, 500}],
IntegerQ[#] && # < 1000 &]
PROG
(Python)
from itertools import count, islice
from sympy import integer_nthroot, divisors
def A363052_gen(startvalue=1): # generator of terms >= startvalue
for m in count(max(startvalue, 1)):
for k in divisors(m**2, generator=True):
p, q = integer_nthroot(k**4+(k*m**2<<2), 2)
if q:
a, b = divmod(p-k**2, k<<1)
if a > 0 and not b:
yield m
break
CROSSREFS
KEYWORD
nonn,easy
AUTHOR
Zhining Yang, May 15 2023
STATUS
approved