login
a(n) = ceiling(n*exp((1+sqrt(5))/2)).
2

%I #13 Sep 08 2022 08:45:34

%S 6,11,16,21,26,31,36,41,46,51,56,61,66,71,76,81,86,91,96,101,106,111,

%T 116,122,127,132,137,142,147,152,157,162,167,172,177,182,187,192,197,

%U 202,207,212,217,222,227,232,238,243,248,253,258,263,268,273,278,283

%N a(n) = ceiling(n*exp((1+sqrt(5))/2)).

%H G. C. Greubel, <a href="/A140232/b140232.txt">Table of n, a(n) for n = 1..10000</a>

%F a(n) = ceiling(n*A139341). - _R. J. Mathar_, Feb 06 2009

%t Ceiling[Exp[GoldenRatio]*Range[60]] (* _G. C. Greubel_, Jun 30 2019 *)

%o (PARI) phi=(1+sqrt(5))/2; vector(60, n, ceil(n*exp(phi)) ) \\ _G. C. Greubel_, Jun 30 2019

%o (Magma) phi:=(1+Sqrt(5))/2; [Ceiling(n*Exp(phi)): n in [1..60]]; // _G. C. Greubel_, Jun 30 2019

%o (Sage) [ceil(n*exp(golden_ratio)) for n in (1..60)] # _G. C. Greubel_, Jun 30 2019

%Y Cf. A001622, A094214, A104457, A098317, A002390, A139339, A139340, A139341, A139342, A139345, A139346, A139347, A139348, A139349, A139350, A140231, A016861.

%K nonn

%O 1,1

%A _Mohammad K. Azarian_, May 13 2008