A221171 - OEIS

%I #15 Dec 13 2017 08:11:08

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

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

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

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

%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.

%p # sigma01f implemented in A221166

%p A221171 := proc(n)

%p         sigma01f(n,7) ;

%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, 7];

%t Table[a[n], {n, 0, 100}] (* _Jean-François Alcover_, Dec 13 2017, after _R. J. Mathar_ *)

%Y Cf. A221166, A221167, A221168, A221169, A221170.

%K nonn

%O 0,9

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