|
|
A037271
|
|
Number of steps to reach a prime under "replace n with concatenation of its prime factors" when applied to n-th composite number, or -1 if no such number exists.
|
|
23
|
|
|
2, 1, 13, 2, 4, 1, 5, 4, 4, 1, 15, 1, 1, 2, 3, 4, 4, 1, 2, 2, 1, 5, 3, 2, 2, 1, 9, 2, 9, 6, 1, 15
(list;
graph;
refs;
listen;
history;
text;
internal format)
|
|
|
OFFSET
|
1,1
|
|
COMMENTS
|
a(33) is presently unknown: starting with 49, no prime has been reached after 110 steps. See A037274 for the latest information.
|
|
LINKS
|
|
|
EXAMPLE
|
Starting with 14 (the seventh composite number) we get 14=2*7, 27=3*3*3, 333=3*3*37, 3337=47*71, 4771=13*367, 13367 is prime; so a(7)=5.
|
|
MATHEMATICA
|
maxComposite = 49; maxIter = 40; concat[n_] := FromDigits[ Flatten[ IntegerDigits /@ Flatten[ Apply[ Table, {#[[1]], {#[[2]]}} & /@ FactorInteger[n], {1}]]]]; composites = Select[ Range[2, maxComposite], ! PrimeQ[#] &]; a[n_] := ( lst = NestWhileList[ concat, composites[[n]], ! PrimeQ[#] &, 1, maxIter]; If[PrimeQ[ Last[lst]], Length[lst] - 1, - 1]); Table[a[n], {n, 1, Length[composites]}] (* Jean-François Alcover, Jul 10 2012 *)
|
|
PROG
|
(Haskell)
a037271 = length . takeWhile ((== 0) . a010051'') .
iterate a037276 . a002808
|
|
CROSSREFS
|
|
|
KEYWORD
|
nonn,nice,base,more,hard
|
|
AUTHOR
|
|
|
STATUS
|
approved
|
|
|
|