cp's OEIS Frontend

A079652 Prime numbers using only the curved digits 0, 3, 6, 8 and 9.

%I A079652 #20 Mar 29 2023 08:58:52
%S A079652 3,83,89,383,389,683,809,839,863,883,983,3083,3089,3389,3803,3833,
%T A079652 3863,3889,3989,6089,6389,6689,6803,6833,6863,6869,6883,6899,6983,
%U A079652 8009,8039,8069,8089,8093,8363,8369,8389,8609,8663,8669,8689,8693,8699,8803,8839
%N A079652 Prime numbers using only the curved digits 0, 3, 6, 8 and 9.
%C A079652 Intersection of A000040 and A072960. - _K. D. Bajpai_, Sep 01 2014
%H A079652 K. D. Bajpai, <a href="/A079652/b079652.txt">Table of n, a(n) for n = 1..11740</a>
%H A079652 Chris K. Caldwell and G. L. Honaker, Jr., <a href="https://primes.utm.edu/curios/page.php?short=30689">30689</a>, Prime Curios!
%H A079652 Chris K. Caldwell and G. L. Honaker, Jr., <a href="https://primes.utm.edu/curios/page.php?short=90863">90863</a>, Prime Curios!
%p A079652 N:= 4: # to get all terms with up to N digits
%p A079652 Digs:= {0,3,6,8,9}:
%p A079652 A:= NULL:
%p A079652 for d from 1 to N do
%p A079652   C:= combinat[cartprod]([Digs minus {0},Digs $(d-1)]);
%p A079652   while not C[finished] do
%p A079652     L:= C[nextvalue]();
%p A079652     x:= add(L[i]*10^(d-i),i=1..d);
%p A079652     if isprime(x) then A:= A,x fi
%p A079652   od
%p A079652 od:
%p A079652 A; # _Robert Israel_, Aug 31 2014
%t A079652 Select[ Range[8850], PrimeQ[ # ] && Union[ Join[ IntegerDigits[ # ], {0, 3, 6, 8, 9}]] == {0, 3, 6, 8, 9} &]
%t A079652 Select[Prime[Range[5000]], Intersection[IntegerDigits[#], {1, 2, 4, 5, 7}] == {} &] (* _K. D. Bajpai_, Sep 01 2014 *)
%t A079652 Select[FromDigits/@Tuples[{0,3,6,8,9},4],PrimeQ] (* _Harvey P. Dale_, Sep 05 2022 *)
%Y A079652 Cf. A072960, A028374.
%Y A079652 Cf. A034470.
%K A079652 base,nonn
%O A079652 1,1
%A A079652 _Shyam Sunder Gupta_, Jan 23 2003