arxiv: v2 [cs.dm] 18 Nov 2010

Transcription

1 Dynamics in parallel of double Boolean automata circuits Mathilde Noual, arxiv:0.90v [cs.dm] 8 Nov 00 August, 08 Université de Lyon, ÉNS-Lyon, LIP, CNRS UMR8, 900 Lyon, France IXXI, Institut rhône-alpin des systèmes complexes, 900 Lyon, France Introduction In this paper, we give some results concerning the dynamics of double Boolean automata circuits (dbac s for short), namely, networks associated to interaction graphs composed of two side-circuits that share a node. More precisely, a double circuit of left-size l N and of right-size r N is a graph that we denote by D l,r. It has n = l r nodes. Nodes that are numbered from 0 to l belong to the left-circuit and the others plus the node 0 (that belongs to both side-circuits) belong to the right-circuit. Node 0 is the only node with in- and out-degree. All other nodes have in- and out-degree. C l 0 l C r l n = l r Figure : Double circuit D l,r.

2 A dbac D l,r = (D l,r, F ) is a network defined by (i) its interaction graph, a double-circuit D l,r, and (ii) a global transition function F that updates the states of all nodes in parallel and that is defined as follows by the local transition functions f i of nodes i < n: x {0, } n, F (x) i = f i (x i ), i / {0, l}, F (x) l = f l (x 0 ), F (x) 0 = f 0 (x l, x lr ) = f0 L(x l ) f0 R(x lr ) where {, }. () All local transition functions are supposed to be non-constant. Thus, i < n, f i, f0 L, f0 R {id, neg} where id : a a and neg : a a, a {0, }. As a consequence all local transition functions are locally monotone. All arcs (i, j) entering a node j whose local transition function is id (resp. neg), with respect to node i, are labeled with a sign (resp. a sign) and called positive arcs (resp. negative arcs). A sidecircuit with an even number of negative arcs (resp. odd number of negative arcs) is called a positive (side-) circuit (resp. a negative (side-) circuit). Given a configuration x = (x 0,..., x n ) {0, } n of a dbac D l,r, we use the following notation: x L = (x 0,..., x l ) and x R = (x 0, x l,..., x n ). A configuration x(t) = (x 0 (t),..., x n (t)) {0, } n such that k N, F k p (x(t)) = x(tk p) = x(t) is said to have period p. If x(t) has period p and does not also have period d < p, then x(t) is said to have exact period p. An attractor of period p N, or p-attractor, is the set of configurations belonging to the orbit of a configuration that has p as exact period. Attractors of period are called fixed points. The graph whose nodes are the configurations x {0, } n of a network and whose arcs represent the transitions (x(t), x(t ) = F (x(t)) is called the transition graph of the network. In [], the authors showed the following results: Proposition. The transition graphs of two dbacs with same side-signs and side-sizes are isomorphic, whatever the definition of f 0 (i.e., whether = or = in the definition () of F above).. Attractor periods of a dbac divide the sizes of the positive side-circuits if there are some and do not divide the sizes of the negative side-circuits if there are some.

3 . If both side-circuits of a dbac D l,r have the same sign, then, attractor periods divide the sum N = l r.. If both side-circuits of a dbac D l,l have the same sign and size, then D l,l behaves as an isolated circuit of that size and sign (i.e., the sub-transition graph generated by the periodic configurations of D l,l is isomorphic to the transition graph of an isolated circuit with the same sign and size).. A dbac has as many fixed points as it has positive side-circuits.. If both side-circuits of a dbac D l,r are positive and p divides l and r, then number of attractors of period p is given by A p (sequence A0 of the OEIS []), namely, the number of attractors of period p of an isolated positive circuit of size a multiple of p. As a result of the first two points of Proposition, we may focus on canonical instances of dbacs. Thus, from now on, we will suppose that = and i 0, f i = id. If the left-circuit is positive (resp. negative), we will suppose that the arc (l, 0) is positive (resp. negative) and f0 L = id (resp. f 0 L = neg) and similarly for the right-circuit. Thus, the only possible negative arcs on the dbacs we will study will be the arcs (l, 0) and (n, 0). The last point of Proposition yields a description of the dynamics of a doubly positive dbac (in terms of combinatorics only but [] also gives a characterisation the configurations of period p, p N for all types of dbacs). Thus, here, we will focus on the cases where the dbacs have at least one negative side-circuit. For a negative-positive or a negative-negative dbac D l,r, we will use the following notations and results. The function µ is the Mobiüs function. It appears in the expressions below because of the Mobiüs inversion formula []. The number of configurations of period p is written: C p (l, r). The number of configurations of period exactly p is written and given by: C p(l, r) = q p µ( p q ) C p(l, r).

4 The number of attractors of period p (where p is an attractor period) is written and given by: A p (l, r) = C p(l, r) p = p µ( p q ) C p(l, r). q p In particular, C (l, r) = C (l, r) =. The total number of attractors is written and given by: T(l, r) = A p (l, r). p attractor period For the sake of clarity, we will write A p (l, r) = A ± p (l, r) (resp. A p (l, r) = A = p (l, r)) and T(l, r) = T ± (l, r) (resp. T(l, r) = T = (l, r)) when D l,r will be a negative-positive dbac (resp. a negative-negative dbac). Positive-Negative We first concentrate on dbacs whose left-circuit is negative and whose rightcircuit is positive. From Proposition, we know that all possible attractor periods l 0 l n = l r l r Figure : Interaction graph of a negative-positive dbac. All arcs are positive (resp. all local transition functions are equal to id) except for the arc (l, 0) (resp. except for the local transition function f0 L = neg). of such dbacs divide r and do not divide l. We also know that these networks have exactly one fixed point. In the sequel, we focus on attractor periods p >.

5 positive l r T l negative gcd(l, r) = T ± l,r T r r r T r Table : Total number of attractors of a negative-positive dbac D l,r (obtained by computer simulations). Each colour corresponds to a value of gcd(l, r). The line T r (resp. the column T l ) gives the total number of attractors of an isolated positive (resp. negative) circuit of size r (resp. l).

6 Characterisation of configurations of period p Let p be a divisor of r = p q that does not divide l = k p d, d = l mod p. We write = gcd(p, l) = gcd(p, d). Let x = x(t) be a configuration of period p: m N, x(tm p) = x(t). Note that to describe the dynamics of D l,r, it suffices to describe the behaviour of node 0 (see []). The configuration x satisfies the following: x 0 (t) = x 0 (t k p) = x l (t k p ) x n (t k p ) = x l k p (t) x 0 (t (k q) p) = x d (t) x 0 (t) Thus, if x 0 (t) = 0, then x d (t) = (and also, because x 0 (t) = x 0 (t d) x 0 (t), if x 0 (t) = 0, then x 0 (t d) = ). From this, we may derive the following characterisation: Proposition Let p N be divisor of r = p q that does not divide l. A configuration x = x(t) has period p if and only if there exists a circular word w {0, } p of size p that does not contain the sub-sequence 0u0, u {0, } d and that satisfies: x L = w k w[0... d ] and x R = w q where w[0... m] = w 0... w m. More precisely, the word w satisfies: i < p, w i = x 0 (t i). Consequently, the property for a configuration to be of period p depends only on d = l mod p and on p (and not on l nor on r). As a result of the characterisation above, we may focus on negative-positive dbacs D l,r such that l < r ( i.e., l = l mod r) and to count the number of attractors of period p we may focus on dbacs D l,r such that l < p (i.e., l = d = l mod p). In other words, from Proposition : Combinatorics A ± p (l, r) = A ± p (l mod r, r) = A ± p (l mod p, r). From Proposition each configuration x of period p is associated with a circular word w {0, } p that does not contain the sub-sequence 0u0, u {0, } d, d = l mod p. It is easy to check that this word w can be written as an interlock of a certain number N of circular words w (), w (),..., w (N) of size T = p/n that do not contain the sub-sequence 00 (see figure ). More precisely, the size of a word w (j) satisfies T d = K p for a certain K such that T and K are minimal. In other words, d p gcd(d,p) p T d = lcm(d, p) = and thus T =. Consequently, N = and we obtain the following lemma which explains why, in each column of Table, every cell of a same colour contains the same number.

7 w () w () w () w w () w w w 0 w () w () 0 w w 0 w w () 0 w w w () 0 w w w () w () w () w () w 9 w () w 8 w w () w w () Figure : The circular word w = w 0... w p = x 0 (t)... x 0 (t p ) mentioned in Proposition that characterises a configuration x(t) of period w = p. In this example, p =, d = l mod p = so that w is made of an interlock of = gcd(d, p) = words w (j) = w (j) 0... w (j) of size p/ =. Lemma Let p be a divisor of r = p q that does not divide l. The number of configurations of period p, C p (l, r), depends only on = gcd(l, p) and on p. Thus, we write: C p (l, r) = C p, p. The number of circular words of size n that do not contain the sub-sequence 00 is counted by the Lucas sequence (sequence A0 of the OEIS []): L() = L() = L(n) = L(n ) L(n ) = φ n φ n = φ n ( φ )n, where φ =.8099 is the golden ratio and φ = φ = Among the properties of φ that will be useful to us in the sequel are the following: φ = φ and φ = φ. Thus, to build a circular word w {0, } p without the sub-sequence 0u0, u {0, } l, one needs to chose among L( p ) words w (j) of size p without the sub-sequence 00. As a result holds Proposition below:

8 Proposition The number of configurations of period p is given by: C p, p = L( p ) p. Consequently, the number of attractors of period p is given by: where = gcd(l, p). A ± p (l, r) = A p, p = p µ( p q ) L( p ) p, Number of configurations of period p Let us develop the expression for C p (l, r) = C p, p : If p is odd, p. p and p C p, p = L( p ) p = L(φ p p q p ( φ) p p ) p = ( p ) p k p k p k k φ p ( φ) p = ( φ) p k ( p ) p k φ p k p ( ) p k p = φ p k ( p ) p k ( φ ) p k p = φ p (( φ ) p p ) p = ( ) p φ p (( ) p p (φ ) p p ) p. cannot be even. Thus, there are three cases only: are odd. Thus, because is necessarily also odd: C p, p = φ p ( (φ ) p p ) p = φ p ((φ ) p p ) p. (). p is even and p is odd. Thus, because is necessarily even: C p, p = φ p ( (φ ) p p ) p = φ p ((φ ) p p ) p. (). p and p are both even. Thus: C p, p = φ p ((φ ) p p ) p. () 8

9 To sum up, we give below Proposition whose last part can be derived from the relation between the Euler totient ψ( ) and the Mobiüs function µ( ), ψ(n) = m n (n/m) µ(m), and from the following equations where q = gcd(q, l): T ± (l, r) = p r A p (l, r) = p r p µ(p q ) C q, q = r C q, q r (p/q) q µ(p q ) q p = r C q, q q r k r q p r q p r k q µ(k) = r ψ( r q ) C q, q. q r Proposition Let = gcd(l, p) where p is a divisor of r that does not divide l. Then, the number of configurations of period p is given by: φ p ((φ ) p p ) p p if is odd, C p, p = φ p ((φ ) p p ) p if is even. In particular, C r,gcd(l,r) counts the total number of periodic configurations of the network. The number of attractors of period p and the total number of attractors are respectively given by: A ± p (l, r) = A ± p, = µ( p p q ) C q, q, and T ± (l, r) = r q p p p r, (p l) ψ( r p ) C p,. Upper bounds From equations (), () and (), one can derive that C p, p and thus A ± p (l, r) = A ± p, are maximal when is minimal (i.e., = ), if p is odd and if p is even, on the contrary, C p, p and A ± p (l, r) are maximal when is maximal (i.e., = p ). Thus, we have: C p, p φ p ((φ ) p p ) p φ p ((φ ) ) p = C p, p = ( φ) p φ p = ( φ φ ) p = p, 9

10 In addition, A ± p, p F p( ) where F p (a) = d p µ( p d ) ad. Now, it can be shown that: a > φ, p, a p < F p (a) = µ( p d ) ad < a p. d p Consequently, p, A p, p p F p( ) < ( )p p F p() = ( )p A p. () where A p refers to the number of attractors of period p of an isolated positive circuit whose size is a multiple of p []. Then, from (), we may derive the following result: Proposition For all p, the number of attractors of period p satisfies: A ± p (l, r) = A p, p < ( )p A p and for p = : A, = = A. Let us denote by T n the number of attractors of a positive circuit of size n. From Proposition, for all r, the total number of attractors of a negative-positive dbac satisfies: T ± l,r < ( )r T r. However, computer simulations (cf Table, last line) show that this bound on T ± l,r is too large. We leave open the problem of finding a better bound. Negative-Negative We now concentrate on doubly negative dbacs. The canonical dbac we will use in the discussion below is defined in Figure. Let p N be a possible attractor period of D l,r (p divides N = l r but divides neither l nor r). Without loss of generality, suppose l mod p > r mod p = d. Because p divides l r, it holds that l mod p = p d. Then, for any configuration x = x(t) {0, } n of period p, we have the following: x 0 (t) = x l (t ) x n (t ) = x 0 (t l) x 0 (t r) = x 0 (t r) x 0 (t l) = x 0 (t d) x 0 (t d) As a consequence, if x 0 (t) = 0, then x 0 (t d) = x 0 (t d) = and if x 0 (t) =, then either x 0 (t d) = 0, or x 0 (t d) = 0. Thus, the circular word w = 0

11 0 l l n = l r l r Figure : Interaction graph of a negative-negative dbac. All arcs are positive (resp. all local transition functions are equal to id) except for the arcs (l, 0) and (n, 0) (resp. except for the local transition functions f L 0 = neg and f R 0 = neg). x 0 (t)... x 0 (t p ) contains neither the sub-sequence 0u0 nor the sub-sequence uu (u, u {0, } d ). Let = gcd(l, r). As in the previous section, w can be written as an interlock of = gcd(d, p) = gcd(, p) words w (j) of size p/ that do not contain the subsequences 00 and. As one may show by induction, the number of such words is counted by the Perrin sequence [], sequence A8 of the OEIS []: P (0) =, P () = 0, P () =, P (n) = P (n ) P (n ) = α n β n β n, where α, β = ( α i α ) and β are the three roots of x x = 0, and α, the only real root of this equation, is called the plastic number []. Using similar arguments to those used in the previous section, we derive Proposition below. This proposition explains why, in Table, all cells of a same diagonal (i.e., when N = l r is kept constant) that have the same colour also contain the same number: the number of attractors depends only on N = l r and on = gcd(l, r) and not on l nor r. Equations in Proposition exploit, in particular, the fact that if p divides l or r then C p, p = 0 (because then P ( p ) = P () = 0). Proposition Let N = l r and let p N be a possible attractor period of D l,r (p divides N but divides neither l nor r). Let also = gcd(l, r) and = gcd(, p).

12 Then, the number of configurations of period p of the doubly negative dbac D l,r depends only on p and. It is given by: C p (l, r) = C p, p = P ( p ) p. The number of p-attractors and the total number of attractors of a doubly negative dbac D l,r are respectively given by: A = p (l, r) = A = p, = p µ( p q ) P ( q ) q, q q p T = (l, r) = T = N, = N ψ( N p ) P ( p ) p. p N The expression for T = N, in Proposition above simplifies into the following if K = N is a prime: T = N, = N q, gcd(q,k)= ψ(q) P (K) q. In particular, if K = N = or K = N =, then because P () = and P () =, the following holds: T = ( N, N )= T= N, N = N q N, gcd(q,)= ψ(q) N q, T = ( N, N )= T= N, N = N q N, gcd(q,)= ψ(q) N q. From the computer simulations we performed (see Tables and ), we observe the following:. Given l, T = (l, r) is maximal when r = l.. Given an integer N which is not a multiple of, T = N, is maximal when is maximal (in particular, if N is even without being a multiple of, then T = N, T= ). N, N. Given an integer N which is a multiple of, T = N, is maximal when = N. We leave the proofs of these three points as an open problem.

13 negative negative l r l r prime l r even l r odd l = r 9 8 l = r Table : Total number of attractors of a negative-negative dbac D l,r (obtained by computer simulations).

14 negative negative l r l = r l = r 9 8 gcd(l, r) = 8 9 l = r T l Table : Total number of attractors of a negative-negative dbac D l,r (obtained by computer simulations). Each colour corresponds to a value of gcd(l, r). The last column gives the total number T l of attractors of an isolated negative circuit.

15 References [] J. Demongeot, M. Noual and S. Sené, Combinatorics of Boolean automata circuits dynamics, submitted to Discrete applied mathematics, 00. [] T.M. Apostol, Introduction to analytic number theory, Springer-Verlag, 9. [] N.J.A. Sloane, 008. The On-Line Encyclopedia of Integer Sequences. [] W.W. Adams and D. Shanks, Strong primality tests that are not sufficient, Math. Comp. 9, 00, 98. [] A.G. Shannona, P.G. Andersonb and A.F. Horadamc, Properties of Cordonnier, Perrin and Van der Laan numbers, International Journal of Mathematical Education in Science and Technology,, 8 8, 00. [] J. Demongeot, M. Noual and S. Sené, On the number of attractors of positive and negative Boolean automata circuits, WAINA 0 proceedings, IEEE Press, 8 89, 00.