%I#44 2023年9月22日05:15:14

%S 1,2,3,6,26,9621011063759125777407242388438281673611，

%电话：9983885163009174645253250203888756503422457140173335164762941，

%电话：3131546394233724135323108502235010653736573217903852593844715806535267031861418150648418565342442057063085601738871276081763499377227299

%N配对{1..2n}的方法的数量，所以每对的和是素数。

%H Martin Fuller，<a href=“/A00341/b000341.txt”>n的表，a（n）表示n=1..36</a>

%H L.E.Greenfield和S.J.Greenfield<a href=“https://cs.uwaterloo.ca/journals/JIS/green.html“>与Bertrand假设相关的组合数论的一些问题，J.Integer Sequences，1998，#98.1.2。

%F a（n）=永久（m），其中n X n矩阵m由m（i，j）=1或0定义，取决于2i+2j-1分别是素数还是复合数_T.D.Noe_，2007年2月10日

%e对于n=4，有6种方法对{1，2，3，4，5，6，7，8}进行配对，以便每对求和为素数：

%e 1+2、3+4、5+8、6+7

%e 1+2、3+8、4+7、5+6

%e 1+4、2+3、5+8、6+7

%e 1+4、2+5、3+8、6+7

%e 1+6、2+3、4+7、5+8

%e 1+6、2+5、3+8、4+7

%e因此a（4）=6_迈克尔·波特（Michael B.Porter），2016年7月19日

%p f：=proc（n）局部M；

%p M：=矩阵（n，n，（i，j）->`if`（i素数（2*i+2*j-1），1,0））；

%p线性代数：-永久（M）

%p端程序：

%p映射（f，[1..20]）；#_罗伯特·伊斯雷尔（Robert Israel），2016年7月19日

%t a[n_]：=永久[Array[Boole[PrimeQ[2*#1+2*#2-1]]&，{n，n}]；表[an=a[n]；打印[an]；an，{n，1，20}]（*_Jean-François Alcover_，2011年10月21日，在_T.D.Noe_之后，2016年2月7日更新*）

%o（PARI）permRWNb（a）=n=材料尺寸（a）[1]；如果（n==1，则返回（a[1,1]）；sg=1；nc=0；in=矢量v（n）；x=英寸；x=a[，n]-总和（j=1，n，a[，j]）/2；p=触头（i=1，n，x[i]）；对于（k=1,2^（n-1）-1，sg=-sg；j=估值（k，2）+1；z=1-2*英寸[j]；单位[j]+=z；nc+=z；x+=z*a[，j]；p+=触头（i=1，n，x[i]，sg）；返回（2*（2*）（n%2）-1）*p）

%o表示（n=1,24，a=矩阵（n，n，i，j，isprime（2*（i+j）-1））；print1（permRWNb（a）“，”）\\Herman Jamke（hermanjamke（AT）fastmail.fm），2007年5月13日

%o（PARI）a（n）=矩阵（n，n，i，j，i素数（2*（i+j）-1））；\\_Martin Fuller，2023年9月22日

%Y参考A005326、A009692。

%K nonn很好

%O 1,2号机组

%A S.J.Greenfield（greenfiel（AT）math.rutgers.edu）

%E来自_David W.Wilson的更多条款_

%E更多条款，来自T.D.Noe_，2007年2月10日

%E来自Herman Jamke（hermanjamke（AT）fastmail.fm）的更多术语，2007年5月13日

%E更多条款，来自Sean A.Irvine_，2010年11月14日