|
|
A130260
|
|
Minimal index k of an even Fibonacci number A001906 such that A001906(k) = Fib(2k) >= n (the 'upper' even Fibonacci Inverse).
|
|
8
|
|
|
0, 1, 2, 2, 3, 3, 3, 3, 3, 4, 4, 4, 4, 4, 4, 4, 4, 4, 4, 4, 4, 4, 5, 5, 5, 5, 5, 5, 5, 5, 5, 5, 5, 5, 5, 5, 5, 5, 5, 5, 5, 5, 5, 5, 5, 5, 5, 5, 5, 5, 5, 5, 5, 5, 5, 5, 6, 6, 6, 6, 6, 6, 6, 6, 6, 6, 6, 6, 6, 6, 6, 6, 6, 6, 6, 6, 6, 6, 6, 6, 6, 6, 6, 6, 6, 6, 6, 6, 6, 6, 6, 6, 6, 6, 6, 6, 6, 6, 6, 6, 6
(list;
graph;
refs;
listen;
history;
text;
internal format)
|
|
|
OFFSET
|
0,3
|
|
COMMENTS
|
Inverse of the even Fibonacci sequence (A001906), since a(A001906(n))=n (see A130259 for another version).
a(n+1) is the number of even Fibonacci numbers (A001906) <=n.
|
|
LINKS
|
|
|
FORMULA
|
a(n) = ceiling(arcsinh(sqrt(5)*n/2)/(2*log(phi))) for n>=0.
a(n) = ceiling(arccosh(sqrt(5)*n/2)/(2*log(phi))) for n>=1.
a(n) = ceiling(log_phi(sqrt(5)*n)/2)=ceiling(log_phi(sqrt(5)*n-1)/2) for n>=1, where phi=(1+sqrt(5))/2.
G.f.: g(x)=x/(1-x)*Sum_{k>=0} x^Fib(2*k).
|
|
EXAMPLE
|
|
|
MATHEMATICA
|
Join[{0}, Table[Ceiling[Log[GoldenRatio, Sqrt[5]*n]/2], {n, 1, 100}]] (* G. C. Greubel, Sep 12 2018 *)
|
|
PROG
|
(PARI) for(n=0, 100, print1(if(n==0, 0, ceil(log(sqrt(5)*n)/(2*log((1+ sqrt(5))/2)))), ", ")) \\ G. C. Greubel, Sep 12 2018
(Magma) [0] cat [Ceiling(Log(Sqrt(5)*n)/(2*Log((1+ Sqrt(5))/2))): n in [1..100]]; // G. C. Greubel, Sep 12 2018
|
|
CROSSREFS
|
|
|
KEYWORD
|
nonn
|
|
AUTHOR
|
|
|
STATUS
|
approved
|
|
|
|