A004810 Numbers that are the sum of 10 positive 10th powers.
10, 1033, 2056, 3079, 4102, 5125, 6148, 7171, 8194, 9217, 10240, 59058, 60081, 61104, 62127, 63150, 64173, 65196, 66219, 67242, 68265, 118106, 119129, 120152, 121175, 122198, 123221, 124244, 125267, 126290, 177154, 178177, 179200, 180223, 181246, 182269, 183292, 184315
Offset: 1
Examples
From _David A. Corneth_, Aug 03 2020: (Start) 72332028 is in the sequence as 72332028 = 1^10 + 1^10 + 1^10 + 2^10 + 2^10 + 2^10 + 4^10 + 4^10 + 5^10 + 6^10. 243962883 is in the sequence as 243962883 = 1^10 + 1^10 + 1^10 + 2^10 + 4^10 + 4^10 + 6^10 + 6^10 + 6^10 + 6^10. 312998872 is in the sequence as 312998872 = 1^10 + 2^10 + 3^10 + 3^10 + 3^10 + 4^10 + 5^10 + 5^10 + 5^10 + 7^10. (End)
Links
- David A. Corneth, Table of n, a(n) for n = 1..10000 (first 1000 terms from T. D. Noe)
Programs
-
Mathematica
M = 1413602992; m = M^(1/10) // Ceiling; Reap[ For[a = 1, a <= m, a++, For[b = a, b <= m, b++, For[c = b, c <= m, c++, For[d = c, d <= m, d++, For[e = d, e <= m, e++, For[f = e, f <= m, f++, For[g = f, g <= m, g++, For[h = g, h <= m, h++, For[i = h, i <= m, i++, For[j = i, j <= m, j++, s = a^10 + b^10 + c^10 + d^10 + e^10 + f^10 + g^10 + h^10 + i^10 + j^10; If[s <= M, Sow[s]]]]]]]]]]]]][[2, 1]] // Union (* Jean-François Alcover, Dec 01 2020 *)