A determinant of Stirling cycle numbers counts unlabeled acyclic single-source automata
David Callan

TL;DR
This paper establishes a novel connection between Stirling cycle number determinants and the enumeration of unlabeled acyclic single-source automata, using combinatorial bijections and involutions.
Contribution
It introduces a new combinatorial interpretation of a determinant of Stirling cycle numbers as counting specific automata structures.
Findings
Determinant counts unlabeled acyclic single-source automata
Bijection between automata and marked lattice paths
Sign-reversing involution used for evaluation
Abstract
We show that a determinant of Stirling cycle numbers counts unlabeled acyclic single-source automata. The proof involves a bijection from these automata to certain marked lattice paths and a sign-reversing involution to evaluate the determinant.
Peer Reviews
No public reviews on file for this paper yet. If you reviewed it on a platform where reviews are public (OpenReview, ICLR, NeurIPS, ICML), you can paste yours below so the community can read it here.
Code & Models
Videos
No videos yet. Explain this paper in a talk, walkthrough, or lecture? Add one.
Taxonomy
Topicssemigroups and automata theory · Chemical Synthesis and Analysis · Molecular Junctions and Nanostructures
A Determinant of Stirling Cycle Numbers Counts Unlabeled Acyclic Single-Source Automata
DAVID CALLAN
Department of Statistics
University of Wisconsin-Madison
1300 University Ave
Madison, WI 53706-1532
March 30, 2007
Abstract
We show that a determinant of Stirling cycle numbers counts unlabeled acyclic single-source automata. The proof involves a bijection from these automata to certain marked lattice paths and a sign-reversing involution to evaluate the determinant.
**1 Introduction ** The chief purpose of this paper is to show bijectively that a determinant of Stirling cycle numbers counts unlabeled acyclic single-source automata. Specifically, let denote the matrix with entry , where is the Stirling cycle number, the number of permutations on with cycles. For example,
[TABLE]
As evident in the example, is formed from copies of each of rows 2 through of the Stirling cycle triangle, arranged so that the first nonzero entry in each row is a 1 and, after the first row, this 1 occurs just before the main diagonal; in other words, is a Hessenberg matrix with 1s on the infra-diagonal. We will show
Main Theorem**.**
The determinant of is the number of unlabeled acyclic single-source automata with transient states on a -letter input alphabet.
Section 2 reviews basic terminology for automata and recurrence relations to count finite acyclic automata. Section 3 introduces column-marked subdiagonal paths, which play an intermediate role, and a way to code them. Section 4 presents a bijection from these column-marked subdiagonal paths to unlabeled acyclic single-source automata. Finally, Section 5 evaluates using a sign-reversing involution and shows that the determinant counts the codes for column-marked subdiagonal paths.
**2 Automata **
A (complete, deterministic) automaton consists of a set of states and an input alphabet whose letters transform the states among themselves: a letter and a state produce another state (possibly the same one). A finite automaton (finite set of states, finite input alphabet of, say, letters) can be represented as a -regular directed multigraph with ordered edges: the vertices represent the states and the first, second, …edge from a vertex give the effect of the first, second, …alphabet letter on that state. A finite automaton cannot be acyclic in the usual sense of no cycles: pick a vertex and follow any path from it. This path must ultimately hit a previously encountered vertex, thereby creating a cycle. So the term acyclic is used in the looser sense that only one vertex, called the sink, is involved in cycles. This means that all edges from the sink loop back to itself (and may safely be omitted) and all other paths feed into the sink.
A non-sink state is called transient. The size of an acyclic automaton is the number of transient states. An acyclic automaton of size thus has transient states which we label and a sink, labeled . Liskovets [1] uses the inclusion-exclusion principle (more about this below) to obtain the following recurrence relation for the number of acyclic automata of size on a -letter input alphabet ():
[TABLE]
A source is a vertex with no incoming edges. A finite acyclic automaton has at least one source because a path traversed backward must have distinct vertices and so cannot continue indefinitely. An automaton is single-source (or initially connected) if it has only one source. Let denote the set of single-source acyclic finite (SAF) automata on a -letter input alphabet with vertices where is the source and is the sink, and set . The two-line representation of an automaton in is the matrix whose columns list the edges in order. For example,
[TABLE]
is in and the source-to-sink paths in include , where the alphabet is .
Proposition 1**.**
The number of SAF automata of size on a -letter input alphabet is given by
[TABLE]
Remark This formula is a bit more succinct than the the recurrence in [1, Theorem 3.2].
Proof Consider the set of acyclic automata with transient vertices in which 1 is a source. Call the interior vertices. For , let
[TABLE]
Then and by Möbius inversion [2] on the lattice of subsets of where is the Möbius function for this lattice. Since if , we have in particular that
[TABLE]
Let so that . When consists entirely of sources, the vertices in and their incident edges form a subautomaton with transient states; there are such. Also, all edges from the vertices comprising go directly into choices. Thus . By definition, is the number of automata in for which 1 is the only source, that is, and the Proposition now follows from (1). ∎
An unlabeled SAF automaton is an equivalence class of SAF automata under relabeling of the interior vertices. Liskovets notes [1] (and we prove below) that has no nontrivial automorphisms, that is, each of the relabelings of the interior vertices of produces a different automaton. So unlabeled SAF automata of size on a -letter alphabet are counted by . The next result establishes a canonical representative in each relabeling class.
Proposition 2**.**
Each equivalence class in under relabeling of interior vertices has size and contains exactly one SAF automaton with the “last occurrences increasing” property: the last occurrences of the interior vertices——in the bottom row of its two-line representation occur in that order.
Proof The first assertion follows from the fact that the interior vertices of an automaton can be distinguished intrinsically, that is, independent of their labeling. To see this, first mark the source, namely 1, with a mark (new label) and observe that there exists at least one interior vertex whose only incoming edge(s) are from the source (the only currently marked vertex) for otherwise a cycle would be present. For each such interior vertex , choose the last edge from the marked vertex to using the built-in ordering of these edges. This determines an order on these vertices; mark them in order . If there still remain unmarked interior vertices, at least one of them has incoming edges only from a marked vertex or again a cycle would be present. For each such vertex, use the last incoming edge from a marked vertex, where now edges are arranged in order of initial vertex with the built-in order breaking ties, to order and mark these vertices . Proceed similarly until all interior vertices are marked. For example, for
[TABLE]
and there is just one interior vertex, namely 4, whose only incoming edge is from the source, and so and 4 becomes a marked vertex. Now all incoming edges to both 3 and 5 are from marked vertices and the last such edges (built-in order comes into play) are and putting vertices 3, 5 in the order 5, 3. So and . Finally, . This proves the first assertion. By construction of the s, relabeling each interior vertex with the subscript of its corresponding produces an automaton in with the “last occurrences increasing” property and is the only relabeling that does so. The example yields
[TABLE]
∎
Now let denote the set of canonical SAF automata in representing unlabeled automata; thus . Henceforth, we identify an unlabeled automaton with its canonical representative.
**3 Column-Marked Subdiagonal Paths **
A subdiagonal -path is a lattice path of steps and , for east and for north, from to that never rise above the line . Let denote the set of such paths.For , it is clear that is nonempty only for and it is known (generalized ballot theorem) that
[TABLE]
A path in can be coded by the heights of its steps above the line ; this gives a a sequence subject to the restrictions and for all .
A column-marked subdiagonal -path is one in which, for each , one of the lattice squares below the th step and above the horizontal line is marked, say with a ‘’. Let \mathsf{C}_{k}^{\textrm{\raisebox{3.01389pt}{*}}}(n,p) denote the set of such marked paths.
*$$*$$*$$*$$*$$*$$*$$*(0,0)(8,4)y=-1$$y=\frac{1}{2}xA path in \mathsf{C}_{2}^{\textrm{\raisebox{2.7125pt}{*}}}(4,3)
A marked path in \mathsf{C}_{k}^{\textrm{\raisebox{3.01389pt}{*}}}(n,n) can be coded by a sequence of pairs \big{(}(a_{i},b_{i})\big{)}_{i=1}^{kn} where is the code for the underlying path and gives the position of the in the th column. The example is coded by .
An explicit sum for |\,\mathsf{C}_{k}^{\textrm{\raisebox{3.01389pt}{*}}}(n,n)\,| is
[TABLE]
because the summand is the number of ways to insert the ‘’s in the underlying path coded by .
It is also possible to obtain a recurrence for |\,\mathsf{C}_{k}^{\textrm{\raisebox{3.01389pt}{}}}(n,p)\,|, and then, using Prop. 1, to show analytically that |\,\mathsf{C}_{k}^{\textrm{\raisebox{3.01389pt}{}}}(n,n)\,|=|\,\mathcal{C}_{k+1}(n)\,|. However, it is much more pleasant to give a bijection and in the next section we will do so. In particular, the number of SAF automata on a 2-letter alphabet is
[TABLE]
sequence \htmladdnormallinkA082161http://www.research.att.com:80/cgi-bin/access.cgi/as/njas/sequences/eisA.cgi?Anum=A082161 in [3].
**4 Bijection from Paths to Automata **
In this section we exhibit a bijection from \mathsf{C}_{k}^{\textrm{\raisebox{3.01389pt}{*}}}(n,n) to . Using the illustrated path as a working example with and ,
*$$*$$*$$*$$*$$*$$*$$*(0,0)(8,4)y=-1$$y=\frac{1}{2}x
first construct the top row of a two-line representation consisting of each 1s, 2s, …, s and number them left to right:
[TABLE]
The last step in the path is necessarily an step. For the second last, third last,… steps in the path, count the number of steps following it. This gives a sequence satisfying and for all . Circle the positions in the two-line representation and then insert (in boldface) in the second row in the circled positions:
[TABLE]
These will be the last occurrences of in the second row. Working from the last column in the path back to the first, fill in the blanks in the second row left to right as follows. Count the number of squares from the up to the path (including the square) and add this number to the nearest boldface number to the left of the current blank entry (if there are no boldface numbers to the left, add this number to 1) and insert the result in the current blank square. In the example the numbers of squares are 2,3,1,2,1,2,1,1 yielding
[TABLE]
This will fill all blank entries except the last. Note that s in the bottom row correspond to sink (that is, ) labels in the second row. Finally, insert into the last remaining blank space to give the image automaton:
[TABLE]
This process is fully reversible and the map is a bijection. ∎
**5 Evaluation of det **
For simplicity, we treat the case , leaving the generalization to arbitrary as a not-too-difficult exercise for the interested reader. Write for . Thus . From the definition of as a sum of signed products, we
show that is the total weight of certain lists of permutations, each list carrying weight . Then a weight-reversing involution cancels all weights and reduces the problem to counting the surviving lists. These surviving lists are essentially the codes for paths in \mathsf{C}_{1}^{\textrm{\raisebox{3.01389pt}{*}}}(n,p), and the Main Theorem follows from §4.
To describe the permutations giving a nonzero contribution to , define the code of a permutation on to be the list with . Since the entry of , is 0 unless , we must have for all . It is well known that there are such permutations, corresponding to compositions of , with codes characterized by the following four conditions: (i) for all , (ii) , (iii) each is immediately followed by zeros in the list, (iv) . Let us call such a list a padded composition of : deleting the zeros is a bijection to ordinary compositions of . For example, is a padded composition of 6. For a permutation with padded composition code c, the nonzero entries in c give the cycle lengths of . Hence sgn , which is the parity of “cycles in ”, is given by .
We have , and so
[TABLE]
where the sum is restricted to padded compositions c of with for all (\htmladdnormallinkA002083http://www.research.att.com:80/cgi-bin/access.cgi/as/njas/sequences/eisA.cgi?Anum=A002083) because unless .
Henceforth, let us write all permutations in standard cycle form whereby the smallest entry occurs first in each cycle and these smallest entries increase left to right. Thus, with dashes separating cycles, 154-2-36 is the standard cycle form of the permutation . We define a nonfirst entry to be one that does not start a cycle. Thus the preceding permutation has 3 nonfirst entries: 5,4,6. Note that the number of nonfirst entries is 0 only for the identity permutation. We denote an identity permutation (of any size) by .
By definition of Stirling cycle number, the product in (2) counts lists of permutations where is a permutation on with cycles, equivalently, with nonfirst entries. So define to be the set all lists of permutations where is a permutation on , # nonfirst entries in is , is the transposition (1,2), each nonidentity permutation is immediately followed by ’s where is the number of nonfirst entries in (so the total number of nonfirst entries is ). Assign a weight to by wt(\pi)=(-1)^{\textrm{\#\,\epsilon\pi}}. Then
[TABLE]
We now define a weight-reversing involution on (most of) . Given , scan the list of its component permutations left to right. Stop at the first one that either (i) has more than one nonfirst entry, or (ii) has only one nonfirst entry, say, and maximum nonfirst entry of the next permutation in the list. Say is the permutation where we stop.
In case (i) decrement (i.e. decrease by 1) the number of ’s in the list by splitting into two nonidentity permutations as follows. Let be the largest nonfirst entry of and let be its predecessor. Replace and its successor in the list (necessarily an ) by the following two permutations: first the transposition and second the permutation obtained from by erasing from its cycle and turning it into a singleton. Here are two examples of this case (recall permutations are in standard cycle form and, for clarity, singleton cycles are not shown).
[TABLE]
and
[TABLE]
The reader may readily check that this sends case (i) to case (ii).
In case (ii), is a transposition with maximum nonfirst entry of . In this case, increment the number of ’s in the list by combining and into a single permutation followed by an : in is a singleton; delete this singleton and insert immediately after in (in the same cycle). The reader may check that this reverses the result in the two examples above and, in general, sends case (ii) to case (i). Since the map alters the number of ’s in the list by 1, it is clearly weight-reversing. The map fails only for lists that both consist entirely of transpositions and have the form
[TABLE]
Such lists have weight 1. Hence is the number of lists \big{(}(a_{i},b_{i})\big{)}_{i=1}^{n} satisfying for , and . After subtracting 1 from each , these lists code the paths in \mathsf{C}_{1}^{\textrm{\raisebox{3.01389pt}{}}}(n,n) and, using §4, \det A(n)=|\,\mathsf{C}_{1}^{\textrm{\raisebox{3.01389pt}{}}}(n,n)\,|=|\,\mathcal{C}_{2}(n)\,|.
The reference list from the paper itself. Each links out to its DOI / PubMed record.
- 1[1] Valery A. Liskovets, Exact enumeration of acyclic deterministic automata, Disc. Appl. Math. , in press, 2006. Earlier version available at \htmladdnormallink http://www.i 3s.unice.fr/fpsac/FPSAC 03/articles.htmlhttp://www.i 3s.unice.fr/fpsac/FPSAC 03/articles.html
- 2[2] J. H. van Lint and R. M. Wilson, A Course in Combinatorics, 2nd ed., Cambridge University Press, NY, 2001.
- 3[3] Neil J. Sloane (founder and maintainer), The On-Line Encyclopedia of Integer Sequences \htmladdnormallink http://www.research.att.com:80/ njas/sequences/index.html?blank=1http://www.research.att.com:80/ njas/sequences/index.html?blank=1
