|
%I #14 Jun 24 2018 14:05:20
%S 1,4,8,14,23,33,48,64,85,110,138,174,219,268,323,387,453,531,612,703,
%T 803,908,1028,1149,1285,1429,1582,1751,1922,2112,2308,2518,2743,2974,
%U 3227,3483,3759,4048,4348,4672,4997,5348,5709,6087,6487,6893,7328,7769,8234
%N Partial sums of union of squares and triangular numbers.
%C a(n) = Sum (A005214(k): 1 <= k <= n).
%H Reinhard Zumkeller, <a href="/A193711/b193711.txt">Table of n, a(n) for n = 1..10000</a>
%t Module[{upto=500,sq,tr},sq=Range[Floor[Sqrt[upto]]]^2;tr=Accumulate[ Range[ Floor[ (Sqrt[1+8*upto]-1)/2]]];Accumulate[Union[Join[tr,sq]]]] (* _Harvey P. Dale_, Jun 24 2018 *)
%o (Haskell)
%o a193711 n = a193711_list !! (n-1)
%o a193711_list = scanl1 (+) a005214_list :: [Integer]
%Y Cf. A000217, A000290, A000292, A000330.
%K nonn
%O 1,2
%A _Reinhard Zumkeller_, Aug 03 2011
%E Program and b-file corrected by _Reinhard Zumkeller_, Aug 04 2011
|
