login
The OEIS is supported by the many generous donors to the OEIS Foundation.

 

Logo
Hints
(Greetings from The On-Line Encyclopedia of Integer Sequences!)
A345509 Numbers that are the sum of ten squares in two or more ways. 5
25, 28, 31, 33, 34, 36, 37, 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, 82, 83, 84, 85, 86, 87, 88, 89, 90, 91, 92, 93, 94, 95, 96 (list; graph; refs; listen; history; text; internal format)
OFFSET

1,1

LINKS

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

EXAMPLE

28 is a term because 28 = 1^2 + 1^2 + 1^2 + 1^2 + 1^2 + 1^2 + 1^2 + 1^2 + 2^2 + 4^2 = 1^2 + 1^2 + 1^2 + 1^2 + 2^2 + 2^2 + 2^2 + 2^2 + 2^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 >= 2])

    for x in range(len(rets)):

        print(rets[x])

CROSSREFS

Cf. A345499, A345508, A345510, A345550.

Sequence in context: A221617 A083806 A239521 * A295748 A254226 A195613

Adjacent sequences:  A345506 A345507 A345508 * A345510 A345511 A345512

KEYWORD

nonn

AUTHOR

David Consiglio, Jr., Jun 20 2021

STATUS

approved

Lookup | Welcome | Wiki | Register | Music | Plot 2 | Demos | Index | Browse | More | WebCam
Contribute new seq. or comment | Format | Style Sheet | Transforms | Superseeker | Recents
The OEIS Community | Maintained by The OEIS Foundation Inc.

License Agreements, Terms of Use, Privacy Policy. .

Last modified October 4 04:26 EDT 2022. Contains 357237 sequences. (Running on oeis4.)