|
%I #10 Aug 05 2021 15:29:12
%S 1296378,1371735,1409400,1614185,1824040,1885248,2016496,2101464,
%T 2302028,2305395,2542968,2562624,2851848,2889216,2974392,2988441,
%U 3185792,3380833,3681280,3689496,3706984,3775680,3906657,4109832,4123008,4142683,4422592,4525632,4783680
%N Numbers that are the sum of three third powers in six or more ways.
%H David Consiglio, Jr., <a href="/A345083/b345083.txt">Table of n, a(n) for n = 1..1000</a>
%e 1296378 is a term because 1296378 = 3^3 + 75^3 + 94^3 = 8^3 + 32^3 + 107^3 = 20^3 + 76^3 + 93^3 = 30^3 + 58^3 + 101^3 = 32^3 + 80^3 + 89^3 = 59^3 + 74^3 + 86^3.
%o (Python)
%o from itertools import combinations_with_replacement as cwr
%o from collections import defaultdict
%o keep = defaultdict(lambda: 0)
%o power_terms = [x**3 for x in range(1, 1000)]
%o for pos in cwr(power_terms, 3):
%o tot = sum(pos)
%o keep[tot] += 1
%o rets = sorted([k for k, v in keep.items() if v >= 6])
%o for x in range(len(rets)):
%o print(rets[x])
%Y Cf. A025334, A343967, A344647, A345084, A345086, A345148.
%K nonn
%O 1,1
%A _David Consiglio, Jr._, Jun 07 2021
|
