OFFSET
0,3
COMMENTS
This sequence is infinite. - Robert G. Wilson v, Dec 09 2014
LINKS
Robert G. Wilson v, Table of n, a(n) for n = 0..46
EXAMPLE
a(0) = 0;
a(1) = 1 -> 10 / 1 = 10;
a(2) = 2 -> 210 / 2 = 105;
a(3) = 3 -> 3210 / 3 = 1070;
Now we cannot use 4 as the next term because 43210 / 4 = 21605 / 2.
a(4) = 5 -> 32105 / 5 = 6421; etc.
MAPLE
with(numtheory); P:=proc(q) local a, k, n; print(0); print(1); a:=10;
for n from 2 to q do if type((n*10^(1+ilog10(a))+a)/n, integer)
then a:=n*10^(1+ilog10(a))+a; print(n);
fi; od; end: P(10^9);
MATHEMATICA
f[lst_List] := Block[{k = lst[[-1]] + 1, id = FromDigits@ Flatten@ IntegerDigits@ Reverse@ lst}, While[ Mod[ id, k] > 0, k++]; Append[lst, k]]; Nest[f, {0}, 36] (* or *)
f[lst_List] := Block[{mn = lst[[-1]], id = FromDigits@ Flatten@ IntegerDigits@ Reverse@ lst}, d = Divisors@ id; Append[lst, Min@ Select[d, # > mn &]]]; Nest[f, {0, 1}, 36] (* Robert G. Wilson v, Dec 08 2014 *)
CROSSREFS
KEYWORD
nonn,base
AUTHOR
Paolo P. Lava, Nov 27 2014
EXTENSIONS
a(37)-a(40) from Robert G. Wilson v, Dec 08 2014
STATUS
approved