login
A345836
Numbers that are the sum of eight fourth powers in exactly four ways.
8
2933, 2948, 3013, 3173, 3188, 3557, 4148, 4163, 4213, 4293, 4388, 4453, 4643, 4772, 4837, 4883, 5012, 5123, 5188, 5203, 5268, 5333, 5363, 5378, 5398, 5428, 5538, 5573, 5603, 5618, 5668, 5733, 5748, 5858, 5923, 6052, 6163, 6227, 6292, 6548, 6578, 6628, 6693
OFFSET
1,1
COMMENTS
Differs from A345579 at term 10 because 4228 = 1^4 + 1^4 + 1^4 + 2^4 + 2^4 + 2^4 + 3^4 + 8^4 = 1^4 + 1^4 + 1^4 + 2^4 + 4^4 + 4^4 + 6^4 + 7^4 = 1^4 + 1^4 + 1^4 + 3^4 + 4^4 + 6^4 + 6^4 + 6^4 = 2^4 + 2^4 + 3^4 + 3^4 + 3^4 + 4^4 + 6^4 + 7^4 = 2^4 + 3^4 + 3^4 + 3^4 + 3^4 + 6^4 + 6^4 + 6^4.
LINKS
EXAMPLE
2948 is a term because 2948 = 1^4 + 1^4 + 1^4 + 2^4 + 2^4 + 4^4 + 4^4 + 7^4 = 1^4 + 1^4 + 1^4 + 2^4 + 3^4 + 4^4 + 6^4 + 6^4 = 2^4 + 2^4 + 2^4 + 3^4 + 3^4 + 3^4 + 4^4 + 7^4 = 2^4 + 2^4 + 3^4 + 3^4 + 3^4 + 3^4 + 6^4 + 6^4.
PROG
(Python)
from itertools import combinations_with_replacement as cwr
from collections import defaultdict
keep = defaultdict(lambda: 0)
power_terms = [x**4 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 == 4])
for x in range(len(rets)):
print(rets[x])
KEYWORD
nonn
AUTHOR
STATUS
approved