A345506 - OEIS

A345506

Numbers that are the sum of seven cubes in ten or more ways.

1704, 1711, 1774, 1800, 1837, 1863, 1889, 1893, 1926, 1938, 1963, 1982, 1989, 2008, 2015, 2019, 2045, 2052, 2053, 2059, 2078, 2097, 2106, 2113, 2143, 2161, 2169, 2171, 2176, 2197, 2204, 2217, 2223, 2224, 2227, 2230, 2234, 2241, 2250, 2260, 2266, 2267, 2276

(list; graph; refs; listen; history; text; internal format)

OFFSET

1,1

LINKS

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

EXAMPLE

1711 is a term because 1711 = 1^3 + 1^3 + 1^3 + 4^3 + 4^3 + 8^3 + 8^3 = 1^3 + 1^3 + 2^3 + 3^3 + 5^3 + 8^3 + 8^3 = 1^3 + 1^3 + 2^3 + 3^3 + 4^3 + 7^3 + 9^3 = 1^3 + 1^3 + 3^3 + 3^3 + 4^3 + 4^3 + 10^3 = 1^3 + 2^3 + 2^3 + 2^3 + 6^3 + 6^3 + 9^3 = 1^3 + 2^3 + 3^3 + 3^3 + 3^3 + 5^3 + 10^3 = 1^3 + 3^3 + 3^3 + 4^3 + 5^3 + 7^3 + 8^3 = 2^3 + 2^3 + 3^3 + 5^3 + 6^3 + 6^3 + 8^3 = 3^3 + 3^3 + 3^3 + 4^3 + 4^3 + 6^3 + 9^3 = 4^3 + 4^3 + 5^3 + 5^3 + 6^3 + 6^3 + 6^3.

PROG

(Python)

from itertools import combinations_with_replacement as cwr

from collections import defaultdict

keep = defaultdict(lambda: 0)

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

for pos in cwr(power_terms, 7):

tot = sum(pos)

keep[tot] += 1

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

for x in range(len(rets)):

print(rets[x])