OFFSET
1,3
COMMENTS
Taking every sixteenth term gives the same sequence: (and similarly for all higher powers of 16 as well): a(n) = a(16*n).
LINKS
A. Karttunen, Table of n, a(n) for n = 1..65536
MATHEMATICA
HC = {L[n_ /; IntegerQ[n/2]] :> {F[n], L[n], L[n + 1], R[n + 2]},
R[n_ /; IntegerQ[(n + 1)/2]] :> {F[n], R[n], R[n + 3], L[n + 2]},
R[n_ /; IntegerQ[n/2]] :> {L[n], R[n + 1], R[n], F[n + 3]},
L[n_ /; IntegerQ[(n + 1)/2]] :> {R[n], L[n + 3], L[n], F[n + 1]},
F[n_ /; IntegerQ[n/2]] :> {L[n], R[n + 1], R[n], L[n + 3]},
F[n_ /; IntegerQ[(n + 1)/2]] :> {R[n], L[n + 3], L[n], R[n + 1]}};
a[1] = F[0]; Map[(a[n_ /; IntegerQ[(n - #)/16]] :=
Part[Flatten[a[(n + 16 - #)/16] /. HC /. HC], #]) &, Range[16]];
Part[FoldList[Mod[Plus[#1, #2], 4] &, 0,
a[#] & /@ Range[4^4] /. {F[n_] :> 0, L[n_] :> 1, R[n_] :> -1}],
2 ;; -1] (* Bradley Klee, Aug 07 2015 *)
CROSSREFS
KEYWORD
nonn
AUTHOR
Antti Karttunen, Aug 01 2009
STATUS
approved