%I#36 2021年10月1日13:17:19
%S 1,2,1234334997125154411574945708172424550598849,
%电话1282759836215548737
%行列式等于其永久数的N×N个二元矩阵的个数。
%C下边界为A088672。
%C类似于A145675和A145676。
%H Christopher Culter,<a href=“https://github.com/Culter/PermDet网站“>用于计算大型项的C++代码</a>
%H数学堆栈交换,<a href=“http://math.stackexchange.com/q/1707432/87023“>使det(A)=perm(A)的n X n个二进制矩阵A的个数是多少</a>
%F a(n)<=2^(n^2),n<=1时相等。
%e a(2)等于12,因为正好有十二个2X2二元矩阵的行列式等于它们的永久值;这些矩阵是:
%e |0 0||1 0||0 1||1 1||0 0||1 0||0 0|1 0|
%e |0 0||0 0||0 0||0 0 ||0 0| |1 0||1 0| |0 1||0 1|
%e、。
%e |0 1||1 1||0 0||1 0|
%e |0 1||0 1||1 1||11 1|
%t求和[KroneckerDelta[Det[Array[Mod[Floor[k/(2^(n*(#1-1)+#2-1))],2]&,{n,n}]],永久[Array[Mod[Floor[k/(2^(n*(#1-1)+#2-1))],2]&,{n,n}]],{k,0,(2^(n^2))-1}]
%o(Python)
%o来自itertools导入产品
%o从sympy导入矩阵
%o定义A192892(n):如果n==0,则返回1,否则求和(1代表乘积([0,1],repeat=n**2)如果(λx:x.det()==x.per())(矩阵(n,n,m))#_Chai Wah Wu_,2021年10月1日
%Y参见A046747、A087983、A089472、A08947、A145675、A145766。
%K难,多,不,好
%0、2
%A _John M.Campbell,2011年7月11日
%2016年4月13日,《基督教文化》中的E a(0)=1和a(5)-a(8)
%E定义和示例由Harvey P.Dale_稍作修改,2017年2月24日
|