login
The OEIS is supported by the many generous donors to the OEIS Foundation.

 

Logo
Hints
(Greetings from The On-Line Encyclopedia of Integer Sequences!)
A089482 Number of real {0,1}-matrices having permanent = 1. 8
1, 1, 6, 150, 13032, 3513720, 2722682160, 5739447495600, 31598877919109760, 440333998013384657280, 15150599165671354541318400, 1261508968034974650352062240000, 250009928097136435131869478983500800, 116299581308873767293693697630883742796800 (list; graph; refs; listen; history; text; internal format)
OFFSET
0,3
COMMENTS
a(6) from Gordon F. Royle.
The following is Max Alekseyev's proof of the formula: Suppose that we have a (0,1)-matrix M with permanent equal to 1. Then in M there is a unique set of n elements, each equal to 1, whose product makes the permanent equal 1. Permute the columns of M so that these n elements become arranged along the main diagonal, denote the resulting matrix by M'. It is clear that each M' correspond to n! different matrices M (here the factor n! comes).
Let M'' be the same as M' except for zeros on the main diagonal. Then the permanent of M'' is zero. Viewing M'' as an adjacency matrix of a directed graph G, we notice that G cannot have a cycle. Indeed, if there is a cycle x_1 -> x_2 -> ... -> x_k -> x_1, then the set of elements (x_1,x_2), (x_2,x_3), ..., (x_k,x_1) together (y_1,y_1), ..., (y_{n-k},y_{n-k}), where { y_1,...,y_{n-k} } is the complement of { x_1, ..., x_k } in the set { 1, 2, ..., n }, form a set of elements of the matrix M' whose product is 1, making the permanent of M' greater than 1.
This works in the reverse direction as well, resulting in the statement: The permanent of M' is 1 if and only if M'' represents the adjacency matrix of some DAG. Therefore there exist A003024(n) distinct matrices M'. - Vladeta Jovovic, Oct 27 2009
LINKS
FORMULA
a(n) = n! * A003024(n). - Vladeta Jovovic, Oct 26 2009
EXAMPLE
a(2)=6 because there are 6 matrices ((1,0),(0,1)), ((0,1),(1,0)), ((0,1),(1,1)), ((1,0),(1,1,)), ((1,1),(0,1)), ((1,1,),(1,0)) with permanent=1.
MAPLE
b:= proc(n) option remember; `if`(n=0, 1, add((-1)^(k+1)*
binomial(n, k)*2^(k*(n-k))*b(n-k), k=1..n))
end:
a:= n-> n!*b(n):
seq(a(n), n=0..14); # Alois P. Heinz, Jun 27 2023
CROSSREFS
Cf. A088672 number of (0, 1)-matrices with zero permanent, A089479 occurrence counts for permanents of all (0, 1)-matrices, A089480 occurrence counts for permanents of non-singular (0, 1)-matrices.
Sequence in context: A233734 A227044 A188420 * A126679 A232689 A165436
KEYWORD
nonn
AUTHOR
Hugo Pfoertner, Nov 05 2003
EXTENSIONS
More terms from Vladeta Jovovic, Oct 26 2009
a(0)=1 prepended by Alois P. Heinz, Jun 27 2023
STATUS
approved

Lookup | Welcome | Wiki | Register | Music | Plot 2 | Demos | Index | Browse | More | WebCam
Contribute new seq. or comment | Format | Style Sheet | Transforms | Superseeker | Recents
The OEIS Community | Maintained by The OEIS Foundation Inc.

License Agreements, Terms of Use, Privacy Policy. .

Last modified December 1 18:10 EST 2023. Contains 367500 sequences. (Running on oeis4.)