%我

%第13、3111333113331311133313113113113131133131331313331113313113，

%电话：311331331311113131333113113133313131133311331331131，

%U 11331113331131131133133311133131131311333131133313113331133131

%N个仅由数字1和3组成的素数，出现频率相等。

%有19个数字对可以产生这样的素数。（1，0），（7，0），（1，3），（1，4），（1，6），（1，7），（1，9），（2，3），（2，9），（3，5），（3，7），（3，8），（4，7），（4，9），（5，9），（6，7），（7，9），（8，9）。-更正人：Robert Israel，2018年7月10日

%C位数为偶数且不能被3整除。--罗伯特·以色列，2018年7月9日

%H Robert Israel，<a href=“/A087511/b087511.txt”>n，a（n）表，n=1..10000</a>

%p排序（select（isprime，[seq（seq（（10^（2*d）-1）/9+2*add（10^i，i=s），s=combinat:-choose（[$0..（2*d-1）]，d）），d=[1,2,4,5,7,8,10]））；\_罗伯特·以色列，2018年7月9日

%t Union[FromDigits/@Select[Flatten[Table[Tuples[{1,3}，k]，{k，10}]，1]，PrimeQ[fromdights[#]]&&Count[#，1]==Count[#，3]&]]（*jayantabasu_年5月19日*）

%o（PARI）d1=1；d2=3；k=0；a=矢量（100）；对于（n=13000，B=二进制（n）；L=长度（B）；s=和（j=1，长度（B），B[j]）；如果（L%2==0&s==L/2，C=向量（L，n，（d2-d1）*B[n]+d1）；p=次方（Pol（C），x，10）；如果（isprime（p），如果（k<100，k++；a[k]=p））；D=矢量（L，n，d2-（d2-d1）*B[n]）；q=次方（Pol（D），x，10）；if（isprime（q），if（k<100，k++；a[k]=q））；）；a=向量（k，n，a[n]）；向量排序（a）

%Y比照A087510、A087512、A087513。

%K基，非N

%O 1,1号

%A·保罗·D·汉娜和阿玛纳特·穆尔蒂，2003年9月11日