|
|
A003346
|
|
Numbers that are the sum of 12 positive 4th powers.
|
|
45
|
|
|
12, 27, 42, 57, 72, 87, 92, 102, 107, 117, 122, 132, 137, 147, 152, 162, 167, 172, 177, 182, 187, 192, 197, 202, 212, 217, 227, 232, 242, 247, 252, 257, 262, 267, 277, 282, 292, 297, 307, 312, 322, 327, 332, 342, 347, 357, 362, 372, 377, 387, 392, 402, 407, 412, 417
(list;
graph;
refs;
listen;
history;
text;
internal format)
|
|
|
OFFSET
|
1,1
|
|
COMMENTS
|
a(88) = 636 = 5^4 + 11 and a(91) = 651 = 5^4 + 2^4 + 10 are the first two terms not congruent to 2 or 7 (mod 10). - M. F. Hasler, Aug 03 2020
|
|
LINKS
|
|
|
EXAMPLE
|
3740 is in the sequence as 3740 = 1^4 + 1^4 + 1^4 + 1^4 + 1^4 + 1^4 + 1^4 + 1^4 + 3^4 + 5^4 + 5^4 + 7^4.
4690 is in the sequence as 4690 = 2^4 + 2^4 + 2^4 + 2^4 + 2^4 + 4^4 + 4^4 + 4^4 + 5^4 + 5^4 + 6^4 + 6^4.
7193 is in the sequence as 7193 = 2^4 + 4^4 + 5^4 + 5^4 + 5^4 + 5^4 + 5^4 + 5^4 + 5^4 + 5^4 + 5^4 + 6^4. (End)
|
|
PROG
|
(PARI) (A003346_upto(N, k=12, m=4)=[i|i<-[1..#N=sum(n=1, sqrtnint(N, m), 'x^n^m, O('x^N))^k], polcoef(N, i)])(500) \\ 2nd & 3rd optional arg allow to get other sequences of this group. See A003333 for alternate code. - M. F. Hasler, Aug 03 2020
(Python)
from itertools import count, takewhile, combinations_with_replacement as mc
def aupto(limit):
qd = takewhile(lambda x: x <= limit, (k**4 for k in count(1)))
ss = set(sum(c) for c in mc(qd, 12))
return sorted(s for s in ss if s <= limit)
|
|
CROSSREFS
|
Other numbers that are the sum of k positive m-th powers:
|
|
KEYWORD
|
nonn,easy
|
|
AUTHOR
|
|
|
STATUS
|
approved
|
|
|
|