%I #99 Jun 04 2024 15:36:18
%S 1,2,3,3,4,4,4,5,5,5,5,5,6,6,6,6,6,6,6,6,7,7,7,7,7,7,7,7,7,7,7,7,7,8,
%T 8,8,8,8,8,8,8,8,8,8,8,8,8,8,8,8,8,8,8,8,9,9,9,9,9,9,9,9,9,9,9,9,9,9,
%U 9,9,9,9,9,9,9,9,9,9,9,9,9,9,9,9,9,9,9,9,10,10,10,10,10,10,10,10,10,10,10
%N n occurs Fibonacci(n) times (cf. A000045).
%C Number of digits in Zeckendorf-binary representation of n. E.g., the Zeckendorf representation of 12 is 8+3+1, which in binary notation is 10101, which consists of 5 digits. - _Clark Kimberling_, Jun 05 2004
%C First position where value n occurs is A000045(n+1), i.e., a(A000045(n)) = n-1, for n >= 2 and a(A000045(n)-1) = n-2, for n >= 3.
%C This is the number of distinct Fibonacci numbers greater than 0 which are less than or equal to n. - _Robert G. Wilson v_, Dec 10 2006
%C The smallest nondecreasing sequence a(n) such that a(Fibonacci(n-1)) = n. - _Tanya Khovanova_, Jun 20 2007
%H Reinhard Zumkeller, <a href="/A072649/b072649.txt">Table of n, a(n) for n = 1..10000</a>
%H A. J. Macfarlane, <a href="https://arxiv.org/abs/2405.18128">On the fibbinary numbers and the Wythoff array</a>, arXiv:2405.18128 [math.CO], 2024. See page 4.
%H Popular Computing (Calabasas, CA), <a href="/A256967/a256967.png">A Coding Exercise (from a suggestion by R. W. Hamming)</a>, Vol. 5 (No. 54, Sep 1977), p. PC55-18.
%F G.f.: (Sum_{n>1} x^Fibonacci(n))/(1-x). - _Michael Somos_, Apr 25 2003
%F From _Hieronymus Fischer_, May 02 2007: (Start)
%F a(n) = floor(log_phi((sqrt(5)*n + sqrt(5*n^2+4))/2)) - 1, where phi is A001622.
%F a(n) = floor(arcsinh(sqrt(5)*n/2)/log(phi)) - 1.
%F a(n) = A108852(n) - 2. (End)
%F a(n) = -1 + floor( log_phi( (n+0.2)*sqrt(5) ) ), where log_phi(x) is the logarithm to the base (1+sqrt(5))/2. - _Ralf Stephan_, May 14 2007
%F Sum_{n>=1} (-1)^(n+1)/a(n) = Pi/(3*sqrt(3)) (A073010). - _Amiram Eldar_, Feb 18 2024
%e 1, 1, then F(2) 2's, then F(3) 3's, then F(4) 4's, ..., then F(k) k's, ...
%p A072649 := proc(n)
%p local j;
%p for j from ilog[(1+sqrt(5))/2](n) while combinat[fibonacci](j+1)<=n do
%p end do;
%p j-1
%p end proc:
%p seq(A072649(n), n=1..120); # _Alois P. Heinz_, Mar 18 2013
%t Table[Table[n, {Fibonacci[n]}], {n, 10}] // Flatten (* _Robert G. Wilson v_, Jan 14 2007 *)
%t a[n_] := Module[{j}, For[j = Floor@Log[GoldenRatio, n], Fibonacci[j+1] <= n, j++]; j-1];
%t Table[a[n], {n, 1, 120}] (* _Jean-François Alcover_, Nov 17 2022, after _Alois P. Heinz_ *)
%o (PARI) a(n) = -1+floor(log(((n+0.2)*sqrt(5)))/log((1+sqrt(5))/2))
%o (PARI) a(n)=local(m); if(n<1,0,m=0; until(fibonacci(m)>n,m++); m-2)
%o (Haskell)
%o a072649 n = a072649_list !! (n-1)
%o a072649_list = f 1 where
%o f n = (replicate (fromInteger $ a000045 n) n) ++ f (n+1)
%o -- _Reinhard Zumkeller_, Jul 04 2011
%o (Python)
%o from sympy import fibonacci
%o def a(n):
%o if n<1: return 0
%o m=0
%o while fibonacci(m)<=n: m+=1
%o return m-2
%o print([a(n) for n in range(1, 101)]) # _Indranil Ghosh_, Jun 09 2017
%o (MIT/GNU Scheme) (define (A072649 n) (let ((b (A072648 n))) (+ -1 b (floor->exact (/ n (A000045 (1+ b))))))) ;; (The implementation below is better)
%o (Scheme) (define (A072649 n) (if (<= n 3) n (let loop ((k 5)) (if (> (A000045 k) n) (- k 2) (loop (+ 1 k)))))) ;; (Use this with the memoized implementation of A000045 given under that entry. No floating point arithmetic is involved). - _Antti Karttunen_, Oct 06 2017
%Y Cf. A000045, A095791, A095792.
%Y Cf. A001622 (golden ratio phi), A073010.
%Y Used to construct A003714. Cf. also A002024, A072643, A072648, A072650.
%Y Cf. A131234.
%Y Partial sums: A256966, A256967.
%K nonn,look
%O 1,2
%A _Antti Karttunen_, Jun 02 2002
%E Typo fixed by _Charles R Greathouse IV_, Oct 28 2009