|
|
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
(list;
graph;
refs;
listen;
history;
text;
internal format)
|
|
|
OFFSET
|
1,1
|
|
COMMENTS
|
Number of terms < 10^n: 0, 0, 0, 2, 31, 341, 2695, ....
|
|
LINKS
|
|
|
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
|
Cf. A037274, A118756, A133957, A133959, A133961, A133963, A133967, A133969, A133971, A133973, A133975, A133977, A133979.
|
|
KEYWORD
|
nonn,base
|
|
AUTHOR
|
|
|
STATUS
|
approved
|
|
|
|