A030079 Primes p such that digits of p appear in p^2.

5, 11, 101, 139, 199, 211, 277, 421, 463, 499, 503, 509, 523, 593, 601, 631, 661, 733, 757, 887, 911, 991, 1021, 1091, 1117, 1163, 1171, 1193, 1201, 1277, 1319, 1451, 1523, 1601, 1621, 1627, 1777, 1783, 1999, 2011, 2027, 2111, 2161, 2203

Offset

1,1

Links

Reinhard Zumkeller, Table of n, a(n) for n = 1..10000

Mathematica

ok[p_] := (idp = Union[ IntegerDigits[p]]; Intersection[idp, Union[ IntegerDigits[p^2]]] == idp); Select[ Prime[Range[400]], ok] (* Jean-François Alcover, Nov 30 2011 *)
Select[Prime[Range[400]], SubsetQ[IntegerDigits[#^2], IntegerDigits[#]]&] (* Requires Mathematica version 10 or later *) (* Harvey P. Dale, Dec 23 2017 *)

Prog

(Haskell)
import Data.List (intersect, nub)
a030079 n = a030079_list !! (n-1)
a030079_list = filter f a000040_list where
   f p = pd == pd `intersect` (nub $ show (p^2)) where
       pd = nub $ show p
-- Reinhard Zumkeller, Nov 03 2011

Crossrefs

Sequence in context: A224270 A123025 A053778 * A066596 A199325 A199305

Adjacent sequences: A030076 A030077 A030078 * A030080 A030081 A030082

Keyword

nonn,base,nice

Author

Patrick De Geest

Status

approved