A183136 -id:A183136 - cp's OEIS frontend

A005614 The binary complement of the infinite Fibonacci word A003849. Start with 1, apply 0->1, 1->10, iterate, take limit.

1, 0, 1, 1, 0, 1, 0, 1, 1, 0, 1, 1, 0, 1, 0, 1, 1, 0, 1, 0, 1, 1, 0, 1, 1, 0, 1, 0, 1, 1, 0, 1, 1, 0, 1, 0, 1, 1, 0, 1, 0, 1, 1, 0, 1, 1, 0, 1, 0, 1, 1, 0, 1, 0, 1, 1, 0, 1, 1, 0, 1, 0, 1, 1, 0, 1, 1, 0, 1, 0, 1, 1, 0, 1, 0, 1, 1, 0, 1, 1, 0, 1, 0, 1, 1, 0, 1, 1, 0, 1, 0, 1, 1, 0, 1, 0, 1, 1, 0

Offset: 0

N. J. A. Sloane

Previous name was: The infinite Fibonacci word (start with 1, apply 0->1, 1->10, iterate, take limit).
Characteristic function of A022342. - Philippe Deléham, May 03 2004
a(n) = number of 0's between successive 1's (see also A003589 and A007538). - Eric Angelini, Jul 06 2005
With offset 1 this is the characteristic sequence for Wythoff A-numbers A000201=[1,3,4,6,...].
Eric Angelini's comment made me think that if 1 is defined to be the number of 0's between successive 1's in a string of 0's and 1's, then this string is 101. Applying the same operation to the digits of 101 leads to 101101, the iteration leads to successive palindromes of lengths given by A001911, up to a(n). - Rémi Schulz, Jul 06 2010
For generalized Fibonacci words see A221150, A221151, A221152, ... - Peter Bala, Nov 11 2013
The limiting mean of the first n terms is phi - 1; the limiting variance is phi (A001622). - Clark Kimberling, Mar 12 2014
Apply the difference operator to every column of the Wythoff difference array, A080164, to get an array of Fibonacci numbers, F(h). Replace each F(h) with h, and apply the difference operator to every column. In the resulting array, every column is A005614. - Clark Kimberling, Mar 02 2015
Binary expansion of the rabbit constant A014565. - M. F. Hasler, Nov 10 2018

			The infinite word is 101101011011010110101101101011...

J.-P. Allouche and J. Shallit, Automatic Sequences, Cambridge Univ. Press, 2003.
G. Melançon, Factorizing infinite words using Maple, MapleTech journal, vol. 4, no. 1, 1997, pp. 34-42, esp. p. 36.

T. D. Noe, Table of n, a(n) for n = 0..10945 (19 iterations)
Ricardo Gómez Aíza, Symbolic dynamical scales: modes, orbitals, and transversals, arXiv:2009.02669 [math.DS], 2020.
F. Axel et al., Vibrational modes in a one dimensional "quasi-alloy": the Morse case, J. de Physique, Colloq. C3, Supp. to No. 7, Vol. 47 (Jul 1986), pp. C3-181-C3-186; see Eq. (10).
Joerg Arndt, Matters Computational (The Fxtbook), 753-754.
E. A. Bender and J. T. Butler, Asymptotic approximations for the number of fanout-free functions, IEEE Trans. Computers, 27 (1978), 1180-1183. (Annotated scanned copy)
M. Bunder and K. Tognetti, On the self matching properties of [j tau], Discrete Math., 241 (2001), 139-151.
J. T. Butler, Letter to N. J. A. Sloane, Dec. 1978.
Glen Joyce C. Dulatre, Jamilah V. Alarcon, Vhenedict M. Florida, and Daisy Ann A. Disu, On Fractal Sequences, DMMMSU-CAS Science Monitor (2016-2017) Vol. 15 No. 2, 109-113.
S. Dulucq and D. Gouyou-Beauchamps, Sur les facteurs des suites de Sturm, Theoret. Comput. Sci. 71 (1990), 381-400.
M. S. El Naschie, Statistical geometry of a Cantor discretum and semiconductors, Computers & Mathematics with Applications, Vol. 29 (Issue 12, June 1995), 103-110.
D. Gault and M. Clint, "Curiouser and curiouser" said Alice. Further reflections on an interesting recursive function, Internat. J. Computer Math., 26 (1988), 35-43.
D. Gault and M. Clint, "Curiouser and curiouser said Alice. Further reflections on an interesting recursive function, Intern. J. Computer. Math., 26 (1988), 35-43. (Annotated scanned copy)
J. Grytczuk, Infinite semi-similar words, Discrete Math. 161 (1996), 133-141.
Clark Kimberling, Intriguing infinite words composed of zeros and ones, Elemente der Mathematik (2021).
Clark Kimberling and K. B. Stolarsky, Slow Beatty sequences, devious convergence, and partitional divergence, Amer. Math. Monthly, 123 (No. 2, 2016), 267-273.
K. L. Kodandapani and S. C. Seth, On combinational networks with restricted fan-out, IEEE Trans. Computers, 27 (1978), 309-318. (Annotated scanned copy)
Ron Knott, The Fibonacci Rabbit Sequence
Wolfdieter Lang, The Tribonacci and ABC Representations of Numbers are Equivalent, arXiv preprint arXiv:1810.09787 [math.NT], 2018.
G. Melançon, Lyndon factorization of sturmian words, Discr. Math., 210 (2000), 137-149.
S. Mneimneh, Fibonacci in The Curriculum: Not Just a Bad Recurrence, in Proceeding SIGCSE '15 Proceedings of the 46th ACM Technical Symposium on Computer Science Education, 253-258.
C. Mongoven, The Rabbit Sequence (a musical composition with A005614).
T. D. Noe, The first 1652 subwords, including leading zeros.
José L. Ramírez, Gustavo N. Rubiano, and Rodrigo de Castro, A Generalization of the Fibonacci Word Fractal and the Fibonacci Snowflake, arXiv preprint arXiv:1212.1368 [cs.DM], 2012-2014.
Jeffrey Shallit, Characteristic words as fixed points of homomorphisms, University of Waterloo Technical Report CS-91-72, 1991.
Jeffrey Shallit, Characteristic words as fixed points of homomorphisms. [Cached copy, with permission]
N. J. A. Sloane, Families of Essentially Identical Sequences, Mar 24 2021 (Includes this sequence)
K. B. Stolarsky, Beatty sequences, continued fractions and certain shift operators, Canad. Math. Bull., 19 (1976), 473-482.
Scott V. Tezlaf, On ordinal dynamics and the multiplicity of transfinite cardinality, arXiv:1806.00331 [math.NT], 2018. See p. 10.
Eric Weisstein's World of Mathematics, Rabbit Constant and Rabbit Sequence.
Index entries for characteristic functions
Index entries for sequences that are fixed points of mappings

Binary complement of A003849, which is the standard form of this sequence.
Two other essentially identical sequences are A096270, A114986.
Subwords: A178992, A171676.
Cf. A000045 (Fibonacci numbers), A001468, A001911, A005206 (partial sums), A014565, A014675, A022342, A036299, A044432, A221150, A221151, A221152.
Cf. A339051 (odd bisection), A339052 (even bisection).
The following sequences are all essentially the same, in the sense that they are simple transformations of each other, with A000201 as the parent: A000201, A001030, A001468, A001950, A003622, A003842, A003849, A004641, A005614, A014675, A022342, A088462, A096270, A114986, A124841. - N. J. A. Sloane, Mar 11 2021

a005614 n = a005614_list !! n
a005614_list = map (1 -) a003849_list
-- Reinhard Zumkeller, Apr 07 2012

[Floor((n+1)*(-1+Sqrt(5))/2)-Floor(n*(-1+Sqrt(5))/2): n in [1..100]]; // Vincenzo Librandi, Jan 17 2019

Digits := 50; u := evalf((1-sqrt(5))/2); A005614 := n->floor((n+1)*u)-floor(n*u);

Nest[ Flatten[ # /. {0 -> {1}, 1 -> {1, 0}}] &, {1}, 10] (* Robert G. Wilson v, Jan 30 2005 *)
Flatten[Nest[{#, #[[1]]} &, {1, 0}, 9]] (* IWABUCHI Yu(u)ki, Oct 23 2013 *)
SubstitutionSystem[{0 -> {1}, 1 -> {1, 0}}, {1, 0}, 9] // Last (* Jean-François Alcover, Feb 06 2020 *)

a(n,w1,s0,s1)=local(w2); for(i=2,n,w2=[ ]; for(k=1,length(w1),w2=concat(w2, if(w1[ k ],s1,s0))); w1=w2); w2
for(n=2,10,print(n" "a(n,[ 0 ],[ 1 ],[ 1,0 ]))) \\ Gives successive convergents to sequence

/* for m>=1 compute exactly A183136(m+1)+1 terms of the sequence */
r=(1+sqrt(5))/2;v=[1,0];for(n=2,m,v=concat(v,vector(floor((n+1)/r),i,v[i]));a(n)=v[n];) /* Benoit Cloitre, Jan 16 2013 */

from math import isqrt
def A005614(n): return (n+isqrt(m:=5*(n+2)**2)>>1)-(n+1+isqrt(m-10*n-15)>>1) # Chai Wah Wu, Aug 17 2022

Define strings S(0)=1, S(1)=10, thereafter S(n)=S(n-1)S(n-2); iterate. Sequence is S(oo). The individual S(n)'s are given in A036299.
a(n) = floor((n+2)*u) - floor((n+1)*u), where u = (-1 + sqrt(5))/2.
Sum_{n>=0} a(n)/2^(n+1) = A014565. - R. J. Mathar, Jul 19 2013
From Peter Bala, Nov 11 2013: (Start)
If we read the present sequence as the digits of a decimal constant c = 0.101101011011010 ... then we have the series representation c = Sum_{n >= 1} 1/10^floor(n*phi). An alternative representation is c = Sum_{n >= 1} 1/10^floor(n/phi) - 10/9.
The constant 9*c has the simple continued fraction representation [0; 1, 10, 10, 100, 1000, ..., 10^Fibonacci(n), ...]. See A010100.
Using this result we can find the alternating series representation c = 1/9 - 9*Sum_{n >= 1} (-1)^(n+1)*(1 + 10^Fibonacci(3*n+1))/( (10^(Fibonacci(3*n - 1)) - 1)*(10^(Fibonacci(3*n + 2)) - 1) ). The series converges very rapidly: for example, the first 10 terms of the series give a value for c accurate to more than 5.7 million decimal places. Cf. A014565. (End)
a(n) = A005206(n+1) - A005206(n). a(2*n) = A339052(n); a(2*n+1) = A339051(n+1). - Peter Bala, Aug 09 2022

Corrected by Clark Kimberling, Oct 04 2000

Name corrected by Michel Dekking, Apr 02 2019

A060143 a(n) = floor(n/tau), where tau = (1 + sqrt(5))/2.

Original entry on oeis.org

0, 0, 1, 1, 2, 3, 3, 4, 4, 5, 6, 6, 7, 8, 8, 9, 9, 10, 11, 11, 12, 12, 13, 14, 14, 15, 16, 16, 17, 17, 18, 19, 19, 20, 21, 21, 22, 22, 23, 24, 24, 25, 25, 26, 27, 27, 28, 29, 29, 30, 30, 31, 32, 32, 33, 33, 34, 35, 35, 36, 37, 37, 38, 38, 39, 40, 40, 41, 42, 42, 43, 43, 44, 45, 45

Offset: 0

Clark Kimberling, Mar 05 2001

Fibonacci base shift right: for n >= 0, a(n+1) = Sum_{k in A_n} F_{k-1}, where n = Sum_{k in A_n} F_k (unique) expression of n as a sum of "noncontiguous" Fibonacci numbers (with index >=2). - Michele Dondi (bik.mido(AT)tiscalenet.it), Dec 30 2001 [corrected, and aligned with sequence offset by Peter Munn, Jan 10 2018]
Numerators a(n) of fractions slowly converging to phi, the golden ratio: let a(1) = 0, b(n) = n - a(n); if (a(n) + 1) / b(n) < (1 + sqrt(5))/2, then a(n+1) = a(n) + 1, else a(n+1) = a(n). a(n) + b(n) = n and as n -> +infinity, a(n) / b(n) converges to (1 + sqrt(5))/2. For all n, a(n) / b(n) < (1 + sqrt(5))/2. - Robert A. Stump (bee_ess107(AT)msn.com), Sep 22 2002
a(10^n) gives the first few digits of phi=(sqrt(5)-1)/2.
Comment corrected, two alternative ways, by Peter Munn, Jan 10 2018: (Start)
(a(n) = a(n+1) or a(n) = a(n-1)) if and only if a(n) is in A066096.
a(n+1) = a(n+2) if and only if n is in A003622.
(End)
From Wolfdieter Lang, Jun 28 2011: (Start)
a(n+1) counts for n >= 1 the number of Wythoff A-numbers not exceeding n.
a(n+1) counts also the number of Wythoff B-numbers smaller than A(n+2), with the Wythoff A- and B-sequences A000201 and A001950, respectively.
a(n+1) = Sum_{j=1..n} A005614(j-1) for n >= 1 (no rounding problems like in the above definition, because the rabbit sequence A005614(n-1) for n >= 1, can be defined by a substitution rule).
a(n+1) = A(n+1)-(n+1) (serving, together with the last equation, as definition for A(n+1), given the rabbit sequence).
a(n+1) = A005206(n), n >= 0.
(End)
Let b(n) = floor((n+1)/phi). Then b(n) + b(b(n-1)) = n [Granville and Rasson]. - N. J. A. Sloane, Jun 13 2014

			a(6)= 3 so b(6) = 6 - 3 = 3. a(7) = 4 because (a(6) + 1) / b(6) = 4/3 which is < (1 + sqrt(5))/2. So b(7) = 7 - 4 = 3. a(8) = 4 because (a(7) + 1) / b(7) = 5/3 which is > (1 + sqrt(5))/2. - Robert A. Stump (bee_ess107(AT)msn.com), Sep 22 2002
From _Wolfdieter Lang_, Jun 28 2011: (Start)
There are a(4) = 2 (positive) Wythoff A-numbers <= 3, namely 1 and 3.
There are a(4) = 2 (positive) Wythoff B-numbers < A(4) = 6, namely 2 and 5.
a(4) = 2 = A(4) - 4 = 6 - 4.
(End)

William A. Tedeschi, Table of n, a(n) for n = 0..10000 [This replaces an earlier b-file computed by Harry J. Smith]
M. Celaya and F. Ruskey (Proposers), Another property of only the golden ratio, Problem 11651, Amer. Math. Monthly, 121 (2014), 550-551.
F. Michel Dekking, Morphisms, Symbolic Sequences, and Their Standard Forms, Journal of Integer Sequences, Vol. 19 (2016), Article 16.1.1.
Vincent Granville and Jean-Paul Rasson, A strange recursive relation, J. Number Theory 30 (1988), no. 2, 238--241. MR0961919(89j:11014). - From _N. J. A. Sloane_, Jun 13 2014
Jeffrey Shallit, The Hurt-Sada Array and Zeckendorf Representations, arXiv:2501.08823 [math.NT], 2025. See p. 2.

Cf. A000045 (Fibonacci numbers), A003622, A022342, A035336.
Terms that occur only once: A001950.
Terms that occur twice: A066096 (a version of A000201).
Numerator sequences for other values, as described in Robert A. Stump's 2002 comment: A074065 (sqrt(3)), A074840 (sqrt(2)).
Apart from initial terms, same as A005206.
First differences: A096270 (a version of A005614).
Partial sums: A183136.

[Floor(2*n/(1+Sqrt(5))): n in [0..80]]; // Vincenzo Librandi, Mar 29 2015

Floor[Range[0,80]/GoldenRatio] (* Harvey P. Dale, May 09 2013 *)

{ default(realprecision, 10); p=(sqrt(5) - 1)/2; for (n=0, 1000, write("b060143.txt", n, " ", floor(n*p)); ) } \\ Harry J. Smith, Jul 02 2009

from math import isqrt
def A060143(n): return (n+isqrt(5*n**2)>>1)-n # Chai Wah Wu, Aug 10 2022

a(n) = floor(phi(n)), where phi=(sqrt(5)-1)/2. [corrected by Casey Mongoven, Jul 18 2008]
a(F_n + 1) = F_{n-1} if F_n is the n-th Fibonacci number. [aligned with sequence offset by Peter Munn, Jan 10 2018]
a(1) = 0. b(n) = n - a(n). If (a(n) + 1) / b(n) < (1 + sqrt(5))/2, then a(n+1) = a(n) + 1, else a(n+1) = a(n). - Robert A. Stump (bee_ess107(AT)msn.com), Sep 22 2002 [corrected by Peter Munn, Jan 07 2018]
A006336(n) = A006336(n-1) + A006336(a(n)) for n>1. - Reinhard Zumkeller, Oct 24 2007
a(n) = floor(n*phi) - n, where phi = (1+sqrt(5))/2. - William A. Tedeschi, Mar 06 2008
Celaya and Ruskey give an interesting formula for a(n). - N. J. A. Sloane, Jun 13 2014

I merged three identical sequences to create this entry. Some of the formulas may need their initial terms adjusting now. - N. J. A. Sloane, Mar 05 2003

More terms from William A. Tedeschi, Mar 06 2008

f(m) is an approximation to A183136(m) = Sum_{k=1..m} floor(k/phi) based on assuming the floor in each term decreases it by 1/2 from what is otherwise a triangular sum; and further offset + 1/(2*phi) in f(m) chosen to improve the accuracy of this approximation.

The actual values of frac(k/phi) can differ from 1/2 each by a net amount which is enough to make m a term of this sequence.

cp's OEIS Frontend

A384947 Positive integers m for which A183136(m) != f(m), where f(m) = floor( (m*(m+1)/2)/phi - m/2 + 1/(2*phi) ) and phi = (1+sqrt(5))/2 is the golden ratio.

Views

Author

Keywords

Comments

Examples

Links

Crossrefs

Programs

Mathematica

A005614 The binary complement of the infinite Fibonacci word A003849. Start with 1, apply 0->1, 1->10, iterate, take limit.

Views

Author

Keywords

Comments

Examples

References

Links

Crossrefs

Programs

Haskell

Magma

Maple

Mathematica

PARI

PARI

Python

Formula

Extensions

A060143 a(n) = floor(n/tau), where tau = (1 + sqrt(5))/2.

Views

Author

Keywords

Comments

Examples

Links

Crossrefs

Programs

Magma

Mathematica

PARI

Python

Formula

Extensions

A183143 [1/r]+[2/r]+...+[n/r], where r=sqrt(3) and []=floor.

Views

Author

Keywords

Comments

Crossrefs

Programs

Maple

Mathematica

PARI

Formula

A183137 a(n) = [1/s] + [2/s] + ... + [n/s], where s = (golden ratio)^2 = (3+sqrt(5))/2 and [] = floor.

Views

Author

Keywords

Comments

Examples

Links

Crossrefs

Programs

Mathematica

Formula

A385177 a(n) = Sum_{k=1..n} ceiling(k/phi), where phi is the golden ratio (A001622).

Views

Author

Keywords

Links

Crossrefs

Programs

Mathematica

Python

Formula

A384947 Positive integers m for which A183136(m) != f(m), where f(m) = floor( (m(m+1)/2)/phi - m/2 + 1/(2phi) ) and phi = (1+sqrt(5))/2 is the golden ratio.