cp's OEIS Frontend

A211792 a(n) = Sum_{y=1..n} Sum_{x=1..n} floor((x^k + y^k)^(1/k)) with k = 3.

%I A211792 #20 Sep 12 2022 11:43:59
%S A211792 1,7,22,51,97,164,258,382,541,741,982,1271,1611,2008,2466,2986,3577,
%T A211792 4241,4982,5807,6715,7714,8808,10000,11297,12701,14217,15848,17600,
%U A211792 19477,21482,23620,25895,28313,30879,33592,36460,39487,42678,46036
%N A211792 a(n) = Sum_{y=1..n} Sum_{x=1..n} floor((x^k + y^k)^(1/k)) with k = 3.
%C A211792 Row 3 of A211798.
%H A211792 David A. Corneth, <a href="/A211792/b211792.txt">Table of n, a(n) for n = 1..10000</a>
%F A211792 a(n) = Sum_{y=1..n} Sum_{x=1..n} floor((x^3 + y^3)^(1/3)).
%F A211792 a(n) = a(n-1) + floor((2*n^3)^(1/3)) + 2*Sum_{i = 1..n-1} floor((n^3 + i^3)^(1/3)) for n >= 2 and a(1) = 1. - _David A. Corneth_, Sep 12 2022
%e A211792 For a(3) we get the floor() values (1+2+3) + (2+2+3) + (3+3+3) = 22.
%t A211792 f[x_, y_, k_] := f[x, y, k] = Floor[(x^k + y^k)^(1/k)]
%t A211792 t[k_, n_] := Sum[Sum[f[x, y, k], {x, 1, n}], {y, 1, n}]
%t A211792 Table[t[1, n], {n, 1, 45}]  (* 2*A002411 *)
%t A211792 Table[t[2, n], {n, 1, 45}]  (* A211791 *)
%t A211792 Table[t[3, n], {n, 1, 45}]  (* A211792 *)
%t A211792 TableForm[Table[t[k, n], {k, 1, 12},
%t A211792                  {n, 1, 16}]] (* A211798 *)
%t A211792 Flatten[Table[t[k, n - k + 1], {n, 1, 12}, {k, 1, n}]]
%o A211792 (PARI) first(n) = { res = vector(n); res[1] = 1; for(i = 2, n, i3 = i^3; s = sum(j = 1, i-1, sqrtnint(i3 + j^3, 3)); res[i] = res[i-1] + sqrtnint(2*i3, 3) + 2*s; ); res } \\ _David A. Corneth_, Sep 12 2022
%Y A211792 Cf. A211791, A211798.
%K A211792 nonn
%O A211792 1,2
%A A211792 _Clark Kimberling_, Apr 26 2012
%E A211792 Definition changed by _Georg Fischer_, Sep 10 2022