login
The OEIS Foundation is supported by donations from users of the OEIS and by a grant from the Simons Foundation.

 

Logo


Hints
(Greetings from The On-Line Encyclopedia of Integer Sequences!)
A331346 Primes using all the square digits {0, 1, 4, 9} and no others. 1
1049, 1409, 4019, 4091, 9041, 10499, 10949, 14009, 49019, 49109, 90149, 90401, 94109, 99041, 99401, 100049, 101149, 101419, 101449, 104009, 104119, 104149, 104491, 104911, 104999, 109049, 109141, 109441, 110419, 110491, 111049, 111409, 114901, 140009, 140191, 140419 (list; graph; refs; listen; history; text; internal format)
OFFSET

1,1

COMMENTS

Subsequence of A061246.

LINKS

Robert Israel, Table of n, a(n) for n = 1..10000

EXAMPLE

a(1) = 1049 is prime containing all the square digits (0, 1, 4, 9) and no others.

a(2) = 1409 is prime containing all the square digits (0, 1, 4, 9) and no others.

MAPLE

f:= proc(n) local L, x;

  L:= convert(n, base, 4);

  if convert(L, set) <> {0, 1, 2, 3} then return NULL fi;

  L:= subs(2=4, 3=9, L);

  x:= add(L[i]*10^(i-1), i=1..nops(L));

  if isprime(x) then x else NULL fi

end proc:

map(f, [$4^3..4^6]); # Robert Israel, Jan 16 2020

MATHEMATICA

Select[FromDigits /@ Tuples[{0, 1, 4, 9}, 6], PrimeQ[#] && Union[IntegerDigits[#]] == {0, 1, 4, 9} &]

PROG

(MAGMA) [p:p in PrimesUpTo(150000)|Set(Intseq(p)) eq {0, 1, 4, 9}]; // Marius A. Burtea, Jan 14 2020

CROSSREFS

Cf. A000040, A061246, A066591, A323387.

Sequence in context: A251375 A281033 A337960 * A164771 A031600 A251272

Adjacent sequences:  A331343 A331344 A331345 * A331347 A331348 A331349

KEYWORD

nonn,base

AUTHOR

K. D. Bajpai, Jan 14 2020

STATUS

approved

Lookup | Welcome | Wiki | Register | Music | Plot 2 | Demos | Index | Browse | More | WebCam
Contribute new seq. or comment | Format | Style Sheet | Transforms | Superseeker | Recent
The OEIS Community | Maintained by The OEIS Foundation Inc.

License Agreements, Terms of Use, Privacy Policy. .

Last modified March 2 09:23 EST 2021. Contains 341746 sequences. (Running on oeis4.)