login
A345827
Numbers that are the sum of seven fourth powers in exactly five ways.
8
6642, 6707, 6772, 6882, 6947, 7922, 7987, 8227, 8962, 9267, 9507, 9747, 10116, 10291, 10722, 10867, 10932, 10962, 11331, 11411, 11571, 12676, 12851, 12916, 13187, 13252, 13891, 13956, 14131, 14211, 14707, 14772, 14802, 14917, 14932, 14947, 15012, 15092, 15316
OFFSET
1,1
COMMENTS
Differs from A345571 at term 16 because 10787 = 1^4 + 1^4 + 1^4 + 6^4 + 6^4 + 8^4 + 8^4 = 1^4 + 2^4 + 2^4 + 2^4 + 3^4 + 8^4 + 9^4 = 1^4 + 2^4 + 4^4 + 4^4 + 6^4 + 7^4 + 9^4 = 1^4 + 3^4 + 4^4 + 6^4 + 6^4 + 6^4 + 9^4 = 2^4 + 2^4 + 3^4 + 3^4 + 7^4 + 8^4 + 8^4 = 3^4 + 4^4 + 4^4 + 6^4 + 7^4 + 7^4 + 8^4.
LINKS
EXAMPLE
6707 is a term because 6707 = 1^4 + 1^4 + 1^4 + 2^4 + 6^4 + 6^4 + 8^4 = 1^4 + 2^4 + 2^4 + 2^4 + 2^4 + 3^4 + 9^4 = 2^4 + 2^4 + 2^4 + 3^4 + 3^4 + 7^4 + 8^4 = 2^4 + 3^4 + 4^4 + 4^4 + 6^4 + 7^4 + 7^4 = 3^4 + 3^4 + 4^4 + 6^4 + 6^4 + 6^4 + 7^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, 7):
tot = sum(pos)
keep[tot] += 1
rets = sorted([k for k, v in keep.items() if v == 5])
for x in range(len(rets)):
print(rets[x])
KEYWORD
nonn
AUTHOR
STATUS
approved