OFFSET
1,1
LINKS
Reinhard Zumkeller, Table of n, a(n) for n = 1..1000
FORMULA
MATHEMATICA
Select[Prime[Range[82000]], Union[IntegerDigits[#]]== Union[ IntegerDigits [#^2]]&] (* Harvey P. Dale, Aug 12 2011 *)
PROG
(PARI) isA030091(n)=isprime(n)&&Set(Vec(Str(n)))==Set(Vec(Str(n^2))) \\ Charles R Greathouse IV, Jun 28 2011
(Haskell)
import Data.List (nub, sort)
import Data.Function (on)
a030091 n = a030091_list !! (n-1)
a030091_list =
filter (\x -> ((==) `on` (nub . sort . show)) x (x^2)) a000040_list
-- Reinhard Zumkeller, Aug 11 2011
CROSSREFS
KEYWORD
nonn,base,nice
AUTHOR
STATUS
approved