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

Hints
(Greetings from The On-Line Encyclopedia of Integer Sequences!)
A133958 Home primes whose homeliness is 2. 12
37, 227, 229, 233, 241, 257, 271, 277, 283, 313, 317, 331, 347, 353, 359, 367, 383, 397, 523, 541, 557, 577, 719, 743, 761, 797, 1117, 1171, 1361, 1367, 1373, 1723, 1741, 1747, 1753, 1759, 1783, 1789, 1973, 1979, 1997, 2113, 2131, 2137, 2179, 2213, 2239 (list; graph; refs; listen; history; internal format)
OFFSET

1,1

COMMENTS

Number of terms < 10^n: 0, 1, 26, 182, 1428, 10395, 84164, ....

EXAMPLE

Only {21 & 37} -> 37, 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 = 2 - 2; lsu = {}; Do[ If[ lst[[n]] == lst[[n + d]] && lst[[n-1]] != lst[[n]] && lst[[n]] != lst[[n + d + 1]], AppendTo[lsu, lst[[n]]]], {n, 188004 - d - 1}]; Take[ Union@ lsu, 50]

CROSSREFS

Cf. A037274, A118756, A133957, A133956, A000014, A133960, A133962, A133964, A133966, A133968, A133970, A133972, A133974, A133976, A133978.

Sequence in context: A141984 A052166 A142010 * A088544 A051463 A142445

Adjacent sequences:  A133955 A133956 A133957 * A133959 A133960 A133961

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 09:26 EST 2012. Contains 205753 sequences.