# Greetings from The On-Line Encyclopedia of Integer Sequences! http://oeis.org/
Search: id:a054347
Showing 1-1 of 1
%I A054347 #53 Aug 30 2022 02:12:30
%S A054347 0,1,4,8,14,22,31,42,54,68,84,101,120,141,163,187,212,239,268,298,330,
%T A054347 363,398,435,473,513,555,598,643,689,737,787,838,891,946,1002,1060,
%U A054347 1119,1180,1243,1307,1373,1440,1509,1580,1652,1726,1802,1879
%N A054347 Partial sums of A000201.
%C A054347 From _Michel Dekking_, Aug 19 2019: (Start)
%C A054347 Limit_{n->oo} a(n)/(n*(n+1)) = phi/2.
%C A054347 Proof: Let {alpha} be the fractional part of a real number alpha and let [alpha] = floor(alpha).
%C A054347 a(n) = [phi] + [2*phi] + ... + [n*phi] = phi + {phi} + 2*phi + {2*phi} + ... + n*phi + {n*phi} = n*(n+1)*phi/2 + [{phi} + {2*phi} + ... + {n*phi}].
%C A054347 When we divide by n*(n+1) this tends to phi/2, since the second term is bounded by n.
%C A054347 (End)
%H A054347 T. D. Noe, Table of n, a(n) for n = 0..10000
%H A054347 M. Griffiths, The Golden String, Zeckendorf Representations, and the Sum of a Series, Amer. Math. Monthly, 118 (2011), 497-507.
%F A054347 a(n) = floor(n*(n+1)/2*phi - n/2) + 0 or +1. - _Benoit Cloitre_, Oct 03 2003
%F A054347 a(n) = floor(n*(n+1)/2*phi - n/2) + 0, +1, or -1 (n = 7920, 18762, 18851, ...), or +2 (n = 12815, 15841, 30358, 30382, ...) if n < 2000000. - _Birkas Gyorgy_, May 06 2011
%t A054347 Accumulate[Table[Floor[GoldenRatio n], {n, 0, 30}]] (* _Birkas Gyorgy_, May 06 2011 *)
%o A054347 (PARI) for(n=0,50, print1(sum(k=0,n, floor(k*(1+sqrt(5))/2)), ", ")) \\ _G. C. Greubel_, Oct 06 2017
%o A054347 (Python)
%o A054347 from math import isqrt
%o A054347 from itertools import count, islice, accumulate
%o A054347 def A054347_gen(): # generator of terms
%o A054347 return accumulate(n+isqrt(5*n**2)>>1 for n in count(0))
%o A054347 A054347_list = list(islice(A054347_gen(),30)) # _Chai Wah Wu_, Aug 29 2022
%Y A054347 Cf. A000201.
%K A054347 nonn,easy
%O A054347 0,3
%A A054347 _N. J. A. Sloane_, May 06 2000
# Content is available under The OEIS End-User License Agreement: http://oeis.org/LICENSE