login
Let phi denote the morphism 0 -> 11, 1 -> 10. This sequence is the limit S(oo) where S(0) = 1; S(n+1) = 1.phi(S(n)).
6

%I #52 Feb 05 2023 09:20:12

%S 1,1,0,1,0,1,1,1,0,1,1,1,0,1,0,1,0,1,1,1,0,1,0,1,0,1,1,1,0,1,1,1,0,1,

%T 1,1,0,1,0,1,0,1,1,1,0,1,1,1,0,1,1,1,0,1,0,1,0,1,1,1,0,1,0,1,0,1,1,1,

%U 0,1,0,1,0,1,1,1,0,1,1,1,0,1,1,1,0,1,0,1,0,1,1,1,0,1,0,1,0,1,1,1,0,1,0,1,0

%N Let phi denote the morphism 0 -> 11, 1 -> 10. This sequence is the limit S(oo) where S(0) = 1; S(n+1) = 1.phi(S(n)).

%C An example of a d-perfect sequence.

%C Concatenation of the bit sequences 1, 10, 1011, 10111010, 1011101010111011, ... used in a construction of A035263 (see Comment there by Benoit Cloitre). - _David Callan_, Oct 08 2005

%C Image, under the coding a,b,d -> 1, c -> 0, of the fixed point, starting with a, of the morphism a -> ab, b -> cd, c -> cd, d -> bb. - _Jeffrey Shallit_, May 15 2016

%H Antti Karttunen, <a href="/A039982/b039982.txt">Table of n, a(n) for n = 0..65537</a>

%H Martin Klazar and Florian Luca, <a href="http://kam.mff.cuni.cz/~klazar/publications.html">On integrality and periodicity of the Motzkin numbers</a>.

%H Martin Klazar and Florian Luca, <a href="http://dx.doi.org/10.1007/s00010-004-2734-x">On integrality and periodicity of the Motzkin numbers</a>, Aequationes Math. 69 (2005), no. 1-2, 68-75.

%H D. Kohel, S. Ling and C. Xing, <a href="http://iml.univ-mrs.fr/~kohel/pub/perfect.pdf">Explicit Sequence Expansions</a>

%H D. Kohel, S. Ling and C. Xing, <a href="http://dx.doi.org/10.1007/978-1-4471-0551-0_23">Explicit Sequence Expansions</a>, Sequences and their Applications, Discrete Mathematics and Theoretical Computer Science 1999, pp 308-317.

%H <a href="/index/Ch#char_fns">Index entries for characteristic functions</a>

%F a(n) = A090344(n) mod 2. - _Christian G. Bower_, Jun 12 2005

%e The first few S(i) are:

%e S(0) = 1

%e S(1) = 1.10 = 110

%e S(2) = 1.101011 = 1101011

%e S(3) = 1.10101110111010 = 110101110111010

%e ...

%t substitutionRule={1->{1, 0}, 0->{1, 1}}; makeSubstitution[seq_]:=Flatten[seq/.substitutionRule]; Flatten[NestList[makeSubstitution, {1}, 5]]

%t NestList[Flatten[ # /. {0 -> {1, 1}, 1 -> {1, 0}}] &, {1}, 6] // Flatten (* _Robert G. Wilson v_, Mar 29 2006 *)

%o (PARI) a(n)=my(A=1+x); for(i=1, n, A=1/(1-x+x*O(x^n))+x^2*A^2+x*O(x^n)); polcoeff(A, n)%2 \\ _Charles R Greathouse IV_, Feb 04 2013

%o (PARI)

%o up_to = 16384;

%o A090344list(up_to) = { my(v=vector(up_to)); v[1] = 1; v[2] = 2; v[3] = 3; for(n=4,up_to,v[n] = ((2*n+2)*v[n-1] -(4*n-6)*v[n-3] +(3*n-4)*v[n-2])/(n+2)); (v); };

%o v090344 = A090344list(up_to);

%o A090344(n) = if(!n,1,v090344[n]);

%o A039982(n) = (A090344(n)%2); \\ _Antti Karttunen_, Sep 27 2018

%o (GAP) b:=[1,1,2];; for n in [4..120] do b[n]:=(1/(n+1))* (2*n*b[n-1]+(3*n-7)*b[n-2]-(4*n-10)*b[n-3]);; od; a:=b mod 2; # _Muniru A Asiru_, Sep 28 2018

%Y Cf. A001006, A035263, A090344.

%K nonn

%O 0,1

%A _N. J. A. Sloane_

%E More terms from _Christian G. Bower_, Jun 12 2005

%E Offset corrected from 1 to 0 by _Antti Karttunen_, Sep 27 2018

%E Entry revised by _N. J. A. Sloane_, Feb 23 2019