|
%I #9 Aug 21 2016 11:47:33
%S 44994949,49444999,49494499,49499449,94449499,4449449999,4449494999,
%T 4494494999,4494499949,4494949499,4494994499,4499449949,4499494949,
%U 4499944949,4944499949,4949444999,4994449949,9444499949,9444994499
%N Primes consisting only of digits 4 and 9 occurring with equal frequency.
%C There are 18 digit pairs which can produce such primes. (1,0),(7,0),(1,3),(1,4),(1,6),(1,7),(1,9),(2,3),(2,9),(3,4),(3,5),(3,7),(3,8),(4,7),(4,9),(5,9),(6,7),(7,9).
%H Harvey P. Dale, <a href="/A087534/b087534.txt">Table of n, a(n) for n = 0..10000</a>
%t Sort[Flatten[Table[Select[FromDigits/@Permutations[PadRight[{},2n,{4,9}]],PrimeQ],{n,5}]]] (* _Harvey P. Dale_, Aug 21 2016 *)
%o (PARI) /* Primes consisting only of digits x and y, occurring with equal frequency. */ x=4; y=9; d1=x; d2=y; k=0; a=vector(100); for(n=1,3000,B=binary(n); L=length(B); s=sum(j=1,length(B),B[j]); if(L%2==0 & s==L/2, C=vector(L,n,(d2-d1)*B[n]+d1); p=subst(Pol(C),x,10); if(isprime(p),if(k<100,k++; a[k]=p)); D=vector(L,n,d2-(d2-d1)*B[n]); q=subst(Pol(D),x,10); if(isprime(q ),if(k<100,k++; a[k]=q))); ); a=vector(k,n,a[n]); vecsort(a)
%Y Cf. A087510, A087511, A087533.
%K base,nonn
%O 0,1
%A _Paul D. Hanna_ and _Amarnath Murthy_, Sep 12 2003
%E Edited by _Charles R Greathouse IV_, Oct 28 2009
|
