OFFSET
0,2
COMMENTS
LINKS
G. C. Greubel, Table of n, a(n) for n = 0..10000
FORMULA
a(n) = floor(phi^2*n), where phi = (1+sqrt(5))/2. - Gary Detlefs, Mar 10 2011
MATHEMATICA
(* First program *)
A002251= Fold[Append[#1, #2 Ceiling[#2/GoldenRatio] -Total[#1]] &, {1}, Range[2, 500]] - 1; (* Birkas Gyorgy's code of A019444, modified *)
A090909= Join[{0}, Select[Partition[A002251, 2, 1], #[[2]] > #[[1]] &][[All, 2]]] (* G. C. Greubel, Sep 12 2023 *)
(* Second program *)
Floor[GoldenRatio^2*Range[0, 80]] (* G. C. Greubel, Sep 12 2023 *)
PROG
(Magma) [Floor((3+Sqrt(5))*n/2): n in [0..80]]; // G. C. Greubel, Sep 12 2023
(SageMath) [floor(golden_ratio^2*n) for n in range(81)] # G. C. Greubel, Sep 12 2023
CROSSREFS
KEYWORD
nonn
AUTHOR
Amarnath Murthy, Dec 14 2003
EXTENSIONS
More terms from R. J. Mathar, Sep 29 2017
STATUS
approved