OFFSET
1,2
COMMENTS
Sorted distinct values of ({A081961} - 1)/4.
LINKS
Karl-Heinz Hofmann, Table of n, a(n) for n = 1..10000
PROG
(Python)
from math import gcd, isqrt
test_all_k_upto = 193
A375750, x, limit = set(), 2, test_all_k_upto * 4 + 1
while x**2 + (lowY := isqrt(2*x**2)-x)**2 < limit:
for y in range(min(x-1, (yy:=isqrt(limit - x**2))-(yy%2 == x%2)), lowY, -2):
if gcd(x, y) == 1: A375750.add(((x**2 + y**2) - 1) // 4)
x += 1
(PARI) \\ Uses function is_a376208 from A376208
is_a376208(n, 1)
CROSSREFS
KEYWORD
nonn
AUTHOR
Hugo Pfoertner, Sep 13 2024
STATUS
approved