login
a(n) = floor(n*phi^13), where phi is the golden ratio, A001622.
20

%I #14 Sep 05 2023 18:49:02

%S 0,521,1042,1563,2084,2605,3126,3647,4168,4689,5210,5731,6252,6773,

%T 7294,7815,8336,8857,9378,9899,10420,10941,11462,11983,12504,13025,

%U 13546,14067,14588,15109,15630,16151

%N a(n) = floor(n*phi^13), where phi is the golden ratio, A001622.

%H G. C. Greubel, <a href="/A004928/b004928.txt">Table of n, a(n) for n = 0..10000</a>

%H Tanya Khovanova, <a href="http://www.tanyakhovanova.com/RecursiveSequences/NonRecursions.html">Non Recursions</a>

%H <a href="/index/Be#Beatty">Index entries for sequences related to Beatty sequences</a>

%t Floor[GoldenRatio^(13)*Range[0, 60]] (* _G. C. Greubel_, Sep 05 2023 *)

%o (Magma) [Floor((521+233*Sqrt(5))*n/2): n in [0..60]]; // _G. C. Greubel_, Sep 05 2023

%o (SageMath) [floor(golden_ratio^(13)*n) for n in range(61)] # _G. C. Greubel_, Sep 05 2023

%Y Cf. A004919, A004920, A004921, A004922, A004923, A004924, A004925.

%Y Cf. A004926, A004927, A004929, A004930, A004931, A004932, A004933.

%Y Cf. A004934, A004935, A004976, A066096, A090909.

%Y Cf. A001622.

%K nonn

%O 0,2

%A _N. J. A. Sloane_