OFFSET
1,1
COMMENTS
This sequence is infinite. All the positive numbers of the form A141046(m) - 1 = 4*m^4 -1 are terms (De Koninck et al., 2024, section 6).
LINKS
Amiram Eldar, Table of n, a(n) for n = 1..5242 (terms below 10^10)
Jean-Marie De Koninck, A. Arthur Bonkli Razafindrasoanaivolala, and Hans Schmidt Ramiliarimanana, Integers with a sum of co-divisors yielding a square, Research in Number Theory, Vol. 10, No. 2 (2024), Article 30; author's copy.
MATHEMATICA
q[k_] := q[k] = If[IntegerQ[Sqrt[k]], IntegerQ[Sqrt[2*Sqrt[k]]], Module[{d = Divisors[k], nh}, nh = Length[d]/2; IntegerQ[Sqrt[d[[nh]] + d[[nh + 1]]]]]]; Select[Range[150000], q[#] && q[#+1] &]
PROG
(PARI) is1(k) = if(issquare(k), issquare(2 * sqrtint(k)), my(d = divisors(k), nh = #d/2); issquare(d[nh] + d[nh + 1]));
lista(kmax) = {my(q1 = is1(1), q2); for(k = 2, kmax, q2 = is1(k); if(q1 && q2, print1(k-1, ", ")); q1 = q2); }
CROSSREFS
KEYWORD
nonn
AUTHOR
Amiram Eldar, Nov 05 2024
STATUS
approved