cp's OEIS Frontend

This is a front-end for the Online Encyclopedia of Integer Sequences, made by Christian Perfect. The idea is to provide OEIS entries in non-ancient HTML, and then to think about how they're presented visually. The source code is on GitHub.

Showing 1-3 of 3 results.

A085656 Number of positive-definite real {0,1} n X n matrices.

Original entry on oeis.org

1, 3, 27, 681, 43369, 6184475, 1688686483, 665444089745
Offset: 1

Views

Author

N. J. A. Sloane, Jul 12 2003

Keywords

Comments

A real matrix M is positive-definite if x M x' > 0 for all nonzero real vectors x. Equivalently, all eigenvalues of M + M' are positive.
M need not be symmetric. For the number of different values of M + M' see A085657. - Max Alekseyev, Dec 13 2005

Examples

			For n = 2 the three matrices are {{{1, 0}, {0, 1}}, {{1, 0}, {1, 1}}, {{1, 1}, {0, 1}}}.

Links

Crossrefs

Cf. A055165, which counts nonsingular {0, 1} matrices and A085506, which counts {-1, 0, 1} matrices with positive eigenvalues.
Cf. A085657, A085658, A086215, A038379 (positive semi-definite matrices), A080858, A083029.

Programs

  • Mathematica
    Table[Count[Tuples[{0, 1}, {n, n}], ?PositiveDefiniteMatrixQ], {n, 4}] (* _Eric W. Weisstein, Jan 03 2021 *)
  • PARI
    { a(n) = M=matrix(n,n,i,j,2*(i==j)); r=0; b(1); r } { b(k) = local(t); if(k>n, t=0; for(i=1,n, for(j=1,i-1, if(M[i,j]==1,t++); )); r+=2^t; return; ); forvec(x=vector(k-1,i,[0,1]), for(i=1,k-1,M[k,i]=M[i,k]=x[i]); if( matdet(vecextract(M,2^k-1,2^k-1),1)>0, b(k+1) ) ) } (Alekseyev)

Extensions

More terms from Max Alekseyev, Dec 13 2005

A114601 Number of n X n symmetric positive definite matrices with 2's on the main diagonal and -1, 0, or 1 elsewhere.

Original entry on oeis.org

1, 3, 23, 393, 13089, 737595, 58969079
Offset: 1

Views

Author

Max Alekseyev, Dec 13 2005, Nov 09 2006

Keywords

Comments

This gives the number of different values of A + A' where A runs through the matrices counted in A086215.

Crossrefs

Cf. A086215, A085657.

Programs

  • PARI
    { b(k) = if(k>m, r++; return); forvec(x=vector(k-1,i,[-1,1]), for(i=1,k-1, M[k,i]=M[i,k]=x[i]); if( matdet(vecextract(M,2^k-1,2^k-1),1)>0, b(k+1) ) ) }
{ a(n) = local(M,r,m); M=matrix(n,n,i,j,2*(i==j)); r=0; m=n; b(1); r }

A127502 Number of n X n positive definite matrices with 1's on the main diagonal and -1's and 0's elsewhere.

Original entry on oeis.org

1, 3, 19, 201, 3001, 55291, 1115003, 21837649, 373215601, 8282131891
Offset: 1

Views

Author

Max Alekseyev, Jan 16 2007

Keywords

Comments

A real matrix M is positive-definite if x M x' > 0 for all nonzero real vectors x. Equivalently, all eigenvalues of M + M' are positive.
M need not be symmetric. For the number of different values of M + M' see A084552.

Examples

			For n = 2 the three matrices are {{{1, 0}, {0, 1}}, {{1, 0}, {-1, 1}}, {{1, -1}, {0, 1}}}.

Crossrefs

Cf. A085656, A085657, A085658, A086215, A038379, A080858, A083029, A127503.

Programs

  • PARI
    { a(n) = M=matrix(n,n,i,j,2*(i==j)); r=0; b(1); r } { b(k) = local(t); if(k> n, t=0; for(i=1,n, for(j=1,i-1, if(M[i,j]==1,t++); )); r+=2^t; return; ); forvec(x=vector(k-1,i,[ -1,0]), for(i=1,k-1,M[k,i]=M[i,k]=x[i]); if( matdet(vecextract(M,2^k-1,2^k-1),1)>0, b(k+1) ) ) } (Alekseyev)
Showing 1-3 of 3 results.

Started in 1964 by Neil J. A. Sloane | Maintained by The OEIS Foundation Inc.

Content is available under The OEIS End-User License Agreement.