OFFSET
1,3
LINKS
Peter Kagey, Table of n, a(n) for n = 1..10000
EXAMPLE
For n = 14, A001481(14) = 25 = 0^2 + 5^2 = 3^2 + 4^2, so a(14) = max{0+5, 3+4} = 7.
PROG
(Python)
from itertools import count, islice
from sympy.solvers.diophantine.diophantine import diop_DN
from sympy import factorint
def A328804_gen(): # generator of terms
return map(lambda n: max((a+b for a, b in diop_DN(-1, n))), filter(lambda n:(lambda m:all(d&3!=3 or m[d]&1==0 for d in m))(factorint(n)), count(0)))
CROSSREFS
KEYWORD
nonn
AUTHOR
Peter Kagey, Oct 27 2019
STATUS
approved