|
%I #36 Dec 07 2019 12:18:28
%S 0,12,35,56,99,132,195,240,323,380,483,552,675,756,899,992,1155,1260,
%T 1443,1560,1763,1892,2115,2256,2499,2652,2915,3080,3363,3540,3843,
%U 4032,4355,4556,4899,5112,5475,5700,6083,6320,6723,6972,7395,7656,8099,8372,8835,9120,9603,9900,10403
%N Number of 2 X 2 matrices with all elements in {0,...,n} with permanent = determinant^n.
%H Indranil Ghosh, <a href="/A280364/b280364.txt">Table of n, a(n) for n = 0..995</a>
%F a(n-1) = (((-2*(n-1)-1)*a(n-1))/(2*(n-1)-1)) + 8*(n-1)^2 + 10*(n-1) + 3 for n>=2, with a(0)=0.
%F Conjectures from _Colin Barker_, Jan 01 2017: (Start)
%F a(n) = 4*n^2 + 8*n + 3 for n>0 and even.
%F a(n) = 4*n^2 + 6*n + 2 for n odd.
%F a(n) = a(n-1) + 2*a(n-2) - 2*a(n-3) - a(n-4) + a(n-5) for n>5.
%F G.f.: x*(12 + 23*x - 3*x^2 - 3*x^3 + 3*x^4) / ((1 - x)^3*(1 + x)^2).
%F (End)
%e For n=2, the matrices are [0,0,0,0], [0,0,0,1], [0,0,0,2], [0,0,1,0], [0,0,1,1], [0,0,1,2], [0,0,2,0], [0,0,2,1], [0,0,2,2], [0,1,0,0], [0,1,0,1], [0,1,0,2], [0,1,1,0], [0,1,1,1], [0,1,1,2], [0,2,0,0], [0,2,0,1], [0,2,0,2], [1,0,0,0], [1,0,0,1], [1,0,1,0], [1,0,1,1], [1,0,2,0], [1,0,2,1], [1,1,0,0], [1,1,0,1], [1,1,1,0], [1,2,0,0], [1,2,0,1], [2,0,0,0], [2,0,1,0], [2,0,2,0], [2,1,0,0], [2,1,1,0], [2,2,0,0].
%e Here each of these matrices M is defined as M=[a,b,c,d], where a=M[1][1], b=M[1][2], c=M[2][1], d=M[2][2]. There are 35 possibilities. So for n=2, a(n)=35.
%t CoefficientList[Series[x (12 + 23 x - 3 x^2 - 3 x^3 + 3 x^4)/((1 - x)^3*(1 + x)^2), {x, 0, 50}], x] (* _Michael De Vlieger_, Jan 01 2017 *)
%o (Python)
%o def t(n):
%o s=0
%o for a in range(0,n+1):
%o for b in range(0,n+1):
%o for c in range(0,n+1):
%o for d in range(0,n+1):
%o if (a*d-b*c)**n==(a*d+b*c):
%o s+=1
%o return s
%o for i in range(0,51):
%o print str(i)+" "+str(t(i))
%o (Python)
%o def a(n):
%o if n==2:
%o return 12
%o return (((-2*(n-1)-1)*a(n-1))/(2*(n-1)-1))+8*(n-1)**2+10*(n-1)+3
%o for i in range(2,1001):
%o print str(i)+" "+str(a(i))
%Y Cf. A280344 (Number of 2 X 2 matrices with all elements in {0,...,n} with determinant = permanent^n).
%K nonn
%O 0,2
%A _Indranil Ghosh_, Jan 01 2017
|
