login

Reminder: The OEIS is hiring a new managing editor, and the application deadline is January 26.

A345540
Numbers that are the sum of eight cubes in ten or more ways.
7
1185, 1243, 1262, 1288, 1295, 1299, 1386, 1393, 1397, 1400, 1412, 1419, 1423, 1448, 1449, 1451, 1458, 1460, 1464, 1467, 1475, 1477, 1497, 1501, 1503, 1504, 1505, 1512, 1513, 1514, 1516, 1521, 1523, 1539, 1540, 1541, 1542, 1553, 1558, 1559, 1560, 1565, 1566
OFFSET
1,1
LINKS
EXAMPLE
1243 is a term because 1243 = 1^3 + 1^3 + 1^3 + 2^3 + 2^3 + 2^3 + 5^3 + 9^3 = 1^3 + 1^3 + 1^3 + 3^3 + 3^3 + 4^3 + 7^3 + 7^3 = 1^3 + 1^3 + 2^3 + 2^3 + 3^3 + 6^3 + 6^3 + 7^3 = 1^3 + 1^3 + 2^3 + 2^3 + 3^3 + 5^3 + 5^3 + 8^3 = 1^3 + 1^3 + 4^3 + 4^3 + 5^3 + 5^3 + 5^3 + 6^3 = 1^3 + 2^3 + 3^3 + 5^3 + 5^3 + 5^3 + 5^3 + 6^3 = 2^3 + 2^3 + 2^3 + 3^3 + 3^3 + 3^3 + 3^3 + 9^3 = 2^3 + 3^3 + 3^3 + 3^3 + 4^3 + 6^3 + 6^3 + 6^3 = 2^3 + 3^3 + 3^3 + 3^3 + 3^3 + 4^3 + 4^3 + 8^3 = 3^3 + 3^3 + 3^3 + 3^3 + 3^3 + 5^3 + 6^3 + 7^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, 8):
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])
CROSSREFS
KEYWORD
nonn
AUTHOR
STATUS
approved