Using Sat solvers for synchronization issues in partial deterministic automata

TL;DR

This paper presents a SAT-based method for finding minimal carefully synchronizing words in partial deterministic automata, demonstrating effectiveness on automata with up to 100 states using modest resources.

Contribution

The paper introduces a novel SAT encoding approach for the synchronization problem in partial automata, enabling practical solutions with limited computational resources.

Findings

Effective for automata with up to 100 states

Satisfactory results with modest computational resources

New encoding approach improves practical solvability

Abstract

We approach the task of computing a carefully synchronizing word of minimum length for a given partial deterministic automaton, encoding the problem as an instance of SAT and invoking a SAT solver. Our experimental results demonstrate that this approach gives satisfactory results for automata with up to 100 states even if very modest computational resources are used.

Equations

$$\delta (S,w):=\{\begin{array}{ll}S& \text{ if }w=\epsilon ,\\ {\bigcup }_{q\in \delta (S,v)}\delta (q,a)& \text{ if w=va with v∈Σ∗ and a∈Σ.}\end{array}$$

$$q.a=\{\begin{array}{ll}qa& \text{if qa is a proper prefix of a word of X},\\ \epsilon & \text{if qa∈X},\\ \text{undefined}& \text{otherwise}.\end{array}$$

$$x_{1,t}\vee \cdots \vee x_{m,t},\neg x_{r,t}\vee \neg x_{s,t},\text{ where }1\le r<s\le m.$$

$$\neg y_{j,t-1}\vee \neg x_{i,t}\vee y_{k,t}.$$

$$\neg y_{j,t-1}\vee \neg x_{i,t}.$$

$$\neg y_{r,\ell }\vee \neg y_{s,\ell },\text{ where }1\le r<s\le n.$$

$$\begin{gathered} n+\ell \left(\frac{m(m-1)}{2}+1\right)+\ell mn+\frac{n(n-1)}{2}= \\ \ell \left(\frac{m(m-1)}{2}+mn+1\right)+\frac{n(n+1)}{2}. \end{gathered}$$

$$\phi (x_{i,t})$$

$$\phi (y_{j,0})$$

$$\phi (y_{j,t})$$

$$\ell (n)\approx 3.92+0.49n-0.005n^2+0.000024n^3.$$

$$q.a:=\{\begin{array}{ll}q+1& \text{if }q=1,2,\\ q& \text{if }q=3,\dots ,n;\end{array}q.b:=\{\begin{array}{ll}\text{undefined}& \text{if }q=1,\\ q+1& \text{if }q=2,\dots ,n-1,\\ 1& \text{if }q=n.\end{array}$$

$$P.a:=\{\begin{array}{ll}\{q.a\mid q\in P\}& \text{provided q.a is defined for all q∈P},\\ \text{undefined}& \text{otherwise}.\end{array}$$

Full text

1 Introduction

A nondeterministic finite automaton (NFA) is a triple $\langle Q,\Sigma,\delta\rangle$, where $Q$ and $\Sigma$ are finite non-empty sets called the state set and the input alphabet respectively, and $\delta$ is a subset of $Q\times\Sigma\times Q$. The elements of $Q$ and $\Sigma$ are called states and letters, respectively, and $\delta$ is referred to as the transition relation111The conventional concept of an NFA includes distinguishing two non-empty subsets of $Q$ consisting of initial and final states. As these play no role in our considerations, the above simplified definition well suffices for the purpose of this paper.. For each pair $(q,a)\in Q\times\Sigma$, we denote by $\delta(q,a)$ the subset $\{q^{\prime}\mid(q,a,q^{\prime})\in\delta\}$ of $Q$; this way $\delta$ can be viewed as a function $Q\times\Sigma\to\mathcal{P}(Q)$, where $\mathcal{P}(Q)$ is the power set of $Q$. When we treat $\delta$ as a function, we refer to it as the transition function.

Let $\Sigma^{*}$ stand for the collection of all finite words over the alphabet $\Sigma$, including the empty word $\varepsilon$. The transition function extends to a function $\mathcal{P}(Q)\times\Sigma^{*}\to\mathcal{P}(Q)$, still denoted $\delta$, in the following inductive way: for every subset $S\subseteq Q$ and every word $w\in\Sigma^{*}$, we set

[TABLE]

(Here the set $\delta(S,v)$ is defined by the induction assumption since $v$ is shorter than $w$.) We say that a word $w\in\Sigma^{*}$ is undefined at a state $q\in Q$ if the set $\delta(q,w)$ is empty; otherwise $w$ is said to be defined at $q$.

When we deal with a fixed NFA, we suppress the sign of the transition relation, introducing the NFA as the pair $\langle Q,\Sigma\rangle$ rather than the triple $\langle Q,\Sigma,\delta\rangle$ and writing $q.w$ for $\delta(q,w)$ and $S.w$ for $\delta(S,w)$.

A partial (respectively, complete) * deterministic* automaton is an NFA $\langle Q,\Sigma\rangle$ such that $|q.a|\leq 1$ (respectively, $|q.a|=1$) for all $(q,a)\in Q\times\Sigma$. We use the acronyms PFA and CFA for the expressions ‘partial deterministic automaton’ and ‘complete deterministic automaton’, respectively.

A CFA $\mathrsfs{A}=\langle Q,\Sigma\rangle$ is called synchronizing if there exists a word $w\in\Sigma^{*}$ whose action leaves the automaton in one particular state no matter at which state in $Q$ it is applied: $q.w=q^{\prime}.w$ for all $q,q^{\prime}\in Q$. Any $w$ with this property is said to be a synchronizing word for the automaton.

Synchronizing automata serve as simple yet adequate models of error-resistant systems in many applied areas (system and protocol testing, information coding, robotics). At the same time, synchronizing automata surprisingly arise in some parts of pure mathematics and theoretical computer science (symbolic dynamics, theory of substitution systems, formal language theory). We refer to the survey [37] and the chapter [18] of the forthcoming ‘Handbook of Automata Theory’ for a discussion of synchronizing automata as well as their diverse connections and applications. From both applied and theoretical viewpoints, the key question is to find the optimal, i.e., shortest reset word for a given synchronizing automaton. Under standard assumptions of complexity theory, this optimization question is known to be computationally hard; see [18, Section 2] for a summary of various hardness results in the area. As it is quite common for hard problems of applied importance, there have been many attempts to develop practical approaches to the question. These approaches have been based on certain heuristics [15, 1, 16] and/or popular techniques, including (but not limiting to) binary decision diagrams [27], genetic and evolutionary algorithms [30, 17], satisfiability solvers [36], answer set programming [10], hierarchical classifiers [28], and machine learning [29].

The present authors [35, 34] have initiated an extension to the realm of NFAs of the approach of [36]. Here we consider a more restricted class, namely, that of PFAs, where studying synchronization issues appears to be much better motivated. While we follow the general strategy of and re-use some technical tricks from [35, 34], our present constructions heavily depend on the specifics of partial automata and have not been obtained via specializing the constructions of those earlier papers.

The rest of the paper is structured as follows. In Sect. 2 we describe and motivate the version of PFA synchronization that we have studied. In Sect. 3 we first outline the approach based on satisfiability solvers and then explain in detail how we encode PFAs and their synchronization problems as instances of the Boolean satisfiability problem. In Sect. 4 we provide samples of our experimental results and conclude in Sect. 5 with a brief discussion of the future work.

We have tried to make the paper, to a reasonable extent, self-contained, except for a few discussions that involve some basic concepts of computational complexity theory. These concepts can be found, e.g., in the early chapters of the textbook [26].

2 Synchonization of NFAs and PFAs

The concept of synchronization of CFAs as defined in Sect. 1 was extended to NFAs in several non-equivalent ways. The following three nowadays popular versions were suggested and analyzed in [11] in 1999 (although, in an implicit form, some of them appeared in the literature much earlier, see, e.g., [4, 9]). For $i\in\{1,2,3\}$, an NFA $\mathrsfs{A}=\langle Q,\Sigma\rangle$ is called $D_{i}$-synchronizing if there exists a word $w\in\Sigma^{*}$ that satisfies the condition $(D_{i})$ from the list below:

1. $(D_{1})$:

$|q.w|=|Q.w|=1$ for all $q\in Q$; 2. $(D_{2})$:

$q.w=Q.w$ for all $q\in Q$; 3. $(D_{3})$:

$\bigcap_{q\in Q}q.w\neq\varnothing$.

Any word satisfying $(D_{i})$ is called $D_{i}$-synchronizing for $\mathrsfs{A}$. The definition readily yields the following properties of $D_{i}$-synchronizing words:

Lemma 1

$a)$* A $D_{1}$- or $D_{3}$-synchronizing word is defined at each state.*

$b)$* A $D_{2}$-synchronizing word is either defined at each state or undefined at each state.*

$c)$* Every $D_{1}$-synchronizing word is both $D_{2}$- and $D_{3}$-synchronizing; every $D_{2}$-synchronizing word defined at each state is $D_{3}$-synchronizing.*

In [35] we adapted the approach based on satisfiability solvers to finding $D_{3}$-synchronizing words of minimum length for NFAs. The first-named author used a similar method in the cases of $D_{1}$- and $D_{2}$-synchronization; results related to $D_{2}$-synchronization were reported in [34].

Yet another version of synchronization for NFAs was introduced in [13] and systematically studied in [21, 22, 23, 24, 25], which terminology we adopt. An NFA $\mathrsfs{A}=\langle Q,\Sigma\rangle$ is called carefully synchronizing if there is a word $w=a_{1}\cdots a_{\ell}$, with $a_{1},\dots,a_{\ell}\in\Sigma$, that satisfies the condition $(C)$, being the conjunction of $(C1)$–$(C3)$ below:

1. $(C1)$:

the letter $a_{1}$ is defined at every state in $Q$; 2. $(C2)$:

the letter $a_{t}$ with $1<t\leq\ell$ is defined at every state in $Q.a_{1}\cdots a_{t-1}$, 3. $(C3)$:

$|Q.w|=1$.

Any $w$ satisfying $(C)$ is called a carefully synchronizing word (c.s.w., for short) for $\mathrsfs{A}$. Thus, when a c.s.w. is applied at any state in $Q$, no undefined transition occurs during the course of application. Every carefully synchronizing word is clearly $D_{1}$-synchronizing but the converse is not true in general; moreover, a $D_{1}$-synchronizing NFA may admit no c.s.w.

In this paper we focus on carefully synchronizing words for PFAs. There are several theoretical and practical reasons for this.

On the theoretical side, it is easy to see that each of the conditions $(C)$, $(D_{1})$, $(D_{3})$ leads to the same notion when restricted to PFAs. As for $D_{2}$-synchronization, if a word $w$ is $D_{2}$-synchronizing for a PFA $\mathrsfs{A}$, then $w$ carefully synchronizes $\mathrsfs{A}$ whenever $w$ is defined at each state. Otherwise $w$ is nowhere defined by Lemma 1b, and such ‘annihilating’ words are nothing but usual synchronizing words for the CFA obtained from $\mathrsfs{A}$ by adding a new sink state and making all transitions undefined in $\mathrsfs{A}$ lead to this sink state. Synchronization of CFAs with a sink state is relatively well understood (see [33]), and therefore, we may conclude that $D_{2}$-synchronization also reduces to careful synchronization in the realm of PFAs. On the other hand, there exists a simple transformation that converts every NFA $\mathrsfs{A}=\langle Q,\Sigma\rangle$ into a PFA $\mathrsfs{B}=\langle Q,\Sigma^{\prime}\rangle$ such that $\mathrsfs{A}$ is $D_{3}$-synchronizing if and only if so is $\mathrsfs{B}$ and the minimum lengths of $D_{3}$-synchronizing words for $\mathrsfs{A}$ and $\mathrsfs{B}$ are equal; see [12, Lemma 8.3.8] and [14, Lemma 3]. These observations demonstrate that from the viewpoint of optimal synchronization, studying carefully synchronizing words for PFAs may provide both lower and upper bounds applicable to arbitrary NFAs and all aforementioned kinds of synchronization.

Probably even more important is the fact that careful synchronization of PFAs is relevant in numerous applications. Due to the page limit, we mention only two examples here.

In industrial robotics, synchronizing automata are widely used to design feeders, sorters, and orienters that work with flows of certain objects carried by a conveyer. The goal is achieved by making the flow encounter passive obstacles placed appropriately along the conveyer belt; see [19, 20] for the origin of this automata approach and [2] for an illustrative example. Now imagine that the objects to be oriented or sorted have a fragile side that could be damaged if hitting an obstacle. In order to prevent any damage, we have to forbid ‘dangerous’ transitions in the automaton modelling the orienter/sorter so that the automaton becomes partial and the obstacle sequences must correspond to carefully synchronizing words. (Actually, the term ‘careful synchronization’ has been selected with this application in mind.)

Our second example relates to so-called synchronized codes222We refer the reader to [3, Chapters 3 and 10] for a detailed account of profound connections between codes and automata.. Recall that a prefix code over a finite alphabet $\Sigma$ is a set $X$ of words in $\Sigma^{*}$ such that no word of $X$ is a prefix of another word of $X$. Decoding of a finite prefix code $X$ over $\Sigma$ can be implemented by a finite deterministic automaton $\mathrsfs{A}_{X}$ whose state $Q$ is the set of all proper prefixes of the words in $X$ (including the empty word $\varepsilon$) and whose transitions are defined as follows: for $q\in Q$ and $a\in\Sigma$,

[TABLE]

In general, $\mathrsfs{A}_{X}$ is a PFA (it is complete if and only if the code $X$ is not contained in another prefix code over $\Sigma$). It can be shown that if $\mathrsfs{A}_{X}$ is carefully synchronizing, the code $X$ has a very useful property: whenever a loss of synchronization between the decoder and the coder occurs (because of a channel error), it suffices to transmit a c.s.w. $w$ of $\mathrsfs{A}_{X}$ such that $w$ sends all states in $Q$ to the state $\varepsilon$ to ensure that the following symbols will be decoded correctly.

We may conclude that the problems of determining whether or not a given PFA is carefully synchronizing and of finding its shortest carefully synchronizing words are both natural and important. The bad news is that these problems turn out to be quite difficult: it is known that the first problem is PSPACE-complete and that the minimum length of carefully synchronizing words for carefully synchronizing PFAs can be exponential as a function of the number of states. (These results were found in [31, 32] and later rediscovered and strengthened in [23].) Thus, one has to use some tools that have proved to be efficient for dealing with computationally hard problems. As mentioned in Section 1, in this paper we make an attempt to employ a satisfiability solver as such a tool.

3 Encoding

For completeness, recall the formulation of the Boolean satisfiability problem (SAT). An instance of SAT is a pair $(V,C)$, where $V$ is a set of Boolean variables and $C$ is a collection of clauses over $V$. (A clause over $V$ is a disjunction of literals and a literal is either a variable in $V$ or the negation of a variable in $V$.) Any truth assignment on $V$, i.e., any map $\varphi\colon V\to\{0,1\}$, extends to a map $C\to\{0,1\}$ (still denoted by $\varphi$) via the usual rules of propositional calculus: $\varphi(\neg x)=1-\varphi(x)$, $\varphi(x\vee y)=\max\{\varphi(x),\varphi(y)\}$. A truth assignment $\varphi$ satisfies $C$ if $\varphi(c)=1$ for all $c\in C$. The answer to an instance $(V,C)$ is YES if $(V,C)$ has a satisfying assignment (i.e., a truth assignment on $V$ that satisfies $C$) and NO otherwise.

By Cook’s classic theorem (see, e.g., [26, Theorem 8.2]), SAT is NP-complete, and by the very definition of NP-completeness, every problem in NP reduces to SAT. On the other hand, over the last score or so, many efficient SAT-solvers, i.e., specialized programs designed to solve instances of SAT have been developed. Modern SAT solvers can solve instances with hundreds of thousands of variables and millions of clauses within a few minutes. Due to this progress, the following approach to computationally hard problems has become quite popular nowadays: one encodes instances of such problems into instances of SAT that are then fed to a SAT-solver333We refer the reader to the survey [8] a detailed discussion of the approach and examples of its successful applications in various areas.. It is exactly the strategy that we want to apply.

We start with the following problem:

CSW (the existence of a c.s.w. of a given length):

Input: a PFA $\mathrsfs{A}$ and a positive integer $\ell$ (given in unary);

Output: YES if $\mathrsfs{A}$ has a c.s.w. of length $\ell$;

Output: NO otherwise.

We have to assume that the integer $\ell$ is given in unary because with $\ell$ given in binary, a polynomial time reduction from CSW to SAT is hardly possible. (Indeed, it easily follows from [23] that the version of CSW in which the integer parameter is given in binary is PSPACE-hard, and the existence of a polynomial reduction from a PSPACE-hard problem to SAT would imply that the polynomial hierarchy collapses at level 1.) In contrast, the version of CSW with the unary integer parameter is easily seen to belong to NP: given an instance $(\mathrsfs{A}=\langle Q,\Sigma\rangle,\ell)$ of CSW in this setting, guessing a word $w\in\Sigma^{*}$ of length $\ell$ is legitimate. Then one just checks whether or not $w$ is carefully synchronizing for $\mathrsfs{A}$, and time spent for this check is clearly polynomial in the size of $(\mathrsfs{A},\ell)$.

Now, given an arbitrary instance $(\mathrsfs{A},\ell)$ of CSW, we construct an instance $(V,C)$ of SAT such that the answer to $(\mathrsfs{A},\ell)$ is YES if and only if so is the answer to $(V,C)$. In the following presentation of our encoding, precise definitions and statements are interwoven with less formal comments explaining the ‘physical’ meaning of variables and clauses.

So, take a PFA $\mathrsfs{A}=\langle Q,\Sigma\rangle$ and an integer $\ell>0$. Denote the sizes of $Q$ and $\Sigma$ by $n$ and $m$ respectively, and fix some numbering of these sets so that $Q=\{q_{1},\dots,q_{n}\}$ and $\Sigma=\{a_{1},\dots,a_{m}\}$.

We start with introducing the variables used in the instance $(V,C)$ of SAT that encodes $(\mathrsfs{A},\ell)$. The set $V$ consists of two sorts of variables: $m\ell$ letter variables $x_{i,t}$ with $1\leq i\leq m$, $1\leq t\leq\ell$, and $n(\ell+1)$ state variables $y_{j,t}$ with $1\leq j\leq n$, $0\leq t\leq\ell$. We use the letter variables to encode the letters of a hypothetical c.s.w. $w$ of length $\ell$: namely, we want the value of the variable $x_{i,t}$ to be 1 if and only if the $t$-th letter of $w$ is $a_{i}$. The intended meaning of the state variables is as follows: we want the value of the variable $y_{j,t}$ to be 1 whenever the state $q_{j}$ belongs to the image of $Q$ under the action of the prefix of $w$ of length $t$, in which situation we say that $q_{j}$ is active after $t$ steps. We see that the total number of variables in $V$ is $m\ell+n(\ell+1)=(m+n)\ell+n$.

Now we turn to constructing the set of clauses $C$. It consists of four groups. The group $I$ of initial clauses contains $n$ one-literal clauses $y_{j,0}$, $1\leq j\leq n$, and expresses the fact that all states are active after 0 steps.

For each $t=1,\dots,\ell$, the group $L$ of letter clauses includes the clauses

[TABLE]

Clearly, the clauses (1) express the fact that the $t$-th position of our hypothetical c.s.w. $w$ is occupied by exactly one letter in $\Sigma$. Altogether, $L$ contains $\ell\left(\frac{m(m-1)}{2}+1\right)$ clauses.

For each $t=1,\dots,\ell$ and each triple $(q_{j},a_{i},q_{k})$ in the transition relation of $\mathrsfs{A}$, the group $T$ of transition clauses includes the clause

[TABLE]

Invoking the basic laws of propositional logic, one sees that the clause (2) is equivalent to the implication $y_{j,t-1}\mathop{\&}x_{i,t}\to y_{k,t},$ that is, (2) expresses the fact that if the state $q_{j}$ has been active after $t-1$ steps and $a_{i}$ is the $t$-th letter of $w$, then the state $q_{k}=q_{j}.a_{i}$ becomes active after $t$ steps. Further, for each $t=1,\dots,\ell$ and each pair $(q_{j},a_{i})$ such that $a_{i}$ is undefined at $q_{j}$ in $\mathrsfs{A}$, we add to $T$ the clause

[TABLE]

The clause is equivalent to the implication $y_{j,t-1}\to\neg x_{i,t},$ and thus, it expresses the requirement that the letter $a_{i}$ should not be occur in the $t$-th position of $w$ if $q_{j}$ has been active after $t-1$ steps. Obviously, this corresponds to the conditions $(C1)$ (for $t=0$) and $(C2)$ (for $t>0$) in the definition of careful synchronization. For each $t=1,\dots,\ell$ and each pair $(q_{j},a_{i})\in Q\times\Sigma$, exactly one of the clauses (2) or (3) occurs in $T$, whence $T$ consists of $\ell mn$ clauses.

The final group $S$ of of synchronization clauses includes the clauses

[TABLE]

The clauses (4) express the requirement that at most one state remains active when the action of the word $w$ is completed, which corresponds to the condition $(C3)$ from the definition of careful synchronization. The group $S$ contains $\frac{n(n-1)}{2}$ clauses.

Summing up, the number of clauses in $C:=I\cup L\cup T\cup S$ is

[TABLE]

In comparison with encodings used in our earlier papers [35, 34], the encoding suggested here produces much smaller SAT instances. Since in the applications the size of the input alphabet is a (usually small) constant, the leading term in (5) is $\Theta(\ell n)$ while the restriction to PFAs of the encodings from [35, 34] has $\Theta(\ell n^{2})$ clauses.

Theorem 2

A PFA $\mathrsfs{A}$ has a c.s.w. of length $\ell$ if and only if the instance $(V,C)$ of SAT constructed above is satisfiable. Moreover, the carefully synchronizing words of length $\ell$ for $\mathrsfs{A}$ are in a 1-1 correspondence with the restrictions of satisfying assignments of $(V,C)$ to the letter variables.

Proof

Suppose that $\mathrsfs{A}$ has a c.s.w. of length $\ell$. We fix such a word $w$ and denote by $w_{t}$ its prefix of length $t=1,\dots,\ell$. Define a truth assignment $\varphi\colon V\to\{0,1\}$ as follows: for $1\leq i\leq m$, $0\leq j\leq n$, $1\leq t\leq\ell$, let

[TABLE]

In view of (6) and (7), $\varphi$ satisfies all clauses in $L$ and respectively $I$. As $w_{\ell}=w$ and $|Q.w|=1$, we see that (8) ensures that $\varphi$ satisfies all clauses in $S$. It remains to analyze the clauses in $T$. For each fixed $t=1,\dots,\ell$, these clauses are in a 1-1 correspondence with the pairs in $Q\times\Sigma$. We fix such a pair $(q_{j},a_{i})$, denote the clause corresponding to $(q_{j},a_{i})$ by $c$ and consider three cases.

Case 1: the letter $a_{i}$ is not the $t$-th letter of $w$. In this case $\varphi(x_{i,t})=0$ by (6), and hence, $\varphi(c)=1$ since the literal $\neg x_{i,t}$ occurs in $c$, independently of $c$ having the form (2) or (3).

Case 2: the letter $a_{i}$ is the $t$-th letter of $w$ but it is undefined at $q_{j}$. In this case the clause $c$ must be of the form (3). Observe that $t>1$ in this case since the first letter of the c.s.w. $w$ must be defined at each state in $Q$. Moreover, the state $q_{j}$ cannot belong to the set $Q.w_{t-1}$ because $a_{i}$ must be defined at each state in this state. Hence $\varphi(y_{j,t-1})=0$ by (8), and $\varphi(c)=1$ since the literal $\neg y_{j,t-1}$ occurs in $c$.

Case 3: the letter $a_{i}$ is the $t$-th letter of $w$ and it is defined at $q_{j}$. In this case the clause $c$ must be of the form (2), in which the literal $y_{k,t}$ corresponds to the state $q_{k}=q_{j}.a_{i}$. If the state $q_{j}$ does not belong to the set $Q.w_{t-1}$, then as in the previous case, we have $\varphi(y_{j,t-1})=0$ and $\varphi(c)=1$. If $q_{j}$ belongs to $Q.w_{t-1}$, then the state $q_{k}$ belongs to the set $(Q.w_{t-1}).a_{i}=Q.w_{t}$, whence $\varphi(y_{k,t})=1$ by (8). We conclude that $\varphi(c)=1$ since the literal $y_{k,t}$ occurs in $c$.

Conversely, suppose that $\varphi\colon V\to\{0,1\}$ is a satisfying assignment for $(V,C)$. Since $\varphi$ satisfies the clauses in $L$, for each $t=1,\dots,\ell$, there exists a unique $i\in\{1,\dots,m\}$ such that $\varphi(x_{i,t})=1$. This defines a map $\chi\colon\{1,\dots,\ell\}\to\{1,\dots,m\}$. Let $w:=a_{\chi(1)}\cdots a_{\chi(\ell)}$. We aim to show that $w$ is a c.s.w. for $\mathrsfs{A}$, i.e., to verify that $w$ fulfils the conditions $(C1)$–$(C3)$ from the definition of a c.s.w. For this, we first prove two auxiliary claims. Recall that a state is said to be active after $t$ steps if it lies in $Q.w_{t}$, where, as above, $w_{t}$ is the length $t$ prefix of the word $w$. (By the length 0 prefix we understand the empty word $\varepsilon$.)

Claim 1. For each $t=0,1,\dots,\ell$, there are states active after $t$ steps.

Claim 2. If a state $q_{k}$ is active after $t$ steps, then $\varphi(y_{k,t})=1$.

We prove both claims simultaneously by induction on $t$. The induction basis $t=0$ is guaranteed by the fact that all states are active after 0 steps and $\varphi$ satisfies the clauses in $I$. Now suppose that $t>0$ and there are states active after $t-1$ steps. Let $q_{r}$ be such a state. Then $\varphi(y_{r,t-1})=1$ by the induction assumption. Let $i:=\chi(t)$, that is, $a_{i}$ is the $t$-th letter of the word $w$. Then $\varphi(x_{i,t})=1$, whence $\varphi$ cannot satisfy the clause of the form (3) with $j=r$. Hence this clause cannot appear in $T$ as $\varphi$ satisfies the clauses in $T$. This means that the letter $a_{i}$ is defined at $q_{r}$ in $\mathrsfs{A}$, and the state $q_{s}:=q_{r}.a_{i}$ is active after $t$ steps. Claim 1 is proved.

Now let $q_{k}$ be an arbitrary state that is active after $t>0$ steps. Since $a_{i}$ is the $t$-th letter of $w$, we have $Q.w_{t}=(Q.w_{t-1}).a_{i}$, whence $q_{k}=q_{j}.a_{i}$ for same $q_{j}\in Q.w_{t-1}$. Therefore the clause (2) occurs in $T$, and thus, it is satisfied by $\varphi$. Since $q_{j}$ is active after $t-1$ steps, $\varphi(y_{j,t-1})=1$ by the induction assumption; besides that, $\varphi(x_{i,t})=1$. We conclude that in order to satisfy (2), the assignment $\varphi$ must fulfil $\varphi(y_{k,t})=1$. This completes the proof of Claim 2.

We turn to prove that the word $w$ fulfils $(C1)$ and $(C2)$. This amounts to verifying that for each $t=1,\dots,\ell$, the $t$-th letter of the word $w$ is defined at every state $q_{j}$ that is active after $t-1$ steps. Let, as above, $a_{i}$ stand for the $t$-th letter of $w$. If $a_{j}$ were undefined at $q_{j}$, then by the definition of the set $T$ of transition clauses, this set would include the corresponding clause (3). However, $\varphi(x_{i,t})=1$ by the construction of $w$ and $\varphi(y_{j,t-1})=1$ by Claim 2. Hence $\varphi$ does not satisfy this clause while the clauses from $T$ are satisfied by $\varphi$, a contradiction.

Finally, consider $(C3)$. By Claim 1, some state is active after $\ell$ steps. On the other hand, the assignment $\varphi$ satisfies the clauses in $S$, which means that $\varphi(y_{j,\ell})=1$ for at most one index $j\in\{1,\dots,n\}$. By Claim 2 this implies that at most one state is active after $\ell$ steps. We conclude that exactly one state is active after $\ell$ steps, that is, $|Q.w|=1$. ∎

4 Experimental results

We have successfully applied the encoding constructed in Sect. 3 to solve CSW instances with the help of a SAT-solver. As in [36, 10, 35, 34], we have used MiniSat 2.2.0 [6, 7]. In order to find a c.s.w. of minimum length for a given PFA $\mathrsfs{A}$, we have considered CSW instances $(\mathrsfs{A},\ell)$ with fixed $\mathrsfs{A}$ and performed binary search on $\ell$. Even though our encoding is different form those we used in [35, 34], it shares with them the following useful feature: when presented in DIMACS CNF format, the ‘primary’ SAT instance that encodes the CSW instance $(\mathrsfs{A},1)$ can be easily scaled to the SAT instances that encode the CSW instances $(\mathrsfs{A},\ell)$ with any value of $\ell$. Due to this feature, one radically reduces time needed to prepare the input data for the SAT-solver. We refer the reader to [34, Sect. 3] for a detailed explanation of the trick and an illustrative example.

We implemented the algorithm outlined above in C++ and compiled with GCC 4.9.2. In our experiments we used a personal computer with an Intel(R) Core(TM) i5-2520M processor with 2.5 GHz CPU and 4GB of RAM. The code can be found at https://github.com/hananshabana/SynchronizationChecker.

As a sample of our experimental findings, we present here our results on synchronization of PFAs with a unique undefined transition. Observe that the problem of deciding whether or not a given PFA is carefully synchronizing remains PSPACE-complete even if restricted to this rather special case [23]. We considered random PFAs with $n\leq 100$ states and two input letters. The condition $(C1)$ in the definition of a carefully synchronizing PFA implies that such a PFA must have an everywhere defined letter. We denoted this letter by $a$ and the other letter, called $b$, was chosen to be undefined at a unique state. Further, it is easy to see that for a PFA $\langle Q,\{a,b\}\rangle$ with $a,b$ so chosen to be carefully synchronizing, it is necessary that $|Q.a|<|Q|$. Therefore, we fixed a state $q_{a}\in Q$ and then selected $a$ uniformly at random from all $n^{n-1}$ maps $Q\to Q\setminus\{q_{a}\}$. Similarly, to ensure there is a unique undefined transition with $b$, we fixed a state $q_{b}\in Q$ (not necessarily different from $q_{a}$) and then selected $b$ uniformly at random from all $(n-1)^{n}$ maps $Q\setminus\{q_{b}\}\to Q$. For each fixed $n$, we generated up to 1000 random PFAs this way and calculated the average length $\ell(n)$ of their shortest carefully synchronizing words. We used the least squares method to find a function that best reflects how $\ell(n)$ depends on $n$, and it turned out that our results are reasonably well approximated by the following expression:

[TABLE]

20$$40$$60$$80$$100$$10$$15$$20$$25Number of states n$$\ell(n)ObservedOur estimation

The next graph shows the relation between the relative standard deviation of our datasets and the number of states. We see that the relative standard deviation gradually decreases as the number of states grows.

20$$40$$60$$80$$100$$0.15$$0.2$$0.25$$0.3$$0.35$$0.4Number of statesRelative standard deviation

We performed similar experiments with random PFAs that have two or three undefined transition. We also tested our algorithm on PFAs from the series $\mathrsfs{P}_{n}$ suggested in [5]. The state set of $\mathrsfs{P}_{n}$ is $\{1,2,\dots,n\}$, $n\geq 3$, on which the input letters $a$ and $b$ act as follows:

[TABLE]

We examined all automata $\mathrsfs{P}_{n}$ with $n=4,5,\dots,11$, and for each of them, our result matched the theoretical value predicted by [5, Theorem 3]. The time consumed ranges from 0.301 sec for $n=4$ to 4303 sec for $n=11$. Observe that in the latter case the shortest c.s.w. has length 116 so that honest binary search started with $(\mathrsfs{P}_{11},1)$ required 11 iterations.

We made also a comparison with the only approach to computing carefully synchronizing words of minimum length that we had found in the literature, namely, the approach based on partial power automata. Given a PFA $\mathrsfs{A}=\langle Q,\Sigma\rangle$, its partial power automaton $\mathcal{P}(\mathrsfs{A})$ has the subsets of $Q$ as the states, the same input alphabet $\Sigma$, and the transition function defined as follows: for each $a\in\Sigma$ and each $P\subseteq Q$,

[TABLE]

It is easy to see that $w\in\Sigma^{*}$ is a c.s.w. of minimum length for $\mathrsfs{A}$ if and only if $w$ labels a minimum length path in $\mathcal{P}(\mathrsfs{A})$ starting at $Q$ and ending at a singleton. Such a path can be found by breadth-first search in the underlying digraph of $\mathcal{P}(\mathrsfs{A})$.

The result of the comparison are presented in the picture below. In this experiment we had to restrict to PFAs with at most 16 states since beyond this number of states, our implementation of the method based on partial power automata could not complete the computation due to memory restrictions (recall that we used rather modest computational resources). However, we think that the exhibited data suffice to demonstrate that the approach based on SAT-solvers shows a by far better performance.

6$$8$$10$$12$$14$$16[math]100$$200$$300Number of states $n$time (sec)SAT Partial power automaton

5 Conclusion and future work

We have presented an attempt to approach the problem of computing a c.s.w. of minimum length for a given PFA via the SAT-solver method. For this, we have developed a new encoding, which, in comparison with encodings used in our earlier papers [35, 34], requires a more sophisticated proof but leads to more economic SAT instances.

We plan to continue our experiments. In particular, it is interesting to compare the minimum lengths of a synchronizing word for a synchronizing DFA and of carefully synchronizing words for PFAs that can be obtained from the DFA by removing one or more of its transitions.

We also plan to extend the SAT-solver approach to so-called exact synchronization of PFAs which is of interest for certain applications.