OFFSET
1,1
COMMENTS
Number of terms < 10^n: 0, 0, 0, 2, 31, 341, 2695, ....
EXAMPLE
(118, 259, 737, 801, 1167, 3389) -> 3389, etc.
MATHEMATICA
lst = {}; f[n_] := FromDigits@ Flatten[ IntegerDigits@ Table[ #[[1]], {#[[2]]}] & /@ FactorInteger@n, 2]; h[n_] := NestWhileList[f@# &, n, !PrimeQ@# &, 1, 28]; Do[p = h[n][[ -1]]; If[ PrimeQ@p && p < 10^7 && p != n, Print[{n, p}]; AppendTo[lst, p]], {n, 2, 1000}];
d = 6 - 2; lsu = {}; Do[ If[ lst[[n]] == lst[[n + d]], AppendTo[ lsu, lst[[n]] ]], {n, 188004 - d - 1}]; Take[ Union@ lsu, 35]
CROSSREFS
KEYWORD
nonn,base
AUTHOR
Robert G. Wilson v, Sep 30 2007
STATUS
approved