login
A133965
Home primes whose homeliness is greater than 5.
12
3389, 3761, 13367, 23251, 23761, 23929, 31193, 31397, 33797, 36389, 37199, 37463, 37547, 37607, 37643, 37717, 37853, 37907, 37951, 37967, 41863, 41887, 41941, 73331, 74759, 74771, 77269, 77711, 77773, 77797, 83383, 112909, 113779, 114773
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]
KEYWORD
nonn,base
AUTHOR
Robert G. Wilson v, Sep 30 2007
STATUS
approved