login

Year-end appeal: Please make a donation to the OEIS Foundation to support ongoing development and maintenance of the OEIS. We are now in our 61st year, we have over 378,000 sequences, and we’ve reached 11,000 citations (which often say “discovered thanks to the OEIS”).

A004938
a(n) = round(n*phi^3), where phi is the golden ratio, A001622.
3
0, 4, 8, 13, 17, 21, 25, 30, 34, 38, 42, 47, 51, 55, 59, 64, 68, 72, 76, 80, 85, 89, 93, 97, 102, 106, 110, 114, 119, 123, 127, 131, 136, 140, 144, 148, 152, 157, 161, 165, 169, 174, 178, 182, 186, 191, 195, 199
OFFSET
0,2
COMMENTS
Nearest integer to n*phi^3.
LINKS
MATHEMATICA
Round[GoldenRatio^3*Range[0, 80]] (* G. C. Greubel, Oct 17 2023 *)
PROG
(PARI) a(n) = round(n*(((sqrt(5)+1))/2)^3); \\ Michel Marcus, Aug 28 2016
(Magma) [Round((2+Sqrt(5))*n): n in [0..80]]; // G. C. Greubel, Oct 17 2023
(SageMath) [round(golden_ratio^3*n) for n in range(81)] # G. C. Greubel, Oct 17 2023
CROSSREFS
Cf. A001622.
Sequence in context: A311815 A311816 A081843 * A299641 A311817 A311818
KEYWORD
nonn
STATUS
approved