login
Parse the infinite string 012012012012... into distinct phrases 0, 1, 2, 01, 20, 12, 012, ...; a(n) = length of n-th phrase.
2

%I #18 Nov 05 2015 10:53:09

%S 1,1,1,2,2,2,3,4,3,4,3,4,5,5,5,6,7,6,7,6,7,8,8,8,9,10,9,10,9,10,11,11,

%T 11,12,13,12,13,12,13,14,14,14,15,16,15,16,15,16,17,17,17,18,19,18,19,

%U 18,19,20,20,20,21,22,21,22,21,22,23,23,23,24,25,24,25,24,25,26,26,26,27,28,27,28,27,28,29,29,29,30,31,30,31,30,31,32,32,32,33,34,33,34

%N Parse the infinite string 012012012012... into distinct phrases 0, 1, 2, 01, 20, 12, 012, ...; a(n) = length of n-th phrase.

%C See A187180 for details.

%H Ray Chandler, <a href="/A187181/b187181.txt">Table of n, a(n) for n = 1..1000</a>

%H <a href="/index/Rec#order_10">Index entries for linear recurrences with constant coefficients</a>, signature (1, 0, 0, 0, 0, 0, 0, 0, 1, -1).

%F After the initial block of three 1's, the sequence is quasi-periodic with period 9, increasing by 3 after each block.

%F From _Colin Barker_, Nov 05 2015: (Start)

%F a(n) = a(n-1) + a(n-9) - a(n-10) for n>12.

%F G.f.: x*(x^11-x^10-x^8+x^7+x^6+x^3+1) / ((x-1)^2*(x^2+x+1)*(x^6+x^3+1)).

%F (End)

%e The sequence begins

%e 1 1 1

%e 2 2 2 3 4 3 4 3 4

%e 5 5 5 6 7 6 7 6 7

%e 8 8 8 9 10 9 10 9 10

%e 11 11 11 12 13 12 13 12 13

%e 14 14 14 15 16 15 16 15 16 ...

%t Join[{1, 1},LinearRecurrence[{1, 0, 0, 0, 0, 0, 0, 0, 1, -1},{1, 2, 2, 2, 3, 4, 3, 4, 3, 4},98]] (* _Ray Chandler_, Aug 26 2015 *)

%o (PARI) Vec(x*(x^11-x^10-x^8+x^7+x^6+x^3+1)/((x-1)^2*(x^2+x+1)*(x^6+x^3+1)) + O(x^100)) \\ _Colin Barker_, Nov 05 2015

%Y See A187180-A187188 for alphabets of size 2 through 10.

%K nonn,easy

%O 1,4

%A _N. J. A. Sloane_, Mar 06 2011