Abstract
In this note we describe how to apply group actions in order to enumerate the isomorphism classes of rhythmic canons. |
Before describing rhythmic canons in cyclic time, we view them as they may appear within a musical composition, in free linear time which has no cyclicity. Like in a melodic canon, one has several voices that may enter one after the other until all voices are present. As in the case of a melodic canon, all voices are just copies of a ground voice that is suitably translated in the time axis. For simplicity we suppose here that all voices are extended in both directions of the time axis ad infinitum. We further suppose that the ground voice is a periodic rhythm that we will call the inner rhythm. Following Vuza’s definition, a periodic rhythm is an infinite subset R of the rationals (marking the attack times, or onsets) with R = R + d for a suitable period d. Furthermore, R is supposed to be locally finite (i.e. the intersection of R with every time segment [a, b] (of finite length b - a > 0) is finite). The period of a periodic rhythm R is the smallest positive rational number d = d(R) satisfying R = R + d. We also mention another important characteristic of a periodic rhythm, its pulsation p = p(R). It is defined as the max , in other words, it is the biggest rational q such that all distances between the attack points of R are integer multiples of q. Obviously, the period d is an integer multiple of the pulsation p.
Let A be the set of all rational numbers a such that the attack times of the voices of the canon can be expressed as R + a. Then A itself is a periodic rhythm, called the outer rhythm of the canon, with period d(A), where d(R) is an integer multiple of d(A). Note that R and A may have different pulsations p(R) and p(A). Hence, the pulsation of a canon is defined as the greatest rational number p such that both p(R) and p(A) are integer multiples of p.
In order to switch from linear time to circular time the fraction n = d(R)/p must be computed. Let r denote a fixed attack time within the inner rhythm R. Then each attack time in any of the voices is of the form r + tp for a suitable integer t, i.e. the whole canon is contained in r + p . Because everything is periodic with period d(R), we can restrict to the factor space (r + p)/d(R) which may be identified with /n, the residue class ring of modulo n.
From now on, we consider the whole canon within Zn := /n and assume that R, A, and V a denote the projections of the inner rhythm, the outer rhythm, and the voices of a canon.
The concept of a canon used in the present paper is described by G. Mazzola in [4] and was presented by him to the author in the following way: A canon is a subset K Zn together with a covering of K by pairwise different subsets, the voices, V iØ for 1 < i < t, where t > 1 is the number of voices of the canon,
such that for all i, j
We prefer to write a canon K as a set of its subsets V i. Two canons K = and L = are called isomorphic if s = t and if there exists a translation T of Zn and a permutation in the symmetric group St such that T (V i) = W (i) for 1 < i < t. Then obviously T (K) = L.
The aim of this note is to determine the number of non-isomorphic canons for given n. First we present the definition of canons in more details and give a short description of group actions, which are the standard tool to describe isomorphism classes of different objects (cf.[2, 3]). Since we are interested in the attack times modulo n, we assume that K is a subset of Zn. The set of translations of Zn is the cyclic group Cn generated by the permutation n := (0, 1, ..., n - 1) which maps i to i + 1 mod n. Cn acts in a natural way both on Zn and on the set of all subsets of Zn:
As usual we identify a subset A of Zn with its characteristic function A : Zn given by
The set of all functions from a set X to a set Y will be indicated as Y X. For arbitrary f Zn let denote the subset f-1() of Z n, whence f = .
Sometimes it is convenient to write functions f from Zn to as vectors of the form (f(0), f(1), ..., f(n - 1)) using the natural order of the elements of Zn, because then the set of all these functions is totally ordered by the lexicographic order. For functions f, g Zn we say f < g if there exists an i Z n such that f(j) = g(j) for 0 < j < i and f(i) < g(i). The group action of Cn on the set of all subsets of Zn described as an action on the set of all characteristic function is the following
As the canonical representative of the orbit Cn(f) = we choose the function f0 Cn(f) such that f0 < g for all g Cn(f). Moreover, we choose 0 as the canonical representative of the orbit Cn( ) = . In general, representatives of orbits Cn(f) of functions f under the action of the cyclic group Cn are called necklaces.
A function f Zn (or the corresponding vector and the set ) is called acyclic if C n(f) consists of n different objects. The canonical representative of the orbit of an acyclic function is usually called a Lyndon word. If f is acyclic, then all elements in the orbit Cn(f) are acyclic as well, and we call Cn(f) an acyclic orbit.
From the first two properties of a canon we derive that all voices V i belong to the same orbit Cn(V 1) and that each voice is acyclic. For that reason, we can describe a canon K = as a pair (L, A) where L is the Lyndon word, which is the canonical representative of the orbit Cn(V 1), and A = is a t-subset of Zn such that V i = nai(L). In other words, L describes the inner rhythm and A the outer rhythm of K.
Now we analyse in which situations a pair (L, A), as described above, is not a canon, i.e. when it fails to fulfill the third property. Assume that (L, A) is not a canon, then the differences K - K generate a proper subgroup of Zn. This subgroup is isomorphic to Zn/d where d > 1 is a divisor of n. In other words, d is a divisor of all differences k - l for all k, l K. We write d | r in order to express that the integer d is a divisor of the integer r. (It is not connected with division in the ring Zn.) It is easy to prove the following
Consider the set S of all pairs (L, A) where L is a Lyndon word and A is a t-subset of Zn. The group of translations of Zn acts on S according to
If (L, A) is a canon, then the orbit Cn(L, A) = (L, Cn(A)) describes the isomorphism class of (L, A). In general, the canonical representative of Cn(L, A) is (L, A0) where A0 is the canonical representative of Cn(A). From Lemma 1 we already know that A0(n - 1) = 1 and L(n - 1) = 1 if L0. This together with Lemma 2 proves
In the next part we show that functions f Zn with the property f(i) = 1 implies i -1 mod d (for d | n, d > 1) can be constructed from functions in Zn/d . Let d be a function defined on , such that d(0) is the vector (0, 0, ..., 0) consisting of d entries of 0, and d(1) = (0, ..., 0, 1) is a vector consisting of d - 1 entries of 0 and 1 in the last position. We write the values of d in the form
If we apply d to each component of a vector f Zr by replacing each component 0 in f by 0d and each component 1 by 0d-11, we get a vector d(f) Zrd and
From Lemma 3 we conclude that (L, A) does not describe a canon if and only if there exists a d > 1 which divides n such that both L and A0 are elements of d( Zn/d ).
Some properties of d are collected in the next
In order to show that the second item is true, assume that d(f o r)(j) = 1 for some j Zrd. According to the definition of d, this is equivalent to j = sd - 1 and (f o r)(s - 1) = 1 for some s . In other words, j = sd - 1 and f(s) = 1, which is equivalent to j = sd - 1 and d(f)((s + 1)d - 1) = 1. This, however, is the same as d(f)(j + d) = 1, whence (d(f) o rdd)(j) = 1. If f = 0 the assertion is always true.
The third part is an immediate consequence of the second.
If f < g, then there exists i Zr such that f(j) = g(j) for all j < i and f(i) < g(i), whence f(i) = 0 and g(i) = 1. For that reason, d(f(i)) = 0d and d(g(i)) = 0d-11. Hence, d(f)(j) = d(g)(j) for j < id-1 and d(f)(id-1) = 0 < 1 = d(g)(id-1), and consequently d(f) < d(g). If, conversely, d(f) < d(g), then there exists an i Zrd such that d(f)(j) = d(g)(j) for all j < i and d(f)(i) < d(g)(i), whence d(g)(i) = 1 and i -1 mod d. Assume that i is of the form sd - 1. Then d(f(j)) = d(g(j)) for j < s, and d(f(s)) = 0d whereas d(g(s)) = 0d-11. Since d is an injective mapping, f(j) = g(j) for j < s and f(s) = 0 < 1 = g(s), which implies that f < g.
From the definition of the canonical representative of an orbit and from the items 4. and 2. it follows immediately that d(f0) < d(f0 o rj) = d(f0) o rddj for all j . Moreover, if n / 0 mod d, then (d(f0) o rdn)(rd - 1) = d(f0)(n - 1)1, since n - 1 / -1 mod d. According to Lemma 2, the representative d(f0) o rdn cannot be the canonical representative of the orbit Crd(d(f0)). Hence, the canonical representative is d(f0).
Then 6. is a trivial consequence of 5. Similar arguments can be applied for proving 7. and 9. Item 8. follows immediately from 7.
In order to compute the number of all Cr-orbits on Zr we apply Pólya’s theory cf. [2, 3, 5, 6]. In the present situation we have to compute the cycle index of the group Cr acting on Zr which is a polynomial in z1, ..., zr over given by
where is the Euler totient function. Replacing each indeterminate zi in C(Cr, Zr) by , which is equal to 2, we compute the number of all Cr-orbits on Zr as
If we replace zi by 1 + zi, where z is an indeterminate over , then the number of C r-orbits of k-sets A Zr is the coefficient of zk in
In conclusion (for the definition of (r) see 6. in Lemma 4)
since we don’t count the orbit of the empty set. Similar methods can be applied to enumerate the number of acyclic Cr-orbits on Zr or, what is equivalent to this, to enumerate all Lyndon words of length r over the alphabet . For r > 1, this number is given as (for the definition of (r) see 8. in Lemma 4)
where is the classical Moebius function. Moreover, the number of Lyndon words of length r over having r1 components 1 and r - r1 components 0 is
And the generating function of Lyndon words of length r over the alphabet with r1 components 1 is of the form
In the case r = 1 there are two Lyndon words, namely f0 = (0) and f1 = (1). The first one is the characteristic function of the empty set in Z1 so we don’t want to count it. For that reason, we must set (1) = 1. Moreover, it should be mentioned that f0 is the unique Lyndon word such that d(f0) is not a Lyndon word for d > 1.
Finally, we have to combine all these results for enumerating the isomorphism classes of canons. Let n > 1 be an integer. For any divisor d > 1 of n, let Mn,d be the set of all pairs (L, Cn(A)), where L is a Lyndon word of length n over , (in the case n = 1 different from 0) such that L(i) = 1 implies i -1 mod d, and A is a non-empty subset of Zn, such that d | k - l for all k, l A. Hence
From Lemma 4 we deduce that d describes a bijection between Mn,d and Mn/d,1, thus
which is the number of possibilities to choose L and Cn(A) according to the desired properties of Mn,d. Finally, let n be the set of isomorphism classes of canons in Zn,
From Lemma 3 we deduce that
where
and
It is clear that Mn,1 contains this union. Moreover, this union is disjoint, since for a canon K = (L, Cn/d(A)) n/d we have
Finally, we have to show that each element of Mn,1 belongs to this union. If (L, Cn(A)) Mn,1, then choose the biggest d such that (L, Cn(A)) belongs to Mn,d. (I.e. (L, Cn(A)) does not belong to Mn,dd' for d' > 1. It is always possible to find such d, since, if (L, Cn(A)) belongs both to Mn,d1 and Mn,d2, then it also belongs to Mn,lcm(d1,d2).) Then there exists (L', C n/d(A')) M n/d,1 such that L = d(L') and A0 = d(A0' ) for the canonical representatives A0 and A0' of C n(A) and Cn/d(A'). Moreover, (L', C n/d(A')) n/d because otherwise there would be some d' > 1 which is a divisor of n/d such that (L', C n/d(A')) M n/d,d' . But then (L, Cn(A)) also belongs to Mn,dd' , which is a contradiction to the choice of d.
Hence,
and by Moebius inversion we get
which finishes the proof.
If and are replaced by the corresponding generating functions
and
then the coefficient of yszt in
gives the number of non-isomorphic canons in Zn which consist of t voices where each voice contains exactly s attack times.
For example the numbers of isomorphism classes of canons in Zn for 1 < n < 10 are:
Finally, we have a closer look at the 13 isomorphism classes of canons in Z4, which can be classified according to their number of voices and attack times per voice by z4y3 + z4y2 + z4y + z3y3 + z3y2 + z3y + 2z2y3 + 2z2y2 + z2y + zy3 + zy2. From this list we can read, for instance, that there are exactly three classes of canons with 4 voices which have 1, 2, or 3 attack times per voice. Or we see that there are five classes of canons with 2 voices, two of them having 3 attack times per voice, two of them having 2 attack times per voice, and only one having 1 attack time in each of its two voices.
We can even compute the canonical representative of each isomorphism class. For instance, there are exactly three Lyndon words L1, L2, L3 of length 4 over and five representatives f1, ..., f5 of the C4-orbits of non empty subsets of Z4:
The Lyndon words describe the distribution of the attack times in one voice (i.e. the inner rhythm of the canon), the necklaces describe the distribution of the voices over the complete period, here in this example over Z4, which is the outer rhythm. Only L1 can be constructed as 4((1)) or 2((0, 1)) from shorter Lyndon words, so all pairs (Li, fj) for i and j describe canons in Z4. Since f1 = 4((1)) = 2((0, 1)) and f3 = 2((1, 1)), and (1), (0, 1), and (1, 1) are necklaces of length 1 or 2 over , the pairs (L1, C4(f1)) = 4((1), C1(1)) and (L1, C4(f3)) = 2((0, 1), C2(1, 1)) do not belong to 4. But the pairs (L1, fj) for j describe again canons in Z4.
From this example we also see how to compute complete lists of isomorphism classes of canons in Zn for arbitrary n. There exist fast algorithms for computing all Lyndon words and all necklaces of length n over . Then each pair (Li, fj) must be tested whether there exists some d > 1 such that (Li, fj) = (d(L'), d(f')) for a Lyndon word L' and a necklace f' of length n/d.
There exist more complicated definitions of canons. A pair (R, A) of inner and outer rhythm defines a rhythmic tiling canon with voices V a for a A in Zn if and only if
It is obvious that periods and pulsations can be determined in cyclic time as well. Rhythmic tiling canons with the additional property,
are called regular complementary canons of maximal category. So far we were not dealing with that kind of canons.
Acknowledgment: The author wants to express his thanks to Guerino Mazzola for pointing his attention to the problem of classification of canons. He is also grateful to Moreno Andreatta, who provided very useful historic and background information about different definitions of canons (cf [1]).
[1] M. Andreatta, Th. Noll, C. Agon, and G. Assayag. The geometrical groove: rhythmic canons between theory, implementation and musical experiment. In Les Actes des 8e Journées d’Informatique Musicale, Bourges 7-9 juin 2001, pages 93-97, 2001.
[2] A. Kerber. Algebraic Combinatorics via Finite Group Actions. B.I. Wissenschaftsverlag, Mannheim, Wien, Zürich, 1991. ISBN 3-411-14521-8.
[3] A. Kerber. Applied Finite Group Actions, volume 19 of Algorithms and Combinatorics. Springer, Berlin, Heidelberg, New York, 1999. ISBN 3-540-65941-2.
[4] G. Mazzola. Topos of music. To appear.
[5] G. Pólya. Kombinatorische Anzahlbestimmungen für Gruppen, Graphen und chemische Verbindungen. Acta Mathematica, 68:145 - 254, 1937.
[6] G. Pólya and R.C. Read. Combinatorial Enumeration of Groups, Graphs and Chemical Compounds. Springer Verlag, New York, Berlin, Heidelberg, London, Paris, Tokyo, 1987. ISBN 0-387-96413-4 or ISBN 3-540-96413-4.
[7] D.T. Vuza. Supplementary sets and regular complementary unending canons (part one). Perspectives of New Music, 29(2):22-49, 1991.
[8] D.T. Vuza. Supplementary sets and regular complementary unending canons (part three). Perspectives of New Music, 30(2):102-125, 1992.
[9] D.T. Vuza. Supplementary sets and regular complementary unending canons (part two). Perspectives of New Music, 30(1):184-207, 1992.
[10] D.T. Vuza. Supplementary sets and regular complementary unending canons (part four). Perspectives of New Music, 31(1):270-305, 1993.
Harald Fripertinger
Institut für Mathematik Karl Franzens Universität Graz Heinrichstr. 36/4 A-8010 Graz, AUSTRIA harald.fripertinger@uni-graz.at |