login
A133962
Home primes whose homeliness is 4.
12
379, 2347, 2383, 2389, 3331, 3359, 3373, 3719, 3767, 4397, 4759, 7331, 7457, 7523, 7547, 7823, 11251, 13883, 17137, 17317, 19157, 19181, 22367, 22397, 23131, 23173, 23197, 23311, 23593, 23677, 23767, 23911, 29101, 31063, 31123, 31181, 31189
OFFSET
1,1
COMMENTS
Number of terms < 10^n: 0, 0, 1, 16, 146, 1085, 8522, ... .
EXAMPLE
Only (42,74,237,379) -> 379, 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 = 4 - 2; lsu = {}; Do[ If[ lst[[n]] == lst[[n + d]] && lst[[n - 1]] != lst[[n]] && lst[[n]] != lst[[n + d + 1]], AppendTo[lsu, lst[[n]]]], {n, 188004 - d - 1}]; Take[Union@ lsu, 45]
KEYWORD
nonn,base
AUTHOR
Robert G. Wilson v, Sep 30 2007
STATUS
approved