A384076 a(n) = pos(M(n)), where M(n) is the n X n circulant matrix with (row 1) = (1,3,5,7, ..., 2n - 1), and pos(M(n)) is the positive part of the determinant of M(n); see A380661.
1, 1, 153, 2668, 250200, 19423560, 2515242520, 404114856640, 84196030473216, 21703670967664000, 6808856052755927808, 2552126898198385479168, 1126590812208410998119424, 578462173661889165983466496, 341831898528862885226121600000
Offset: 1
Keywords
Examples
The rows of M(4) are (1,3,5,7), (7,1,3,5), (5,7,1,3), (3,5,7,1); determinant(M(4)) = -4716; permanent(M(4)) = 2668, so neg(M(4)) = (-2048 - 7384)/2 = -4716 and pos(M(4)) = (-2048 + 7384)/2 = 2668.
Programs
-
Mathematica
z = 19; v[n_] := Table[2 k + 1, {k, 0, n - 1}]; u[n_] := Table[RotateRight[#, k - 1], {k, 1, Length[#]}] &[v[n]]; p = Table[Simplify[Permanent[u[n]]], {n, 1, z}] (* A384074 *) d = Table[Simplify[Det[u[n]]], {n, 1, z}] (* A193678,, with alternating signs *) neg = (d - p)/2 (* A384075 *) pos = (d + p)/2 (* A384076 *)