A221166 - OEIS

%I #21 Nov 21 2017 03:08:52

%S 0,1,0,3,0,3,2,3,0,3,0,1,0,1,2,1,0,1,0,3,0,1,0,1,2,1,2,3,2,1,2,1,0,1,

%T 0,3,0,1,0,1,2,1,0,1,0,3,0,3,2,3,0,3,0,1,0,3,0,3,2,3,2,1,2,3,2,3,0,3,

%U 0,1,0,3,0,3,2,3,0,3,0,1,0,1,2,1,0,1,0,3,0,3

%N The infinite generalized Fibonacci word p^[2].

%H José L. Ramírez, Gustavo N. Rubiano, and Rodrigo de Castro, <a href="http://arxiv.org/abs/1212.1368">A Generalization of the Fibonacci Word Fractal and the Fibonacci Snowflake</a>, arXiv preprint arXiv:1212.1368 [cs.DM], 2012-2014.

%e The infinite Fibonacci word f^[2] is A003849. If we apply the morphism {1,0}->{0,2} we have 2, 0, 2, 2, 0 ,2 ... Prepending a 1 and replacing the sequence with the partial sums plus 1 (mod 4), applying operator sigma_1, we have 1, 3, 3, 1, 3, 3, 1, 1, 3, 1. Finally prepending 0 and replacing the that sequence with the partial sums (mod 4), applying operator sigma_0, we have the a(n). - _R. J. Mathar_, Jul 09 2013

%p # fibi and fibonni defined in A221150

%p fmorph := proc(n,i)

%p     if fibonni(n,i) = 0 then

%p         2;

%p     else

%p         0 ;

%p     end if;

%p end proc:

%p sigma1f := proc(n,i)

%p     if n = 0 then

%p         1;

%p     else

%p         1 + modp(add(fmorph(j,i),j=0..n-1),4) ;

%p     end if;

%p end proc:

%p sigma01f := proc(n,i)

%p     if n = 0 then

%p         0;

%p     else

%p         modp(add(sigma1f(j,i),j=0..n-1),4) ;

%p     end if;

%p end proc:

%p A221166 := proc(n)

%p     sigma01f(n,2) ;

%p end proc: # _R. J. Mathar_, Jul 09 2013

%t fibi[n_, i_] := fibi[n, i] = Which[n == 0, {0}, n == 1, Append[Table[0, {j, 1, i - 1}], 1], True, Join[fibi[n - 1, i], fibi[n - 2, i]]];

%t fibonni[n_, i_] := fibonni[n, i] = Module[{fn, Fn}, For[fn = 0, True, fn++, Fn = fibi[fn, i]; If[Length[Fn] >= n + 1 && Length[Fn] > i + 3, Return[ Fn[[n + 1]]]]]];

%t fmorph[n_, i_] := If[fibonni[n, i] == 0, 2, 0];

%t sigma1f[n_, i_] := If[n == 0, 1, 1+Mod[Sum[fmorph[j, i], {j, 0, n-1}], 4]];

%t sigma01f[n_, i_] := If[n == 0, 0, Mod[Sum[sigma1f[j, i], {j, 0, n-1}], 4]];

%t a[n_] := sigma01f[n, 2]; Table[a[n], {n, 0, 100}] (* _Jean-François Alcover_, Nov 21 2017, after _R. J. Mathar_ *)

%Y Cf. A221166-A221171.

%K nonn

%O 0,4

%A _N. J. A. Sloane_, Jan 04 2013

%E Changed name from p^[1] to p^[2] because p^[1] could not be reproduced. - _R. J. Mathar_, Jul 09 2013