cp's OEIS Frontend

A359559 a(n) is the determinant of an n X n Hermitian Toeplitz matrix whose first row consists of 1, 2*i, ..., n*i, where i denotes the imaginary unit.

%I A359559 #22 Jan 25 2023 18:37:08
%S A359559 1,1,-3,-16,-36,-40,20,184,400,432,-112,-1472,-3136,-3328,576,9856,
%T A359559 20736,21760,-2816,-59392,-123904,-129024,13312,333824,692224,716800,
%U A359559 -61440,-1785856,-3686400,-3801088,278528,9207808,18939904,19464192,-1245184,-46137344,-94633984
%N A359559 a(n) is the determinant of an n X n Hermitian Toeplitz matrix whose first row consists of 1, 2*i, ..., n*i, where i denotes the imaginary unit.
%H A359559 Wikipedia, <a href="https://en.wikipedia.org/wiki/Toeplitz_matrix">Toeplitz Matrix</a>
%F A359559 A359614(n) <= a(n) <= A359615(n).
%F A359559 Conjectured formulas: (Start)
%F A359559 O.g.f.: (1 - 5*x + 9*x^2 - 12*x^3 + 10*x^4 - 4*x^5)/(1 - 2*x + 2*x^2)^3.
%F A359559 a(n) = 6*a(n-1) - 18*a(n-2) + 32*a(n-3) - 36*a(n-4) + 24*a(n-5) - 8*a(n-6) for n > 5.
%F A359559 E.g.f.: (2 + exp(x)*((1 + x)*(2 + x)*cos(x) - (1 + x + x^2)*sin(x)))/4. (End)
%e A359559 a(3) = -16:
%e A359559   [   1,  2*i, 3*i;
%e A359559    -2*i,    1, 2*i;
%e A359559    -3*i, -2*i,   1 ]
%t A359559 Join[{1},Table[Det[ToeplitzMatrix[Join[{1},I Range[2,n]]]],{n,36}]]
%o A359559 (PARI) a(n) = matdet(matrix(n, n, i, j, if (i==j, 1, if (i<j, I*(j-i+1), I*(j-i-1))))); \\ _Michel Marcus_, Jan 20 2023
%o A359559 (Python)
%o A359559 from sympy import Matrix, I
%o A359559 def A359559(n): return Matrix(n,n,[i-j+(1 if i>j else -1) if i!=j else I for i in range(n) for j in range(n)]).det()*(1,-I,-1,I)[n&3] # _Chai Wah Wu_, Jan 25 2023
%Y A359559 Cf. A001792 (symmetric Toeplitz matrix), A143182.
%Y A359559 Cf. A359560 (permanent), A359561, A359562.
%Y A359559 Cf. A359614 (minimal), A359615 (maximal).
%K A359559 sign
%O A359559 0,3
%A A359559 _Stefano Spezia_, Jan 06 2023