cp's OEIS Frontend

A363497 a(n) = Sum_{k=0..n} floor(sqrt(k))^3.

%I A363497 #45 Jun 15 2023 07:44:10
%S A363497 0,1,2,3,11,19,27,35,43,70,97,124,151,178,205,232,296,360,424,488,552,
%T A363497 616,680,744,808,933,1058,1183,1308,1433,1558,1683,1808,1933,2058,
%U A363497 2183,2399,2615,2831,3047,3263,3479,3695,3911,4127,4343,4559,4775,4991,5334
%N A363497 a(n) = Sum_{k=0..n} floor(sqrt(k))^3.
%C A363497 Partial sums of the third powers of the terms of A000196.
%H A363497 Karl-Heinz Hofmann, <a href="/A363497/b363497.txt">Table of n, a(n) for n = 0..10000</a>
%F A363497 a(n) = (n+1)*m^3 - (1/60)*m*(m+1)*(3*m-1)*(12*m^2+7*m-4), where m = floor(sqrt(n)).
%t A363497 Table[(n + 1) #^3 - (1/60) # (# + 1) (3 # - 1) (12 #^2 + 7 # - 4) &[Floor@ Sqrt[n]], {n, 0, 50}] (* _Michael De Vlieger_, Jun 10 2023 *)
%o A363497 (PARI) a(n) = sum(k=0, n, sqrtint(k)^3); \\ _Michel Marcus_, Jun 06 2023
%o A363497 (Python)
%o A363497 from math import isqrt
%o A363497 A363497 = [0]
%o A363497 for n in range(1,50): A363497.append(A363497[-1] + isqrt(n)**3)
%o A363497 print(A363497) # _Karl-Heinz Hofmann_, Jun 14 2023
%o A363497 (Python)
%o A363497 from math import isqrt
%o A363497 def A363497(n):return (m:=isqrt(n))**3*(n+1)-(m*(m+1)*(3*m-1)*(12*m**2+7*m-4))//60
%o A363497 # _Karl-Heinz Hofmann_, Jun 14 2023
%Y A363497 Sums of powers of A000196: A022554 (1st), A174060 (2nd), this sequence (3rd), A363498 (4th), A363499 (5th).
%K A363497 nonn,easy
%O A363497 0,3
%A A363497 _Hans J. H. Tuenter_, Jun 05 2023