|
%I #30 Dec 10 2015 14:11:26
%S 1,2,11,144,4149,251622,31340799,7913773980,4024015413705,
%T 4106387069191890,8395359475529822355,34357677843892688699400,
%U 281336437060919094044274525,4608419756389534634440592965950,150992374805715685629827976712244775
%N Determinant of the n-th principal submatrix of A204242.
%H Colin Barker, <a href="/A204243/b204243.txt">Table of n, a(n) for n = 1..80</a>
%H Wikipedia, <a href="https://en.wikipedia.org/wiki/Matrix_determinant_lemma">Matrix determinant lemma</a>
%F a(n) = (1 - Sum_{k=2..n} 1/(2^k-1)) * Product_{k=2..n} (2^k-1) = 2*A005329(n) - A203011(n). - _Robert Israel_, Nov 30 2015
%p f:= n -> (1 - add(1/(2^i-1),i=2..n))*mul(2^i-1,i=2..n):
%p seq(f(n),n=1..30); # _Robert Israel_, Nov 30 2015
%t f[i_, j_] := 0; f[1, j_] := 1; f[i_, 1] := 1; f[i_, i_] := 2^i - 1;
%t m[n_] := Table[f[i, j], {i, 1, n}, {j, 1, n}]
%t TableForm[m[8]] (* 8x8 principal submatrix *)
%t Flatten[Table[f[i, n + 1 - i],
%t {n, 1, 12}, {i, 1, n}]] (* A204242 *)
%t Table[Det[m[n]], {n, 1, 15}] (* A204243 *)
%t Permanent[m_] :=
%t With[{a = Array[x, Length[m]]},
%t Coefficient[Times @@ (m.a), Times @@ a]];
%t Table[Permanent[m[n]], {n, 1, 15}] (* A203011 *)
%o (PARI) vector(20, n, matdet(matrix(n, n, i, j, if(i==1, 1, if(j==1, 1, if(i==j, 2^i-1)))))) \\ _Colin Barker_, Nov 27 2015
%Y Cf. A005329, A203011, A204238, A204242.
%K nonn
%O 1,2
%A _Clark Kimberling_, Jan 13 2012
|
