OFFSET
1,1
COMMENTS
A self-generating sequence.
By calculation it looks as though lim_{m->oo} (1/m)*Sum_{n=1..m} a(n) = sqrt(3).
The 3's are at positions 4, 10, 12, 16, 22, 26, 28, 34, 36, 40, 46, 48, 54, 58, 60, 64, 70, 74, 76, 82, 84, 88, 94, 98,..., which is 2*A059009(n).
From Michel Dekking, Sep 11 2020: (Start)
The conjectured limit above is incorrect, this limit equals 7/4 = 1.75, instead of sqrt(3) = 1.732... This follows from the CLAIM below.
Let nu be the following extension of the 'exchanged' Thue-Morse morphism to a morphism on 3 letters:
nu(1) = 21, nu(2) = 13, nu(3) = 13.
CLAIM: a(567...) = nu(a(234...)).
Here a(567...) denotes the word associated to the sequence a(5),a(6),....
Proof: If a(n)=1 then a(2n+1)=2, and it also follows that a(2n+2)=1. If a(n)>1 then a(2n+2)=3, and it also follows that a(2n+1)=1.
This can also be expressed as nu(a(n))=a(2n+1)a(2n+2). This gives the claim.
(End)
From Michel Dekking, Sep 12 2020: (Start)
The second conjecture above is correct. For a proof, we consider nu^2:
nu^2(1) = 1321, nu^2(2) = 2113, nu^2(3) = 1321.
Let delta be the decoration 0 -> 21, 1 -> 13.
CLAIM: delta(b) = a, where b := A059448.
A proof is involved, as (a(n)) is not a fixed point of nu^2.
First we note that nu(a(n))=a(2n+1)a(2n+2) implies that for all n
nu^2(a(2n-1) )= a(8n-1)a(8n)...a(8n+6).
Let mu be the square of the Thue-Morse morphism:
mu(0) = 0110, mu(1) = 1001.
We use the following property, noted by Shallit in the comments of A059448:
b = 010 mu(b).
First we prove the following lemma on compositions of morphisms.
LEMMA delta mu = nu^2 delta.
Proof: delta(mu(0)) = delta(0110) = 21131321 = nu^2(21) = nu^2(delta(0)), and delta(mu(1)) = delta(1001) = 13212113 = nu^2(13) = nu^2(delta(1)).
To prove the CLAIM, start with
delta(b(123)) = delta(010) = 211321 = a(123456).
We continue with induction, using a recursion formula for b from the Comments of A059448, the LEMMA, and the induction hypothesis:
delta(b(4n)...b(4n+3)) = delta(mu(b(n)) = nu^2(delta(b(n))) =
nu^2(a(2n-1)a(2n)) = a(8n-1)....a(8n+6).
Looking at the positions of 3 in (a(n)), and the positions of 1 in (b(n)), one sees that the claim implies the second conjecture.
(End)
LINKS
Michel Dekking, Table of n, a(n) for n = 1..210
J.-P. Allouche and M. Mendes France, Automata and Automatic Sequences, in: Axel F. and Gratias D. (eds), Beyond Quasicrystals. Centre de Physique des Houches, vol 3. Springer, Berlin, Heidelberg, pp. 293-367, 1995; DOI https://doi.org/10.1007/978-3-662-03130-8_11.
J.-P. Allouche and M. Mendes France, Automata and Automatic Sequences, in: Axel F. and Gratias D. (eds), Beyond Quasicrystals. Centre de Physique des Houches, vol 3. Springer, Berlin, Heidelberg, pp. 293-367, 1995; DOI https://doi.org/10.1007/978-3-662-03130-8_11. [Local copy]
EXAMPLE
a(1) = 2;
for n = 2, a(1) has been assigned: a(1) = 2, thus a(4) = 3;
for n = 3, a(2) is unassigned, thus a(2) = 1 and a(5) = 2;
for n = 4, a(3) is unassigned, thus a(3) = 1 and a(7) = 2; etc.
CROSSREFS
KEYWORD
nonn
AUTHOR
Ctibor O. Zizka, Jan 26 2019
EXTENSIONS
New first element: a(1)=2 instead of a(1)=1. Michel Dekking, Sep 12 2020
STATUS
approved