A345508 Numbers that are the sum of ten squares in one or more ways.

10, 13, 16, 18, 19, 21, 22, 24, 25, 26, 27, 28, 29, 30, 31, 32, 33, 34, 35, 36, 37, 38, 39, 40, 41, 42, 43, 44, 45, 46, 47, 48, 49, 50, 51, 52, 53, 54, 55, 56, 57, 58, 59, 60, 61, 62, 63, 64, 65, 66, 67, 68, 69, 70, 71, 72, 73, 74, 75, 76, 77, 78, 79, 80, 81

Offset

1,1

Links

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

Example

13 is a term because 13 = 1^2 + 1^2 + 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, 10):
    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

Cf. A003333, A345498, A345509.

Sequence in context: A135548 A102502 A309304 * A153045 A055984 A282108

Adjacent sequences: A345505 A345506 A345507 * A345509 A345510 A345511

Keyword

nonn

Author

David Consiglio, Jr., Jun 19 2021

Status

approved