A133965 Home primes whose homeliness is greater than 5.

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, ....

Links

Table of n, a(n) for n=1..34.

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.

Sequence in context: A116352 A153746 A243217 * A031611 A133966 A043671

Adjacent sequences: A133962 A133963 A133964 * A133966 A133967 A133968

Keyword

nonn,base

Author

Robert G. Wilson v, Sep 30 2007

Status

approved