login
This site is supported by donations to The OEIS Foundation.
Logo

Hints
(Greetings from The On-Line Encyclopedia of Integer Sequences!)
A133963 Home primes whose homeliness is greater than 4. 12
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 (list; graph; refs; listen; history; internal format)
OFFSET

1,1

COMMENTS

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

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: A205357 A045078 A114543 * A133964 A033919 A055521

Adjacent sequences:  A133960 A133961 A133962 * A133964 A133965 A133966

KEYWORD

nonn

AUTHOR

Robert G. Wilson v (rgwv(AT)rgwv.com), Sep 30 2007

Lookup | Welcome | Wiki | Register | Music | Plot 2 | Demos | Index | Browse | More | WebCam
Contribute new seq. or comment | Format | Transforms | Puzzles | Hot | Classics
Recent Additions | More pages | Superseeker | Maintained by The OEIS Foundation Inc.

Content is available under The OEIS End-User License Agreement .

Last modified February 15 02:50 EST 2012. Contains 205694 sequences.