 # Sets of Stochastic Matrices with Converging Products: Bounds and Complexity

An SIA matrix is a stochastic matrix whose sequence of powers converges to a rank-one matrix. This convergence is desirable in various applications making use of stochastic matrices, such as consensus, distributed optimization and Markov chains. We study the shortest SIA products of sets of matrices. We observe that the shortest SIA product of a set of matrices is usually very short and we provide a first upper bound on the length of the shortest SIA product (if one exists) of any set of stochastic matrices. We also provide an algorithm that decides the existence of an SIA product. When particularized to automata, the problem becomes that of finding periodic synchronizing words, and we develop the consequences of our results in relation with the celebrated Cerny conjecture in automata theory. We also investigate links with the related notions of positive-column, Sarymsakov, and scrambling matrices.

## Authors

##### This week in AI

Get the week's most popular data science and artificial intelligence research sent straight to your inbox every Saturday.

## 1 Introduction

A stochastic matrix is an SIA matrix (Stochastic111A matrix is stochastic if it is nonnegative and the sum of the elements on each row is 1. Indecomposable Aperiodic) if the limit

 limt→+∞Pt

exists and all of its rows are equal, i.e., the limit is a rank-one matrix. These matrices play a fundamental role in the theory of Markov chains as they correspond to chains converging to unique limiting distribution . They also arise in discrete-time consensus systems, which are the systems representing groups of agents trying to agree on some value by iterative averaging . Such systems can be modelled by the update equation

 x(t+1)=P(t)x(t), (1)

where

is the vector of the values of the agents at time

and is a stochastic transition matrix representing how agents compute their new values.

Markov chains, and consensus systems, are often time-inhomogeneous, i.e., the transition matrix can be different at different steps , where typically belongs to a finite set of stochastic matrices . Such chains arise in diverse settings including flocking , simulated annealing in optimization , Monte Carlo methods  and queueing theory .

The analysis of convergence is much more difficult in the inhomogenous case than in the time-invariant case. One of the tools to analyze Markov chains with multiple modes is the classical result of J. Wolfowitz stating that any product of transition matrices from converges to a rank-one matrix if and only if every product of matrices in is SIA . Sets satisfying this condition are called quasidefinite , consensus , or weakly-ergodic matrix sets . Informally, the chain represented by such forgets its distant past. This theorem has found many applications [20, 6] and guarantees convergence of the modelled stochastic system in every possible scenario or switching between matrices in . In the present paper we study sets of matrices that have at least one converging sequence. While the former condition is likely to appear in a general or adversarial setting, the latter is more relevant in a controlled environment.

We will say that a set of stochastic matrices is SIA if there is a product of matrices from with repetitions allowed that is an SIA matrix. It is not hard to show (see Section 2) that has an infinite product converging to a rank-one matrix if and only if is an SIA set. Therefore, SIA sets correspond to stochastic systems that can be brought to convergence by controlling the switching signal. Another motivation to study SIA sets comes from the study of inhomogeneous Markov chains or consensus systems with random switching 

. If at every step the matrix is chosen independently and with non-zero probability from

, then the system converges with probability one if222Due to the fact that any finite product of matrices from appears infinitely many times with probability one and only if is an SIA set.

For an SIA set , we denote by the length of the shortest SIA product of matrices from , and we call this quantity the SIA-index of . The SIA-index is a natural parameter of SIA sets. It is similar and related to the classical and well studied notion of exponent – the length of the shortest product that is entrywise positive, if one exists (see [7, Section 3.5] for a survey of the single matrix case, and [18, 27] for more recent work on matrix sets). Similar quantities have been defined for different matrix classes. For instance, the scrambling-index is defined as the length of the shortest scrambling333These classes will be defined in Section 2. product of a matrix set, if one exists. In this article, we will define indices for all the considered classes: the classes of Sarymsakov, SIA and positive-column matrices33footnotemark: 3.

In the context of a system whose switching sequence can be controlled, an SIA product corresponds to a switching sequence that can be repeated to make the system converge to a rank-one matrix. Thus, shorter SIA products correspond to simpler controllers. Furthermore, the length of the SIA product has an influence on the converging rate. Indeed, if is an SIA product, the sequence converges to a rank-one matrix at an average rate of , where

is the second largest eigenvalue of

.

Finally, the SIA-index brings new insights to synchronizing automata and the Černý conjecture  stating that for any set of stochastic444Typically, the conjecture is posed for deterministic finite state automata, i.e., matrices having exactly one 1 on each row, but the presented version is equivalent . matrices either there is a product of length at most of matrices from having a positive column or there is no such product at all. The conjecture has been open for half a century and the best bound obtained so far is cubic in  . As we will soon see, a good upper bound on the SIA-index of matrices would improve the state of the art on the Černý conjecture. In particular, any subquadratic bound would bring a breakthrough.

Our contribution. We introduce and study the SIA-index of sets of stochastic matrices. We show that the largest value of the SIA-index among all SIA sets of stochastic matrices, , is . Moreover, we show that grows at least as . We conjecture that the actual growth rate is closer to the provided lower bound and support this statement by performing an exhaustive search for small values of on a computer grid.

We show that a SIA set has a scrambling, a Sarymsakov and a positive-column product, and, conversely, a set that has a scrambling, a Sarymsakov or a positive-column product is a SIA set. As a consequence, the same polynomial-time procedure can be used to decide whether a given set of stochastic matrices has any of the aforementioned products. We show the SIA-index is NP-hard to compute, even if all matrices in have a positive diagonal, and to approximate within a factor of . Matrices with positive diagonals appear in consensus applications and many problems are computationally easier for these matrices .

Finally, we show that the st power of any SIA automata matrix has a positive column and that the rd power of any SIA matrix has a positive column.

We mention that the majority of the presented results only depend on the pattern of nonzero elements in the matrices. Therefore, they also can be applied for row-allowable matrices, i.e., nonnegative matrices that have at least one positive element on each row.

Paper Organization. In Section 2 we give an overview of commonly used matrix classes: Positive-column, scrambling, Sarymsakov and SIA matrices (Subsection 2.1) and show the relations between them (Subsection 2.2). We study the associated indices and bounds on these indices (Subsection 2.3) and derive an upper bound on the power at which an SIA matrix has a positive column (Subsection 2.4). Section 3 is devoted to the properties of the SIA-index. In Subsection 3.1, lower and upper bounds on the largest SIA-index are provided, and we present an experiment that computes the largest SIA-index for small values of . In Subsection 3.2 we discuss the procedure to decide whether a given set of matrices is SIA and the hardness of computing and approximating the SIA-index.

## 2 Common Matrix Classes and their Indices

In this section, we define different matrix classes that have been introduced in the context of Markov chains and probabilistic automata [30, 25, 33]: positive-column matrices, scrambling matrices, Sarymsakov matrices and SIA matrices. We show that these classes are included in one another and that sufficiently large powers of matrices from the larger classes belong to the smaller classes (for example, any sufficiently large power of a Sarymsakov matrix is scrambling).

We also define the associated indices – the lengths of the shortest products belonging to the corresponding classes. We relate these indices to each other and present upper bounds on them.

### 2.1 Definition of the Matrix Classes

We start by defining the rather intuitive notion of a positive-column matrix. These matrices are called Markov matrices in certain sources [30, Definition 4.7], but we avoid this terminology as the term Markov matrix is sometimes used to denote a stochastic matrix. Positive-column matrices have a direct interpretation in many applications. For instance, in the consensus system (1), an element of the transition matrix represents the influence of the agent corresponding to the column on that of the row, and a positive-column transition matrix thus corresponds to a situation in which an agent influences all the others. Additionally, positive-column matrices are related to convergence, as the sequence of powers of a positive-column stochastic matrix converges to a rank-one matrix.

###### Definition 1 (Positive-column matrix).

A positive-column matrix is a nonnegative matrix that has a positive column. We will denote the set of positive-column matrices by .

We now define scrambling matrices, a class that has been extensively investigated due to their applications to Markov chains [30, 25].

###### Definition 2 (Scrambling matrix ).

A nonnegative matrix is called scrambling if for any pair of rows , there is a column such that and . We will denote the set of scrambling matrices by .

Scrambling matrices have two interesting properties:

• the class is closed under multiplication: the product of two scrambling matrices is scrambling;

• for any there is a length such that any product of or more scrambling matrices has a positive column.

The Sarymsakov class of matrices is a larger class that has the same desirable properties .

###### Definition 3 (Sarymsakov matrix ).

We define the consequent function of a given nonnegative matrix . For any , the consequent function is equal to

 FP(S)={j:∃i∈S s.t. aij>0}.

A nonnegative matrix is called a Sarymsakov matrix if for any two disjoint nonempty subsets and ,

 FP(S)∩FP(S′)≠∅ (2)

or

 |FP(S)∪FP(S′)|>|S∪S′|. (3)

We will denote the set of Sarymsakov matrices by .

This somewhat formal definition of a Sarymsakov matrix can be interpreted in terms of the graph associated with the matrix , that is, the directed graph whose adjacency matrix is . The consequent function returns the set of out-neighbors of a set of nodes . Condition (2) corresponds to intersecting sets of neighbours and Condition (3) corresponds to “expansion” of sets of neighbours: the union of contains more elements than the union of the sets .

###### Definition 4 (SIA Matrix).

A matrix is called SIA, i.e., Stochastic, Indecomposable and Aperiodic, if it is stochastic and

 Q=limn→∞Pn

exists and all the rows of are the same. We will denote the set of SIA matrices by .

SIA matrices are called regular matrices in .

### 2.2 Relation between the four Matrix Classes

The four matrix classes that we have defined are included in one another:

 SPC⊂SSCR⊂SSAR⊂SSIA. (4)

The first inclusion follows from the definitions and the second and third inclusions are proved in [30, Section Bibliography and Discussion to §§4.3–4.4] and [34, Section II].

We present an alternative characterization of SIA matrices.

###### Proposition 1 (Characterization of SIA Matrices).

A stochastic matrix is SIA if and only if there is such that has a positive column.

###### Proof.

Only if. Let be SIA. We prove that there is such that has a positive column. The exists by definition of an SIA matrix. If there is no such that has a positive column, then has no positive column. On the other hand, all the rows of are the same and the sum of the elements on each row is equal to 1, so that has a positive column and we have a contradiction.

If. Let be such that is positive-column. By inclusions (4), is SIA because positive-column matrices are SIA. Hence

 limt→+∞(Pp)t

exists and has all rows equal, and therefore

 limt→+∞Pt

exists and has all rows equal and is SIA. ∎

This simple characterization has several consequences. It means that only the pattern of zero/nonzero elements determines if a matrix is SIA or not. Another consequence of Proposition 1 is the next proposition.

###### Proposition 2.

For any set of stochastic matrices the four statements are equivalent.

• There exists an SIA product of matrices from

• there exists a Sarymsakov product of matrices from

• there exists a scrambling product of matrices from

• there exists a positive-column product of matrices from .

###### Proof.

This is a consequence of Proposition 1 and of the inclusions (4). ∎

We define the notion of an SIA set as a set of stochastic matrices that has an SIA product (and therefore a Sarymsakov, a scrambling and a positive-column product).

We now present a particular class of stochastic matrices that has been studied extensively in the context of automata theory . We will call them automaton matrices in this article.

###### Definition 5 (automaton matrix).

A matrix is called automaton matrix if it is stochastic and has exactly one 1 on each row.

In the case of automaton matrices, the notions of Sarymsakov, scrambling and positive-column coincide, as shown in the next proposition.

###### Proposition 3.

Let be an automaton matrix. The three properties are equivalent:

1. is Sarymsakov

2. is scrambling

3. is positive-column.

###### Proof.

The inclusions in Equation (4) imply (iii) (ii) and (ii) (i), so that (i) (iii) remains to be proved. Let be Sarymsakov and let . Condition (3) cannot be satisfied for sets and because

 |FP({i})∪FP({j})|>|{i}∪{j}|=2

would imply

 |FP({i})|>1 or |FP({j})|>1,

and at least one of the rows and would have more than one positive element, which is impossible by the definition of an automaton matrix. Therefore, Condition (2) is satisfied for any pair of singletons , meaning that for any rows , the positive element in row is in the same column as the positive element in row , so that the matrix is in fact positive-column.

### 2.3 Indices and Bounds

In the present section we will associate a natural combinatorial parameter with each of the previously discussed classes and provide bounds on them that depend only on the size of the matrices.

###### Definition 6 (X-index).

Let be a class of matrices. The -index of a set of square matrices is the smallest such that there is a product of length of matrices from belonging to .

The -index can be seen as the length of the shortest witness proving that a given set of matrices does actually belong to the class . If no product of belongs to , we will agree that the -index is undefined. In this section, we will focus on SIA-, Sarymsakov-, scrambling- and positive-column-indices of SIA sets and we will denote the corresponding index of a set by , , , . We note that our terminology is in agreement with the case of a single matrix, where the scrambling-index has received a lot of attention .

Observe that for every SIA set , we have the following inequalities as a direct consequence of the inclusions (4):

 sia(S)≤sar(S)≤scr(S)≤pc(S). (5)

One of the basic question arising in regard with these indices is how large the index of a set of stochastic matrices can be? Such questions have received a lot of attention for the exponent [7, Section 3.5] that can be seen as the -index, where is the class of positive matrices. Likewise, the positive-column-index has been studied extensively in the particular case of automaton matrices. It is called reset threshold in this context.

###### Theorem 1.

Let , , , be the largest values of the corresponding indices among all SIA sets of matrices. For any dimension , we have

 sia(n)≤sar(n)=scr(n)=pc(n). (6)
###### Proof.

First, notice that we have the following inequalities a direct consequence of (5):

 sia(n)≤sar(n)≤scr(n)≤pc(n).

It remains to prove that . Let (resp. ) be the largest positive-column-index (resp. Sarymsakov-index) among all sets of automaton matrices.

 pc(n) =pcA(n) (7) =sarA(n) (8) ≤sar(n). (9)

The first equality (7) is Theorem 4 of . The second (8) is a consequence of the fact that an automaton matrix is Sarymsakov if and only if it is positive-column (see Proposition 3). Inequality (9) holds true trivially, since automaton matrices form by definition a subset of stochastic matrices. ∎

### 2.4 From SIA to a Positive Column

In this section, we will provide upper bounds on the powers at which SIA matrices are guaranteed to have a positive column. First, we will prove in Proposition 4 that the st power of any SIA automata matrix has a positive column. Then, we will show that the power of any arbitrary SIA matrix has a positive column and present a finer bound that depends on the number of columns that are eventually positive (Theorem 2.4).

We will say that a matrix and a graph are associated if the matrix is the adjacency matrix of the graph.

###### Proposition 4.

Let be an SIA automaton matrix. Its st power has a positive column. Moreover, the value cannot be decreased (in general): for any there exists an SIA automaton matrix whose power has no positive column.

###### Proof.

By Proposition 1, there exists such that . The entry of matrix counts the number of walks of length from node to node in the graph associated to , and the positive column means that there is a walk of length from any node to the node corresponding to the positive column.

Claim: has no cycle other than a self-loop. If there is a cycle of length , the submatrix corresponding to the nodes of the cycle is a permutation matrix. Up to reordering of the nodes, the matrix is then equal to

 (AB0C),

with a permutation matrix. Then is not positive-column and we have a contradiction. Therefore there is no cycle of length . Furthermore, has a self-loop because there is no other cycle and all node have one outgoing edge.

It is now clear that the positive column of is the column corresponding to the node that has a self-loop. In , there is a walk of length at most from any node to , and because node has a self-loop, there is a walk of length exactly .

For the tightness part, the following matrix is an example of SIA automaton matrix whose nd power has no positive column and its has one.

 ⎛⎜ ⎜ ⎜ ⎜ ⎜ ⎜ ⎜⎝111⋱10⎞⎟ ⎟ ⎟ ⎟ ⎟ ⎟ ⎟⎠.

This last proposition has possible consequences for the Černý conjecture because it means that is a bound on the reset threshold555Recall that the reset threshold is the positive-column index of a set of automaton matrices. of synchronizing automata and therefore, any subquadratic bound would improve the state-of-the-art upper bound on the Černý conjecture.

Proposition 4 gives a tight upper bound on the power at which an SIA automata matrix has a positive column. We will now generalize this result to all SIA matrices. In order to do so, we will rely on the local exponents of primitive matrices. Recall that a square matrix is primitive if (entrywise) for some natural . The smallest such is known as the exponent of  [7, Section 3.5]. The local exponents are refinements of this characteristic : the th local exponent of a primitive matrix is the smallest power having at least positive rows . Observe that the first local exponent of is the positive-column index of the transpose of and the th local exponent is . We will make use of the following theorem:

###### Theorem 2 ([8, Theorem 3.4]).

The largest value of the th local exponent among primitive matrices, with , is equal to .

Before stating our result, we define the notion of a column that is positive in sufficiently large powers of the matrix.

###### Definition 7 (Eventually Positive Columns).

Let be a stochastic matrix. We say that a column is eventually positive if there is a power such that the th column of is positive.

We can already notice that the if has a positive th column, then has a positive th column as well, and so has a positive th column for any .

###### Theorem 3.

Let be an SIA matrix, and let be the number of eventually positive columns of . The power of has a positive column.

###### Proof.

We can assume that the first columns are eventually positive, and let us partition in blocks.

 P=(ABCD),

with having a size of and blocks , and having sizes , and respectively.

Claim: is primitive and . Let be such that has positive columns. If we partition the matrix in blocks, as we did for , we obtain

 Pp=(EFGH),

again with being . We have that and are entrywise positive because we assumed that the first columns of are positive.

 Pp+1=(EA+FCEB+FDGA+HCGB+HD).

We prove now that . Assume to the contrary that , the element at position in the matrix , is positive. The th column of is equal to

 (EB)j=∑iEibij,

with being the th column of . Because is positive, the th column of is positive, and similarly, the th column of is positive. Therefore, the th column of is positive, and we have a contradiction because only the first columns of are eventually positive. Hence . Since , we have that and we can conclude that is primitive.

Let us now compute . Let us define and as below:

 Pn−k=(An−k0JDn−k)

and

 Pk2−3k+3=(Ak2−3k+30KDk2−3k+3)

The matrix has a positive column thanks to Theorem 2 applied to the transpose of . Let us assume that the th column is positive.

We can notice that both and have a positive element on each row. has a positive element on each row because has a positive element on each row and the first rows of are . And has a positive element on each row because there is a path of length from each node to some node in the graph associated to .

Therefore and have a positive th column and therefore

 Pk2−4k+3+n=(An−kAk2−3k+30JAk2−3k+3+Dn−kKDk2−4k+3+n)

has a positive column. ∎

###### Corollary 1.

The rd power of a SIA matrix has a positive column.

###### Proof.

This is a consequence of the previous theorem and

 ∀n∈N,k∈{1,…,n},k2−4k+3+n≤n2−3n+3.

## 3 The SIA-index

In this section we study the SIA-index. We first show that, for any dimension , the largest SIA-index, , is reached by a set of automaton matrices. We provide general upper and lower bounds on that depend only on . We describe the results of our numerical experiments to compute exactly for small values of . Finally, we conclude with algorithmic problems related to the computation of the SIA index.

The following result highlights the importance of automaton matrices in the context of SIA matrices.

###### Proposition 5.

For any , there is a set of matrices such that and every matrix of is an automaton matrix.

###### Proof.

In order to construct we will use techniques from . Let be a set of stochastic matrices such that (such a set exists by definition of ). Let

 Bn={B∈SAUTO|∃A∈An,B≲A}, (10)

where is the set of automaton matrices, and means that any zero element in must be zero in : .

By definition of , we have . We prove .

Observation 1: The set has an SIA product. The set has a positive-column product. We can use Lemma 1 of , stating that for any set of stochastic matrices that has a positive-column product, the set , defined as in (10) has a positive-column product as well. The positive-column product of the set is thus SIA by inclusions (4).

Observation 2: Any SIA product of has a length larger than or equal to . Indeed, if is an SIA product, replacing each matrix by a matrix such that yields an SIA product of matrices from .

Therefore, . ∎

### 3.1 Bounds on the SIA-index

#### 3.1.1 Upper bounds

Propositions 4 and 5 tell us that there are tight connections between the SIA sets of matrices and synchronizing automata theory. Such connections were already utilized in , where the authors proved that the largest positive column index, , is equal to the largest reset threshold among automata of size . The best known upper bound on the reset threshold of automata of size is [26, 15]. We can combine it with inequality (6) to derive an upper bound on the SIA-index:

 sia(n)≤n3−n6.

In fact, we believe that is much smaller.

###### Conjecture 1.

The SIA-index of a set of stochastic matrices is bounded by .

In the next subsections we will support this conjecture by providing results of computational experiments and analysis of sets that are extremal for the Černý conjecture. We believe that Conjecture 1 offers a new angle on the Černý conjecture and can bring new insights. First, by Proposition 4 the SIA-index of automaton matrices is the shortest such that there exists a positive-column product , where is a product of length . In other words, if Conjecture 1 is true, then a synchronizing automaton of size has a reset threshold at most , which is a significant improvement of the state of the art. Second, the SIA-index tends to be surprisingly small for automata with large reset thresholds highlighting the structural properties of these particular cases: Rystsov’s automata , Černý automata , and other slowly synchronizing automata  have small SIA-indices.

#### 3.1.2 Numerical results

We now present the results of our computational experiments that support our Conjecture 1. Since the bound on SIA-index for automaton matrices is equal to the bound on the SIA-index for stochastic matrices (Proposition 5), we only investigate automaton matrices. We have computed on a computer cluster the SIA-index of all automata made of two matrices up to , and up to for initially connected automata, a notion that we will define soon. The results are summarized in Table 1. We have done the same for all triplets of automaton matrices up to (Table 1), and we obtain exactly the same maximum SIA-indices. The maximum SIA-index grows approximately like , as shown in Figure 2. Examples of sets of pairs of and matrices that have an SIA-index of and are depicted on Figure 3. Figure 2: Maximum SIA-index for pairs of matrices (blue) and the curve y=2n (red). (a) Pair of 8×8 matrices that has an SIA-index of 15. This is the only such pair up to relabelling of the nodes.

Our methodology was the following. For the exhaustive tests, we have enumerated all sets of automaton matrices and for each set we have computed its SIA-index. In order to compute the SIA-index, we have enumerated all matrix products corresponding to Lyndon words (defined below) and looked for the shortest SIA product of this form. Proposition 6 guarantees that the correct SIA-index is computed.

Let us now define Lyndon words. We will treat now the elements of a matrix set as symbols (or letters) that are linearly ordered in an arbitrary fashion, e.g. . Clearly, products of matrices from are in one-to-one correspondence with words – the sequences of symbols . The lexicographic order on the words is defined as follows: if and only if either for some word ; or and for and some words .

###### Definition 8 (Lyndon word).

A cyclic shift of a word is a word of the form if . A non-empty word is Lyndon if it is strictly smaller in the lexicographic order than all of its cyclic shifts.

###### Proposition 6.

Let be an SIA set of matrices such that . There is a Lyndon word of length such that the corresponding product of matrices from is SIA.

###### Proof.

Let be a word corresponding to an SIA product of length . Observe first that any cyclic shift of leads to an SIA product as well: if the powers of the matrix associated with converge to the rank one matrix , then the matrices associated with converge to rank one matrix as . Thus, we can always choose a word corresponding to an SIA product that is smaller or equal than any of its cyclic shifts in the lexicographic order. Secondly, if has the form for some , then clearly corresponds to an SIA product. Since , we immediately conclude that , i.e. is aperiodic. Now we can invoke the classical result about the Lyndon words stating that an aperiodic word that is not larger than any of its cyclic shifts is actually Lyndon [21, Proposition 5.1.2] and the proposition follows. ∎

We have noticed in our tests that all extremal examples are initially connected (in fact, even strongly connected). Therefore, we have decided to analyze larger values of restricted to the case of initially connected automaton matrices. This has allowed us to perform the tests up to instead of .

###### Definition 9 (Initially Connected).

A set of automaton matrices is called initially connected or IC if in the graph associated with the matrix there exists a node such that there is a path from to any node of the graph. In particular, if the graph associated to is strongly connected, the set is initially connected.

#### 3.1.3 Lower bounds

Proposition 4 implies that automaton matrices with a small SIA-index have a small reset threshold as well. Therefore, we focused on automata that are known to be tight for the Černý conjecture. The results are summarized in Fig.4. A list of automata that are known to be tight for the Černý conjecture can be found in .

The Černý family is the only known infinite series of automaton matrices with the shortest positive-column product of length   . For , the Černý set is equal to

 Cn=⎧⎪ ⎪ ⎪⎨⎪ ⎪ ⎪⎩A=⎛⎜ ⎜ ⎜⎝0100001000011000⎞⎟ ⎟ ⎟⎠,B=⎛⎜ ⎜ ⎜⎝0100010000100001⎞⎟ ⎟ ⎟⎠⎫⎪ ⎪ ⎪⎬⎪ ⎪ ⎪⎭.

We will see in the following proposition that its shortest SIA product is equal to , which implies . For small values of , the values computed in Table 1 provide slightly better lower bounds.

###### Proposition 7.

The Černý set of matrices of dimension has an SIA-index of .

###### Proof.

Observe first that has a positive-column and is SIA. Thus, . Furthermore, since the shortest positive-column product of has length , see e.g. , and due to Proposition 4 we have . Therefore, .

It remains to show now that the case is impossible. Assume to the contrary that is an SIA product of length . Therefore, has a positive column (by Proposition 4). By [19, Proposition 4] every positive-column product of has at least of occurrences of and occurrences of , thus, has exactly one occurrence of and occurrences of . Applying Proposition 6 we further conclude that the word corresponds to an SIA product as well, which is not the case.

Now we will analyze a set of matrices derived from the Wielandt series of matrices that have the largest possible exponent among matrices [7, Chapter 3.5]. Matrix sets of this kind often appear in the study of combinatorial characteristics of matrix sets, e.g. generalizations of the exponents  or in the study of positive-column indicies . We define the Wielandt set of automaton matrices as

 Wn=⎧⎪ ⎪ ⎪ ⎪ ⎪ ⎪ ⎪ ⎪⎨⎪ ⎪ ⎪ ⎪ ⎪ ⎪ ⎪ ⎪⎩A=⎛⎜ ⎜ ⎜ ⎜ ⎜ ⎜ ⎜⎝11⋱101⎞⎟ ⎟ ⎟ ⎟ ⎟ ⎟ ⎟⎠,B=⎛⎜ ⎜ ⎜ ⎜ ⎜ ⎜ ⎜⎝11⋱110⎞⎟ ⎟ ⎟ ⎟ ⎟ ⎟ ⎟⎠⎫⎪ ⎪ ⎪ ⎪ ⎪ ⎪ ⎪ ⎪⎬⎪ ⎪ ⎪ ⎪ ⎪ ⎪ ⎪ ⎪⎭,

where the omitted elements are zeros.

###### Proposition 8.

The Wielandt set of matrices has an SIA-index of .

###### Proof.

Recall that the shortest positive-column product of has length  [3, Theorem 2]. Since cannot be strictly smaller than by Proposition 4, we immediately conclude that . This bound is tight, since is the desired SIA product. ∎

### 3.2 Complexity Results

In the present section we address algorithmic problems related to SIA sets and the SIA-index. We assume that the reader is familiar with the computational complexity theory.

#### 3.2.1 Checking the Existence of an SIA Product

Proposition 2 states that for a set , the existences of a positive-column product, of a scrambling product, of a Sarymsakov product, and of an SIA-product are equivalent. Therefore, deciding algorithmically the existence of an SIA-product can be done in polynomial time with techniques developed in [27, Section 5] and in  to decide the existence of a positive-column product.

#### 3.2.2 Computing the SIA-index

Now we are going to show that the problem of finding the SIA-index of a given set of matrices is computationally hard. Furthermore, it is hard even for matrices with positive diagonals – a special case frequently arising in consensus applications. More precisely, we will present a reduction from the Boolean satisfiability problem to the problem of computing the SIA-index of sets of matrices with the following zero pattern:

 ⎛⎜ ⎜ ⎜ ⎜ ⎜⎝+⋅+⋮⋱⋅+⎞⎟ ⎟ ⎟ ⎟ ⎟⎠,

where denotes a positive element and denotes an element that can be either positive or zero. The following lemma establishes that such a matrix is SIA if and only if it is positive-column.

###### Lemma 1.

Let be a stochastic matrix that has a positive diagonal and all elements that belong neither to the diagonal nor to the first column are 0. is SIA if and only if it is positive-column.

###### Proof.

The ”if” part is a consequence of inclusions (4). We prove the ”only if” part. If is SIA, there is such that is positive-column. If , is positive-column and the proof is complete. So we are left with the case . It is clear that only the first column of can be positive and that , so we prove that for any , .

For any , the element is positive only if or is positive. If is positive, the proof is complete. If not, we use and we repeat the same argument, so that either or and we can continue iteratively until we arrive at the conclusion that either or . ∎

###### Theorem 4.

For a given integer and a set of stochastic matrices with positive diagonal, the problem of deciding whether is NP-complete.

###### Proof.

We use a reduction from 3-SAT similar to that of Theorem 8 in .

A 3-SAT instance consists of variables , clauses of the form where each is a literal, that is, either a variable or the negation of a variable. The problem is to determine whether the formula is satisfiable.

Given a 3-SAT formula with variables, we construct a set of matrices that has an SIA product of length smaller or equal to if and only if the 3-SAT formula is satisfiable.

We define a set of matrices. The matrices have size . The matrix , representing literal , has ones on the diagonal, a 1 in position and a 1 in position for every for which assigning to TRUE satisfies clause . The matrix has ones on the diagonal, a 1 in position and a 1 in position if assigning to FALSE satisfies clause .

Theorem 6 of  proves that the formula is satisfiable if and only if has a positive-column product of length at most . By Lemma 1, has a positive-column product of length at most if and only if it has an SIA-product of length at most , which concludes the proof of NP-hardness.

Additionally, it is clear that the problem belongs to NP. ∎

By the same argument we can conclude that the problems of deciding whether there exists a Sarymsakov, scrambling, or positive-column product of length at most are NP-complete as well. The previous NP-completeness result holds for matrices with positive diagonals. If we remove the restriction that the matrices have positive diagonals, we can obtain a stronger infeasibility result.

###### Theorem 5.

For every , it is NP-hard to approximate the SIA-index of sets of automaton matrices of size within a factor .

###### Proof.

Our proof is based on the reduction given in [17, Section 3] showing inapproximability of the reset thresholds, i.e., positive-column indices of automaton matrices. Recall that the classical set cover problem is formalized as follows. Given a set of elements and a collection of subsets of whose union is equal to , i.e., . The set cover problem is to compute a sub-collection of the smallest possible size, whose union still equals to . This problem is computationally hard, namely, for every , there is no approximation algorithm computing the set cover, unless  .

The reduction from the set cover problem to our problem is done as follows. For every set we construct a -matrix of size : if and only if either and , or and . Clearly, is stochastic and automatic. Let . It is easy to see that a product of matrices from is SIA if and only if the st column of is positive. Observe now that the latter condition is equivalent to

 U=⋃i=1,…,ℓTi.

Therefore, is equal to the size of the smallest cover of and the result follows.

## Conclusion

In the present paper we studied SIA sets of stochastic matrices and compared them to other classical matrix classes. We introduced the SIA-index for such sets and related it to Lyndon words and well-studied combinatorial characteristics such as local exponents and the scrambling index. Bounds on the SIA-index in terms of matrix sizes were provided and it was shown that is hard to compute for a given matrix set .

Several problems remain open. First, there is a significant gap between the lower and upper bounds on . As we have seen, this problem is related to the long-standing Černý conjecture from automata theory and might be difficult to resolve completely. Another open problem is to find the best approximation ratio achievable by polynomial-time algorithms computing the SIA-index of a given matrix set. Since our contribution is based on synchronizing automata theory, the methods used in  can potentially be used to establish the exact ratio. We have proved that for each matrix size there is an SIA set of automaton matrices having the largest possible SIA-index among stochastic matrices of size . We wonder whether representing the largest SIA-index as the solution of an optimization problem on the space of stochastic matrix sets can lead to another proof of this result. Such formulation can potentially unify and generalize similar results appearing in [18, 5, 10].

## References

•  M. Akelbek and S. Kirkland. Primitive digraphs with the largest scrambling index. Linear Algebra and its Applications, 430(4):1099 – 1110, 2009.
•  A. S. Alfa and B. H. Margolius. Two classes of time-inhomogeneous markov chains: Analysis of the periodic case. Annals of Operations Research, 160(1):121–137, 2008.
•  D. S. Ananichev, M. V. Volkov, and V. V. Gusev. Primitive digraphs with large exponents and slowly synchronizing automata. Journal of Mathematical Sciences, 192(3):263–278, 2013.
•  S. Anily and A. Federgruen. Ergodicity in parametric nonstationary markov chains: An application to simulated annealing methods. Operations Research, 35(6):867–874, 1987.
•  V. D. Blondel, R. M. Jungers, and A. Olshevsky. On primitivity of sets of matrices. Automatica, 2014.
•  V. D. Blondel and A. Olshevsky. How to decide consensus? A combinatorial necessary and sufficient condition and a proof that consensus is decidable but NP-hard. SIAM Journal on Control and Optimization, 52(5):2707 – 2726, 2014.
•  R. Brualdi and H. Ryser. Combinatorial Matrix Theory. Cambridge University Press, 1991.
•  R. A. Brualdi and Bolian Liu. Generalized exponents of primitive directed graphs. Journal of Graph Theory, 14(4):483–499, 1990.
•  J. Černý. Poznámka k homogénnym eksperimentom s konečnými automatami. Matematicko-fysikalny Casopis SAV, 14:208–216, 1964.
•  P.-Y. Chevalier, J. M. Hendrickx, and R. M. Jungers. Reachability of consensus and synchronizing automata. IEEE Control and Decision Conference, 2015.
•  D. Coppersmith and C. W. Wu. Conditions for weak ergodicity of inhomogeneous markov chains. Statistics & Probability Letters, 78(17):3082 – 3085, 2008.
•  O.L.V. Costa, M.D. Fragoso, and R.P. Marques. Discrete-Time Markov Jump Linear Systems. Applied probability. Springer, 2005.
•  I. Dinur and D. Steurer. Analytical approach to parallel repetition. In

ACM Symposium on Theory of Computing (STOC)

, pages 624–633, 2014.
•  D. Eppstein. Reset sequences for monotonic automata. SIAM Journal on Computing, 19(3):500–510, 1990.
•  P. Frankl. An extremal problem fro two families of sets. European Journal of Combinatorics, 3:125–127, 1982.
•  P. Gawrychowski and D. Straszak. Strong inapproximability of the shortest reset word. In Mathematical Foundations of Computer Science 2015, pages 243–255. Springer Berlin Heidelberg, 2015.
•  M. Gerbush and B. Heeringa. Approximating minimum reset sequences. In Implementation and Application of Automata, volume 6482 of LNCS, pages 154–162. Springer, 2011.
•  B. Gerencsér, V. V. Gusev, and R. M. Jungers. Primitive sets of nonnegative matrices and synchronizing automata. CoRR, abs/1602.07556, 2016.
•  V. V. Gusev. Lower bounds for the length of reset words in eulerian automata. International Journal of Foundations of Computer Science, 24(02):251–262, 2013.
•  A. Jadbabaie and J. Lin. Coordination of groups of mobile autonomous agents using nearest neighbor rules. IEEE Transactions on Automatic Control, 48(6):988–1001, 2003.
•  M. Lothaire. Combinatorics on Words. Cambridge University Press, 1997.
•  F. Maire, R. Douc, J. Olsson, et al.

Comparison of asymptotic variances of inhomogeneous Markov chains with application to Markov chain Monte Carlo methods.

The Annals of Statistics, 42(4):1483–1510, 2014.
•  A. Paz. Graph-theoretic and algebraic characterizations of some Markov processes. Israel Journal of Mathematics, 1(3):169–180, Sep 1963.
•  A. Paz. Definite and quasidefinite sets of stochastic matrices. Proceedings of the American Mathematical Society, 16(4):634–641, 1965.
•  A. Paz. Introduction to Probabilistic Automata. Academic Press, New York, 1971.
•  J.-É. Pin. On two combinatorial problems arising from automata theory. Annals of Discrete Mathematics, 17:535–548, 1983.
•  V. Yu. Protasov and A. S. Voynov. Sets of nonnegative matrices without positive products. Linear Algebra and its Applications, 437:749 – 765, 2012.
•  I. Rystsov. Reset words for commutative and solvable automata. Theoretical Computer Science, 172(1):273 – 279, 1997.
•  T. A. Sarymsakov. Inhomogeneous Markov chains (in Russian). Teor. Verojatnost. i Primen., 6:194 – 201, 1961.
•  E. Seneta. Non-negative Matrices and Markov Chains. Springer-Verlag, New York, 1981.
•  B. L. Shader and S. Suwilo. Exponents of nonnegative matrix pairs. Linear Algebra and its Applications, 363:275 – 293, 2003.
•  M. V. Volkov. Synchronizing automata and the Černý conjecture. In Proceedings of the International Conference on Language and Automata Theory and Applications, pages 11–27. Springer, 2008.
•  J. Wolfowitz. Products of indecomposable, aperiodic, stochastic matrices. Proceedings of the American Mathematical Society, 15:733–736, 1963.
•  W. Xia, J. Liu, M. Cao, K. H. Johansson, and T. Başar. Products of generalized stochastic sarymsakov matrices. IEEE Control and Decision Conference, 2015.