(Translated by https://www.hiragana.jp/)
# Greetings from The On-Line Encyclopedia of Integer Sequences! http://oeis.org/
Search: id:a026465
Showing 1-1 of 1
%I A026465 #121 Aug 26 2023 08:56:40
%S A026465 1,2,1,1,2,2,2,1,1,2,1,1,2,1,1,2,2,2,1,1,2,2,2,1,1,2,2,2,1,1,2,1,1,2,
%T A026465 1,1,2,2,2,1,1,2,1,1,2,1,1,2,2,2,1,1,2,1,1,2,1,1,2,2,2,1,1,2,2,2,1,1,
%U A026465 2,2,2,1,1,2,1,1,2,1,1,2,2,2,1,1,2,2,2,1,1,2,2,2,1,1,2,1,1,2,1,1,2,2,2,1,1
%N A026465 Length of n-th run of identical symbols in the Thue-Morse sequence A010060 (or A001285).
%C A026465 It appears that the sequence can be calculated by any of the following methods:
%C A026465 (1) Start with 1 and repeatedly replace 1 with 1, 2, 1 and 2 with 1, 2, 2, 2, 1;
%C A026465 (2) a(1) = 1, all terms are either 1 or 2 and, for n > 0, a(n) = 1 if the length of the n-th run of 2's is 1; a(n) = 2 if the length of the n-th run of consecutive 2's is 3, with each run of 2's separated by a run of two 1's;
%C A026465 (3) replace each 3 in A080426 with 2. - _John W. Layman_, Feb 18 2003
%C A026465 Number of representations of n as a sum of Jacobsthal numbers (1 is allowed twice as a part). Partial sums are A003159. With interpolated zeros, g.f. is (Product_{k>=1} (1 + x^A078008(k)))/2. - _Paul Barry_, Dec 09 2004
%C A026465 In other words, the consecutive 0's or 1's in A010060 or A010059. - Robin D. Saunders (saunders_robin_d(AT)hotmail.com), Sep 06 2006
%C A026465 From _Carlo Carminati_, Feb 25 2011: (Start)
%C A026465 The sequence (starting with the second term) can also be calculated by the following method:
%C A026465 Apply repeatedly to the string S_0 = [2] the following algorithm: take a string S, double it, if the last figure is 1, just add the last figure to the previous one, if the last figure is greater than one, decrease it by one unit and concatenate a figure 1 at the end. (This algorithm is connected with the interpretation of the sequence as a continued fraction expansion.) (End)
%C A026465 This sequence, starting with the second term, happens to be the continued fraction expansion of the biggest cluster point of the set {x in [0,1]: F^k(x) >= x, for all k in N}, where F denotes the Farey map (see A187061). - _Carlo Carminati_, Feb 28 2011
%C A026465 Starting with the second term, the fixed point of the substitution 2 -> 211, 1 -> 2. - _Carlo Carminati_, Mar 03 2011
%C A026465 It appears that this sequence contains infinitely many distinct palindromic subsequences. - _Alexander R. Povolotsky_, Oct 30 2016
%C A026465 From _Michel Dekking_, Feb 13 2019: (Start)
%C A026465 Let tau defined by tau(0) = 01, tau(1) = 10 be the Thue-Morse morphism, with fixed point A010060. Consecutive runs in A010060 are 0, 11, 0, 1, 00, 1, 1, ..., which are coded by their lengths 1, 2, 1, 1, 2, ... Under tau^2 consecutive runs are mapped to consecutive runs:
%C A026465 tau^2(0) = 0110, tau^2(1) = 1001,
%C A026465 tau^2(00) = 01100110, tau^2(11) = 10011001.
%C A026465 The reason is that (by definition of a run!) runs of 0's and runs of 1's alternate in the sequence of runs, and this is inherited by the image of these runs under tau^2.
%C A026465 Under tau^2 the runs of length 1 are mapped to the sequence 1,2,1 of run lengths, and the runs of length 2 are mapped to the sequence 1,2,2,2,1 of run lengths. This proves John Layman's conjecture number (1): it follows that (a(n)) is fixed point of the morphism alpha
%C A026465 alpha: 1 -> 121, 2 -> 12221.
%C A026465 Since alpha(1) and alpha(2) are both palindromes, this also proves Alexander Povolotsky's conjecture.
%C A026465 (End)
%H A026465 N. J. A. Sloane, Table of n, a(n) for n = 1..10000
%H A026465 J.-P. Allouche, Andre Arnold, Jean Berstel, Srecko Brlek, William Jockusch, Simon Plouffe and Bruce E. Sagan, A sequence related to that of Thue-Morse, Discrete Math., Vol. 139, No. 1-3 (1995), pp. 455-461.
%H A026465 G. Allouche, J.-P. Allouche and J. Shallit, Kolam indiens, dessins sur le sable aux îles Vanuatu, courbe de Sierpinski et morphismes de monoïde, Ann. Inst. Fourier (Grenoble), Vol. 56, No. 7 (2006), pp. 2115-2130.
%H A026465 Claudio Bonanno, Carlo Carminati, Stefano Isola and Giulio Tiozzo, Dynamics of continued fractions and kneading sequences of unimodal maps, arXiv:1012.2131 [math.DS], 2010-2012.
%H A026465 Srećko Brlek, Enumeration of factors in the Thue-Morse word, Discrete Applied Math., Vol. 24, No. 1-3 (1989), pp. 83-96.
%H A026465 Julien Cassaigne, Limit values of the recurrence quotient of Sturmian sequences, Theoret. Comput. Sci., Vol. 218, No. 1 (1999), pp. 3-12.
%H A026465 Artūras Dubickas, On the distance from a rational power to the nearest integer, Journal of Number Theory, Vol. 117, No. 1 (March 2006), pp. 222-239.
%H A026465 Artūras Dubickas, On a sequence related to that of Thue-Morse and its applications, Discrete Math., Vol. 307, No. 9-10 (2007), pp. 1082-1093. MR2292537 (2008b:11086).
%H A026465 Cor Kraaikamp, Thomas A. Schmidt and Wolfgang Steiner, Natural extensions and entropy of alpha-continued fractions, arXiv:1011.4283 [math.DS], 2010-2012.
%H A026465 Claude Lenormand, Deux transformations sur les mots, Preprint, 5 pages, Nov 17 2003. Apparently unpublished. This is a scanned copy of the version that the author sent to me in 2003. - _N. J. A. Sloane_, Sep 09 2018. See page 2.
%H A026465 Kevin Ryde, PARI/GP Code
%H A026465 Jeffrey Shallit, Automaticity IV: Sequences, sets, and diversity, J. Théor. Nombres Bordeaux, Vol. 8, No. 2 (1996), pp. 347-367. See page 354.
%F A026465 a(1) = 1; for n > 1, a(n) = A003159(n) - A003159(n-1). - _Benoit Cloitre_, May 31 2003
%F A026465 G.f.: Product_{k>=1} (1 + x^A001045(k)). - _Paul Barry_, Dec 09 2004
%F A026465 Asymptotic mean: lim_{m->oo} (1/m) * Sum_{k=1..m} a(k) = 3/2. - _Amiram Eldar_, Jan 16 2022
%p A026465 # From _Carlo Carminati_, Feb 25 2011:
%p A026465 ## period-doubling routine:
%p A026465 double:=proc(SS)
%p A026465 NEW:=[op(S), op(S)]:
%p A026465 if op(nops(NEW),NEW)=1
%p A026465 then NEW:=[seq(op(j,NEW), j=1..nops(NEW)-2),op(nops(NEW)-1,NEW)+1]:
%p A026465 else NEW:=[seq(op(j,NEW), j=1..nops(NEW)-1),op(nops(NEW)-1,NEW)-1,1]:
%p A026465 fi:
%p A026465 end proc:
%p A026465 # 10 loops of the above routine generate the first 1365 terms of the sequence
%p A026465 # (except for the initial term):
%p A026465 S:=[2]:
%p A026465 for j from 1 to 10 do S:=double(S); od:
%p A026465 S;
%p A026465 # From _N. J. A. Sloane_, Dec 31 2013:
%p A026465 S:=[b]; M:=14;
%p A026465 for n from 1 to M do T:=subs({b=[b,a,a], a=[b]}, S);
%p A026465 S := map(x->op(x),T); od:
%p A026465 T:=subs({a=1,b=2},S): T:=[1,op(T)]: [seq(T[n],n=1..40)];
%t A026465 Length /@ Split@ Nest[ Flatten@ Join[#, # /. {1 -> 2, 2 -> 1}] &, {1}, 7]
%t A026465 NestList[ Flatten[# /. {1 -> {2}, 2 -> {1, 1, 2}}] &, {1}, 7] // Flatten (* _Robert G. Wilson v_, May 20 2014 *)
%o A026465 (Haskell)
%o A026465 import Data.List (group)
%o A026465 a026465 n = a026465_list !! (n-1)
%o A026465 a026465_list = map length $ group a010060_list
%o A026465 -- _Reinhard Zumkeller_, Jul 15 2014
%o A026465 (PARI) See links.
%Y A026465 Cf. A010060, A001285, A101615, A026490 (run lengths).
%Y A026465 A080426 is an essentially identical sequence with another set of constructions.
%Y A026465 Cf. A104248 (bisection odious), A143331 (bisection evil), A003159 (partial sums).
%Y A026465 Cf. A187061, A363361 (as continued fraction).
%K A026465 nonn,eigen
%O A026465 1,2
%A A026465 _Clark Kimberling_
%E A026465 Corrected and extended by _John W. Layman_, Feb 18 2003
%E A026465 Definition revised by _N. J. A. Sloane_, Dec 30 2013
# Content is available under The OEIS End-User License Agreement: http://oeis.org/LICENSE