OFFSET

1,1

LINKS

Sean A. Irvine, Table of n, a(n) for n = 1..1000

EXAMPLE

11 is a term because 11 = 1^2 + 1^2 + 1^2 + 1^2 + 1^2 + 1^2 + 1^2 + 2^2.

PROG

(Python)

from itertools import combinations_with_replacement as cwr

from collections import defaultdict

keep = defaultdict(lambda: 0)

power_terms = [x**2 for x in range(1, 1000)]

for pos in cwr(power_terms, 8):

tot = sum(pos)

keep[tot] += 1

rets = sorted([k for k, v in keep.items() if v >= 1])

for x in range(len(rets)):

print(rets[x])

CROSSREFS

KEYWORD

nonn

AUTHOR

David Consiglio, Jr., Jun 19 2021

STATUS

approved