OFFSET
1,1
COMMENTS
Square entries of A024364.
LINKS
FORMULA
a(n) = (2*u*v)^2, where u=sqrt(j/2) and v=sqrt(j+k) {for coprime pairs (j,k),j>k with odd k such that pairs (u,v),u<v are coprime with v odd}.
MAPLE
isA024364 := proc(an) local r::integer, s::integer ; for r from floor((an/4)^(1/2)) to floor((an/2)^(1/2)) do for s from r-1 to 1 by -2 do if 2*r*(r+s) = an and gcd(r, s) < 2 then RETURN(true) ; fi ; if 2*r*(r+s) < an then break ; fi ; od ; od : RETURN(false) ; end : isA120089 := proc(an) RETURN( issqr(an) and isA024364(an)) ; end: for n from 2 to 1200 do if isA120089(n^2) then printf("%d, ", n^2) ; fi ; od ; # R. J. Mathar, Jun 08 2006
MATHEMATICA
A078926[n_] := Sum[Boole[n < d^2 < 2n && CoprimeQ[d, n/d]], {d, Divisors[n/2^IntegerExponent[n, 2]]}];
Reap[For[k = 2, k <= 10^6, k += 2, If[A078926[k/2] > 0 && IntegerQ@Sqrt@k, Print[k]; Sow[k]]]][[2, 1]] (* Jean-François Alcover, Oct 25 2023 *)
CROSSREFS
KEYWORD
nonn
AUTHOR
Lekraj Beedassy, Jun 07 2006
EXTENSIONS
Corrected and extended by R. J. Mathar, Jun 08 2006
STATUS
approved