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

%I #38 Aug 25 2023 08:50:17

%S 0,29,58,87,116,145,174,203,232,261,290,319,348,377,406,435,464,493,

%T 522,551,580,609,638,667,696,725,754,783,812,841,871,900,929,958,987,

%U 1016,1045,1074,1103,1132,1161,1190,1219,1248,1277,1306,1335,1364

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

%H Karl V. Keller, Jr., <a href="/A004922/b004922.txt">Table of n, a(n) for n = 0..1000</a>

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

%t Table[Floor[n ((1 + Sqrt[5])/2)^7], {n, 0, 50}] (* _Vincenzo Librandi_, Jul 22 2015 *)

%o (Python)

%o from sympy import sqrt

%o phi = (1 + sqrt(5))/2

%o for n in range(0,101): print(int(n*phi**7), end=', ') # _Karl V. Keller, Jr._, Jul 22 2015

%o (Magma) [Floor(n*((1 + Sqrt(5))/2)^7): n in [0..50]]; // _Vincenzo Librandi_, Jul 22 2015

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

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

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

%Y Cf. A001622, A004962 (ceiling).

%K nonn

%O 0,2

%A _N. J. A. Sloane_