Permanents and determinants of (0,1) matrices
Jaap Spies
j.spies at hccnet.nl
Mon Nov 3 17:32:23 CET 2003
Gordon Royle wrote:
>>n=6: Computing the determinants of all 2^36 6X6 (0,1) matrices took
>>~6 days CPU on a 500 MHz Digital Alphastation (using Gauss elimination).
>>I'll not try the permanents ... maybe someone else in this group with
>>more powerful resources could attack this problem?
>
>
>
> Actually, I think that large computer power is not needed..
>
> Any nxn 0/1 matrix is the bipartite adjacency matrix of a bipartite graph
> with two parts of size n.
>
> Each different labelling of this graph that preserves the two parts (ie does
> not interchange them) will give a different 0/1 matrix with the same
> permanent as the original matrix.
>
> However the number of different labellings of such a graph is
>
> n! x n! / g
>
> where g is the size of the automorphism group of the *bicoloured* graph (ie.
> only count automorphisms that do not exchange the two parts).
>
> The advantage of this is that one graph with a trivial group will account
> for over half a million matrices in one hit, and of course we expect most of
> the graphs to have trivial group....
>
> So, calculate the bicoloured graphs on 6+6 vertices - this takes about 1
> second to compute using Brendan McKay's "genbg" program, as there are only
> 251610 bicoloured graphs.
>
> Find their groups (not permitting the parts to exchange): about 3 seconds to
> compute.
>
> Find their permanents: about 20 seconds using a 10-line naive C program that
> simply goes through all possible permutations to find the permanent.
>
> Each graph with group of size g will yield
>
> 6! x 6! / g
>
> matrices, all with the same permanent.
>
> Add it all up and in less than 30 seconds cpu time we get the values
> below......
>
> Permanents of 6x6 0/1 matrices
> (both singular and non-singular)
> --------------------------------
[...]
This is impressive, but unfortunately only part of the story.
The production of a list like Hugo's needs more power, next to elegance,
I suppose.
Jaap
More information about the SeqFan
mailing list