%I#14 2022年1月11日22:00:06
%S 23232233323223222232332322322333322322233332322322322323232323223232322322333,
%电话:2232332233232223333232233322323232322323233232233223223232232323,
%电话:323222333332323323223323232232232333232232222223333333
%N只由数字2和数字3组成的素数,出现频率相同。
%C有18个数字对可以产生这样的素数。(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 Michael S.Branicky,n的表,a(n)表示n=1..10000</a>
%o(Python)
%o来自sympy import isprime
%o来自sympy.utilities.iterables导入multiset_permutations
%o定义auptodigs(最大数字):
%o alst=[]
%o表示范围(2,最大数字+1,2)中的d:
%o ms=“2”*(d//2)+“3”*(d//2-1)
%o表示多置换(ms,d-1)中的p:
%o t=int(“.join(p)+”3“)
%o如果是质数(t):
%o另外附加(t)
%o返回alst
%o打印(auptodigs(10))#_Michael S.Branicky_,2022年1月11日
%o(PARI)seq(n,d1=2,d2=3)={my(L=List());for(d=1,oo,forperm(向量(2*d,i,if(i<=d,d1,d2)),v,my(q=fromdigits(Vec(v))))
%Y参见A087510、A087511和A087514。
%K基,nonn
%O 1,1号机组
%A·保罗·D·汉纳和阿玛娜斯·穆尔西,2003年9月12日
%E偏移改为1,a(19)由_Georg Fischer修正,2022年1月11日
|