A133963 Home primes whose homeliness is greater than 4.

773, 1129, 3137, 3251, 3389, 3593, 3761, 3797, 7541, 7853, 13367, 23167, 23251, 23557, 23719, 23743, 23761, 23773, 23929, 23971, 31193, 31397, 31973, 33191, 33331, 33769, 33797, 33863, 36389, 37199, 37307, 37463, 37547, 37573, 37607

Offset

1,1

Comments

Number of terms < 10^n: 0, 0, 1, 10, 75, 758, 5994, ....

Links

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

Example

(10, 25, 55, 511, 773) -> 773, 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]], AppendTo[ lsu, lst[[n]] ]], {n, 188004 - d - 1}]; Take[ Union@ lsu, 40]

Crossrefs

Cf. A037274, A118756, A133957, A133959, A133961, A133965, A133967, A133969, A133971, A133973, A133975, A133977, A133979.

Sequence in context: A285755 A264255 A240843 * A133964 A033919 A055521

Adjacent sequences: A133960 A133961 A133962 * A133964 A133965 A133966

Keyword

nonn,base

Author

Robert G. Wilson v, Sep 30 2007

Status

approved