OFFSET
0,2
COMMENTS
Partial sums of A001950.
A001950 is the upper Beatty sequence for the constant phi^2, where phi = (1 + sqrt(5))/2 and the sequence is generated by floor(n*phi). A054347 = partial sums of the lower Beatty sequence (A000201).
From Michel Dekking, Aug 19 2019: (Start)
Using that floor(k*phi^2) = floor(k*phi)+k, for k=1,...,n, we obtain a(n)/A054347(n) = (A054347(n)+(n*(n+1)/2))/((n*(n+1))))/(A054347(n)/(n*(n+1)/2) [Ambiguous, unbalanced parens - Editors of OEIS] -> (phi/2 + 1/2)/(phi/2), which equals phi.
(End)
FORMULA
a(n) = Sum_{k=1..n} floor(k*phi^2).
a(n) = floor((n*(n+1)/2)*phi^2 - n/2) + (0 or 1). - Benoit Cloitre, Sep 27 2003
EXAMPLE
A001950(1) = 2, then 5, 7, 10, 13, ...; partial sums are 2, 7, 14, 24, 37, ...
MATHEMATICA
a[0] = 0; a[n_] := a[n] = (a[n - 1] + Floor[n*(1 + Sqrt[5])^2/4]); Table[ a[n], {n, 1, 50}] (* Robert G. Wilson v, Sep 27 2003 *)
Accumulate[Floor[GoldenRatio^2 Range[0, 50]]] (* Harvey P. Dale, Aug 11 2021 *)
PROG
(Python)
from math import isqrt
from itertools import islice, count, accumulate
def A088207_gen(): # generator of terms
return accumulate((n+isqrt(5*n**2)>>1)+n for n in count(0))
CROSSREFS
KEYWORD
nonn
AUTHOR
Gary W. Adamson, Sep 23 2003
EXTENSIONS
More terms from Robert G. Wilson v and Benoit Cloitre, Sep 27 2003
STATUS
approved