OFFSET
1,1
COMMENTS
a(25) = 446185740, and a(n) for n >= 17 and n != 25 is larger, if it exists. - Amiram Eldar, Feb 07 2026
PROG
(PARI) A024360(n) = {my(f = factor(n), n2 = n^2, s, b, c, d, g); s = if(n == 1 || (n%4 == 2), 0, 2^(omega(f)-1)); for(k = 1, #f~, f[k, 2] *= 2); fordiv(f, d, g = n2\d; if(d <= g && (d+g) % 2 == 0, c = (d+g)\2; b = g-c; if(n < b && gcd(b, c) == 1, s--))); s; }
list(len) = {my(v = vector(len), c = 0, n = 1, i); while(c < len, i = A024360(n); if(i > 0 && i <= len && v[i] == 0, c++; v[i] = n); n++); v; } \\ Amiram Eldar, Feb 07 2026
CROSSREFS
KEYWORD
nonn,more
AUTHOR
Lorenzo Sauras Altuzarra, Feb 02 2026
EXTENSIONS
a(6)-a(10) from Michel Marcus, Feb 02 2026
a(11)-a(16) from Amiram Eldar, Feb 07 2026
STATUS
approved
