login
a(n) is the determinant of a symmetric Toeplitz matrix M(n) whose first row consists of prime(1), prime(2), ..., prime(n).
3

%I #25 Jan 31 2023 05:35:07

%S 1,2,-5,12,-19,-22,1143,-4284,14265,-46726,-84405,1306096,32312445,

%T 522174906,4105967871,5135940112,-642055973735,-2832632334858,

%U 14310549077571,380891148658140,4888186898996125,-49513565563840210,383405170118692791,-2517836083641473036,-3043377347606882055

%N a(n) is the determinant of a symmetric Toeplitz matrix M(n) whose first row consists of prime(1), prime(2), ..., prime(n).

%C Conjecture: abs(a(n)) is prime only for n = 1, 2, and 4.

%H Mathematics Stack Exchange, <a href="https://math.stackexchange.com/questions/3736861/determinant-of-a-toeplitz-matrix">Determinant of a Toeplitz matrix</a>

%H Wikipedia, <a href="https://en.wikipedia.org/wiki/Toeplitz_matrix">Toeplitz Matrix</a>

%F A350955(n) <= a(n) <= A350956(n).

%e For n = 1 the matrix M(1) is

%e 2

%e with determinant a(1) = 2.

%e For n = 2 the matrix M(2) is

%e 2, 3

%e 3, 2

%e with determinant a(2) = -5.

%e For n = 3 the matrix M(3) is

%e 2, 3, 5

%e 3, 2, 3

%e 5, 3, 2

%e with determinant a(3) = 12.

%p A356490 := proc(n)

%p local T,c ;

%p if n =0 then

%p return 1 ;

%p end if;

%p T := LinearAlgebra[ToeplitzMatrix]([seq(ithprime(c),c=1..n)],n,symmetric) ;

%p LinearAlgebra[Determinant](T) ;

%p end proc:

%p seq(A356490(n),n=0..15) ; # _R. J. Mathar_, Jan 31 2023

%t k[i_]:=Prime[i]; M[ n_]:=ToeplitzMatrix[Array[k, n]]; a[n_]:=Det[M[n]]; Join[{1},Table[a[n],{n,24}]]

%o (PARI) a(n) = matdet(apply(prime, matrix(n,n,i,j,abs(i-j)+1))); \\ _Michel Marcus_, Aug 12 2022

%o (Python)

%o from sympy import Matrix, prime

%o def A356490(n): return Matrix(n,n,[prime(abs(i-j)+1) for i in range(n) for j in range(n)]).det() # _Chai Wah Wu_, Aug 12 2022

%Y Cf. A005843 (trace of M(n)), A309131 (k-superdiagonal sum of M(n)), A356483 (hafnian of M(2*n)), A356491 (permanent of M(n)).

%Y Cf. A350955, A350956.

%K sign

%O 0,2

%A _Stefano Spezia_, Aug 09 2022