cp's OEIS Frontend

This is a front-end for the Online Encyclopedia of Integer Sequences, made by Christian Perfect. The idea is to provide OEIS entries in non-ancient HTML, and then to think about how they're presented visually. The source code is on GitHub.

A130241 Maximal index k of a Lucas number such that Lucas(k) <= n (the 'lower' Lucas (A000032) Inverse).

This page as a plain text file.
%I A130241 #26 Sep 08 2022 08:45:30
%S A130241 1,1,2,3,3,3,4,4,4,4,5,5,5,5,5,5,5,6,6,6,6,6,6,6,6,6,6,6,7,7,7,7,7,7,
%T A130241 7,7,7,7,7,7,7,7,7,7,7,7,8,8,8,8,8,8,8,8,8,8,8,8,8,8,8,8,8,8,8,8,8,8,
%U A130241 8,8,8,8,8,8,8,9,9,9,9,9,9,9,9,9,9,9,9,9,9,9,9,9,9,9,9,9,9,9,9,9
%N A130241 Maximal index k of a Lucas number such that Lucas(k) <= n (the 'lower' Lucas (A000032) Inverse).
%C A130241 Inverse of the Lucas sequence (A000032), nearly, since a(Lucas(n))=n for n>=1 (see A130242 and A130247 for other versions). For n>=2, a(n)+1 is equal to the partial sum of the Lucas indicator sequence (see A102460). Identical to A130247 except for n=2.
%H A130241 G. C. Greubel, <a href="/A130241/b130241.txt">Table of n, a(n) for n = 1..5000</a>
%F A130241 a(n) = floor(log_phi((n+sqrt(n^2+4))/2)) = floor(arcsinh((n+1)/2)/log(phi)) where phi=(1+sqrt(5))/2.
%F A130241 a(n) = A130242(n+1) - 1 for n>=2.
%F A130241 a(n) = A130247(n) except for n=2.
%F A130241 G.f.: g(x) = 1/(1-x) * Sum{k>=1, x^Lucas(k)}.
%F A130241 a(n) = floor(log_phi(n+1/2)) for n>=2, where phi is the golden ratio.
%e A130241 a(10)=4, since Lucas(4)=7<=10 but Lucas(5)=11>10.
%t A130241 Join[{1}, Table[Floor[Log[GoldenRatio, n + 1/2]], {n, 2, 50}]] (* _G. C. Greubel_, Dec 24 2017 *)
%o A130241 (PARI) for(n=1,50, print1(floor(log((2*n+1)/2)/log((1+sqrt(5))/2)), ", ")) \\ _G. C. Greubel_, Sep 09 2018
%o A130241 (Magma) [Floor(Log((2*n+1)/2)/Log((1+Sqrt(5))/2)): n in [2..50]]; // _G. C. Greubel_, Sep 09 2018
%o A130241 (Python)
%o A130241 from itertools import count, islice
%o A130241 def A130241_gen(): # generator of terms
%o A130241     a, b = 1, 3
%o A130241     for i in count(1):
%o A130241         yield from (i,)*(b-a)
%o A130241         a, b = b, a+b
%o A130241 A130241_list = list(islice(A130241_gen(),40)) # _Chai Wah Wu_, Jun 08 2022
%Y A130241 For partial sums see A130243. Other related sequences: A000032, A130242, A130245, A130247, A130249, A130255, A130259. Indicator sequence A102460. Fibonacci inverse see A130233 - A130240, A104162.
%K A130241 nonn
%O A130241 1,3
%A A130241 _Hieronymus Fischer_, May 19 2007, Jul 02 2007