login
A214641
Indices of a in the sequence closed under a -> abc, b -> ab, c -> b . Start with a.
4
1, 4, 7, 10, 12, 14, 17, 20, 23, 25, 28, 30, 33, 36, 39, 41, 43, 46, 49, 52, 54, 57, 60, 63, 65, 68, 71, 74, 76, 78, 81, 84, 87, 89, 92, 94, 97, 100, 103, 105, 107, 110, 113, 116, 118, 121, 124, 127, 129, 131, 134, 137, 140, 142, 145, 148, 151, 153, 155, 158
OFFSET
1,2
COMMENTS
Indices of 1 in A214640.
LINKS
EXAMPLE
a -> abc -> abcabb -> abcabbabcabab -> abcabbabcabababcabbabcababcab -> ...
MAPLE
S:= proc(n) option remember; `if`(n=0, [1], subs(x=[1, 2, 3][],
y=[1, 2][], z=2, subs(1=x, 2=y, 3=z, S(n-1))))
end:
a:= proc(n) option remember; local i, k; i:=1;
for k from 1+`if`(n=1, 0, a(n-1)) do
while nops(S(i))<k do i:=i+1 od;
if S(i)[k] = 1 then return k fi
od
end:
seq(a(n), n=1..100); # Alois P. Heinz, Jul 28 2012
MATHEMATICA
SubstitutionSystem[{"a" -> {"a", "b", "c"}, "b" -> {"a", "b"}, "c" -> {"b"}}, {"a"}, 7] // Last // Position[#, "a"]& // Flatten (* Jean-François Alcover, Nov 04 2020 *)
CROSSREFS
KEYWORD
nonn
AUTHOR
Philippe Deléham, Jul 23 2012
STATUS
approved