Języki, automaty i obliczenia/Wykład 5: Algorytmy konstrukcji automatu minimalnego

Z Studia Informatyczne
Przejdź do nawigacjiPrzejdź do wyszukiwania

W tym wykładzie podamy algorytmy konstrukcji automatu minimalnego i twierdzenia dowodzące ich poprawności.

Algorytmy konstrukcji automatu minimalnego

Dla języka rozpoznawanego L konstrukcję automatu minimalnego można rozpocząć, startując z opisu języka danego na przykład przez wyrażenie regularne lub też jakiegoś automatu rozpoznającego ten język. W niniejszym wykładzie przedstawimy algorytmy konstrukcji automatu minimalnego obejmujące oba wspomniane punkty startu. Jako pierwszy, nawiązując do rezulatów przedstawionych w poprzednim wykładzie, prezentujemy algorytm, dla którego punktem wyjścia jest język L. Prezentację poprzedzimy wprowadzeniem pewnej operacji na słowach zwanej pochodną J.Brzozowskiego.

Definicja 1.1.

Niech LA* będzie dowolnym językiem, a uA* dowolnym słowem. Pochodną Brzozowskiego (residuum) z języka L względem słowa u nazywamy język

u1L={wA*:uwL}.

Podczas obliczeń pochodnych Brzozowskiego (residuów języka L) można wykorzystać poniższe równości.

Niech L1,L2A* będą dowolnymi językami, aA dowolną literą, a u,vA* dowolnymi słowami. Prawdziwe są następujące równości:

Parser nie mógł rozpoznać (nieznana funkcja „\nonumber”): {\displaystyle \begin{array}{rll}\nonumber u^{-1}(L_1 \cup L_2) & = & u^{-1}L_1 \cup u^{-1}L_2, \\ \nonumber u^{-1}(L_1 \backslash L_2) & = & u^{-1}L_1 \backslash u^{-1}L_2, \\ \nonumber u^{-1}(L_1 \cap L_2) & = & u^{-1}L_1 \cap u^{-1}L_2, \\ \nonumber a^{-1}(L_1L_2) & = & (a^{-1}L_1)L_2 \mbox{, jeśli } 1 \notin L_1, \\ \nonumber a^{-1}(L_1L_2) & = & (a^{-1}L_1)L_2 \cup a^{-1}L_2 \mbox{, jeśli } 1 \in L_1, \\ \nonumber a^{-1}L^* & = & (a^{-1}L)L^*, \\ \nonumber v^{-1}(u^{-1}L) & = & (uv)^{-1}L. \end{array}}

Przykład 1.1.

Obliczmy wszystkie pochodne dla języka L=a+b+. Okazuje się, że są tylko cztery różne pochodne liczone względem a, b, ab i słowa pustego 1. Mianowicie:
a1L=a*b+,
b1L=,
ab1L=b*,
11L=L.
Dla wszystkich innych słów otrzymujemy uzyskane powyżej języki, co wynika z własności pochodnych (patrz wyżej wypisane równości) i z następujacych obliczeń:
n(an)1L=a*b+,
n(bn)1L=,
n(abn)1L=b*.

Zauważmy również, nawiązując raz jeszcze do rezulatów przedstawionych w poprzednim wykładzie, że prawdziwa jest następująca równoważność wiążąca wprowadzone pojęcie pochodnej Brzozowskiego z prawą kongruencją syntaktyczną:

uPrLvu1L=v1L.

Rozpisując z definicji lewą stronę tej równoważności, otrzymujemy, iż dla dowolnego słowa zA* słowo uzL wtedy i tylko wtedy, gdy vzL. A to równoważnie oznacza (znów z definicji), że u1L=v1L.

Z uzasadnionej równoważności oraz twierdzenia 3.4 o prawej kongruencji syntaktycznej z poprzedniego wykładu wnioskujemy równoważność rozpoznawalności języka L i skończonej ilości różnych pochodnych Brzozowskiego tego języka.

Pierwszy z przedstawianych algorytmów będzie konstruował automat minimalny, wyznaczając prawą kongruencję automatową poprzez zastosowanie metody pochodnych Brzozowskiego. Metoda ta umożliwia przeprowadzanie odpowiednich obliczeń bezpośrednio na wyrażeniach regularnych. Ogólny opis algorytmu jest następujący. Stany konstruowanego automatu minimalnego etykietowane są zbiorami odpowiadającymi pewnym językom. Mając dany język L, ustanawiamy stan początkowy automatu jako L, wpisujemy go na listę i obliczamy a1L dla każdej litery aA. Jeśli wśród obliczonych wartości znajduje się język niewystępujący na liście, dodajemy go do listy. Obliczenia pochodnych Brzozowskiego wykonujemy, dopóki będziemy uzyskiwać nowe języki (nowe stany). Funkcja przejść konstruowanego automatu minimalnego zdefiniowana jest następująco:

f(X,a)=a1X,

gdzie

X

jest pewnym językiem z listy

.

Obliczone języki określają stany automatu minimalnego.

Uwaga [Dla zainteresowanych]

Obliczone języki określające stany automatu minimalnego to elementy monoidu syntaktycznego języka L.

Automatem minimalnym dla automatu 𝒜=(S,A,f,s0,T) będzie zatem automat

𝒜L=(SL,A,fL,s0L,TL),

gdzie:

  • SL={u1L: uA*},
  • s0L=L,
  • TL={u1L: uL},
  • fL(u1L,a)=a1(u1L)=(ua)1L.

Jeśli zdefiniujmy odwzorowanie ν:SSL, kładąc:

ν(s)=u1L, gdzie s=f(s0,u),

to można dowieść, że ν jest dobrze określone, jest epimorfizmem oraz ν(s0)=s0L - porównaj twierdzenie 3.1 z wykładu 4. Prawdą jest też, iż sT wtedy i tylko wtedy, gdy ν(s)TL oraz że następujący diagram komutuje:

tutaj diagram

Formalny zapis algorytmu przedstawiony jest poniżej.

Algorytm Minimalizuj1 - algorytm minimalizacji wykorzystujący pochodne Brzozowskiego


  1  Wejście: 𝒜=(S,A,f,s0,T) - automat taki, że L=L(𝒜)
  2  Wyjście: automat minimalny 𝒜=(SL,A,fL,sL,TL) dla 𝒜
  3  SL{L};
  4  włóż(,L);
  5  while  do 
  6    M zdejmij();   
  7    for each aA do 
  8      Na1M;
  9      if NSL= then
 10       SLSL{N};
 11       włóż(,N);
 12     end if
 13    end for
 14  end  while  
 15  for each MSL do
 16    for each aA do
 17      fL(M,a)a1M;
 18    end for
 19  end for
 20  sLL;
 21  TL{u1L: uL};
 22  return 𝒜;


Funkcja zdejmij(), występująca w linii 6., zdejmuje z kolejki pierwszy element i zwraca go jako swoją wartość. Procedura włóż(,L), występująca w liniach 4. oraz 11., wstawia na koniec kolejki element L.

Przykład 1.2.

Dla języka L=a+b+ z przykładu 1.1 (patrz przykład 1.1.) w wyniku działania powyższego algorytmu otrzymamy czterostanowy automat 𝒜L=(SL,A,f,L,T), gdzie
SL={L,a*b+,,b*},

a funkcja przejść zadana jest grafem:

<flash>file=ja-lekcja05-w-rys1.swf|width=250|height=250</flash>

<div.thumbcaption>Rysunek 1

Prezentowane poniżej twierdzenie uzasadnia kolejny algorytm konstrukcji automatu minimalnego. Tym razem punktem wyjścia jest dowolny automat rozpoznający język L.

Uwaga [Dla zainteresowanych]

Algorytm oblicza również monoid syntaktyczny języka L.

Analogicznie do konstrukcji relacji ρi, przedstawionej w wykładzie 4, możemy określić ciąg odpowiednich relacji na zbiorze stanów dowolnego automatu rozpoznającego język L. Relacje te służą do efektywnego określenia automatu minimalnego, równoważnego zadanemu.

Twierdzenie 1.1.

Niech 𝒜=(S,A,f,s0,T) będzie dowolnym automatem i niech L=L(𝒜). Przez 𝒜S×S oznaczmy relację równoważności zawierającą dwie klasy równoważności T i ST. Przez ρi dla i oznaczmy zstępujący ciąg relacji określony następująco:

ρ1=𝒜, a dla i=2,... przyjmijmy

ρi={(s,t)S×S:aA{1}f(s,a)ρi1f(t,a)}..

Wtedy ρi jest największą prawą kongruencją automatową zawartą w relacji ρ1=𝒜 i automat minimalny ma postać

𝒜L=(S/ρi,A,f*,[s0],T/ρi).

Dowód tego twierdzenia przebiega analogicznie do dowodu twierdzenia 3.3 z wykładu 4 (patrz twierdzenie 3.3. wykład 4).

Algorytm działajacy w oparciu o powyższe twierdzenie na podstawie zadanego automatu A=(S,A,f,s0,T), konstruuje efektywnie automat minimalny dla 𝒜, obliczając ciąg relacji ρ1. Proces konstrukcji przebiega w sposób iteracyjny. Zaczynamy od zdefiniowania relacji Parser nie mógł rozpoznać (nieznana funkcja „\mathacal”): {\displaystyle \approx _\mathacal {A}} , która posiada tylko dwie klasy abstrakcji: do pierwszej z nich należą wszystkie stany końcowe, a do drugiej - wszystkie pozostałe stany. Tak więc

s,tS s𝒜t(sTtT)(sSTtST).

Definiujemy pierwszą z relacji ρ, czyli relację ρ1 jako równą 𝒜, a następnie, dla każdej obliczonej już relacji ρi1, obliczamy relację ρi w następujący sposób:

s1ρis2(s1ρi1s2)(aA f(s1,a)ρi1f(s2,a)).

Nowo obliczona relacja ρi albo podzieli jedną lub kilka klas abstrakcji relacji ρi1, albo będzie identyczna z relacją ρi1. Jeśli zajdzie ta druga sytuacja, to znaczy, że dla każdego j>i w oczywisty sposób spełniona jest równość ρj=ρi, czyli ciąg relacji ustabilizuje się. W tym momencie algorytm kończy swoje działanie i klasy abstrakcji relacji ρi będą reprezentować stany automatu minimalnego.

Algorytm Minimalizuj2 - algorytm minimalizacji automatu wykorzystujący stabilizujący się ciąg relacji



  1  Wejście: 𝒜=(S,A,f,s0,T) - automat taki, że L=L(𝒜).
  2  Wyjście: automat minimalny 𝒜=(S,A,f,s0,T) dla 𝒜.
  3  ρ1𝒜;
  4  i1;
  5  repeat
  6    Parser nie mógł rozpoznać (nieznana funkcja „\slash”): {\displaystyle \displaystyle \slash \slash}
 oblicz ρi: s1ρis2(s1ρi1s2)(aA f(s1,a)ρi1f(s2,a));
  7    ii+1;
  8    empty(ρi) 
  9    for each (s1,s2)S×S do
 10      flagtrue;
 11      for each aA 
 12        if not f(s1,a)ρi1f(s2,a) then
 13          flagfalse;
 14        end if
 15      end for
 16      if flag=true and s1ρi1s2 then
 17        ρiρi{(s1,s2)};
 18      end if
 19    end for
 20  until ρi=ρi1
 21  Parser nie mógł rozpoznać (nieznana funkcja „\slash”): {\displaystyle \displaystyle S' \leftarrow S \slash \overline{\rho}_i}
;
 22  for each Parser nie mógł rozpoznać (nieznana funkcja „\slash”): {\displaystyle \displaystyle [s]_{\overline{\rho}_i} \in S \slash \overline{\rho}_i}
 do
 23    for each aA do
 24      f([s]ρi,a)[f(s,a)]ρi;
 25    end for
 26  end for
 27  s0[s0]ρi;
 28  T{[t]ρi: tT};
 29  return 𝒜=(S,A,f,s0,T);



Przykład 1.3.

Zminimalizujemy automat 𝒜=(S,A,f,s0,T), dla którego
S={s0,s1,s2,s3,s4,s5},A={a,b},T={s1,s2,s4} ,

a funkcja przejść określona jest przy pomocy grafu.

<flash>file=ja-lekcja05-w-rys2.swf|width=250|height=250</flash>

<div.thumbcaption>Rysunek 2

Konstruujemy ciąg relacji ρi.

Na początku ρ1 dzieli S na dwie klasy abstrakcji; pierwsza zawiera stany końcowe, a druga -- wszystkie pozostałe, czyli uzyskujemy dwa zbiory {s1,s2,s4} oraz {s0,s3,s5}.

Obliczmy ρ2 (pierwszy przebieg pętli w liniach 5.-20. algorytmu). Aby dwa elementy (stany) s i t były ze sobą w relacji ρ2 muszą być ze sobą w relacji ρ1 oraz musi zachodzić

aAf(s,a)ρ1f(t,a).

Czyli kolejna, nowa relacja może ewentualnie podzielić już istniejące zbiory zdefiniowane przez poprzednią relację. Nie może więc zajść taka sytuacja, że w jednej klasie abstrakcji relacji ρi+1 znajdą się elementy z różnych klas abstrakcji relacji ρi.

Rozważmy najpierw zbiór {s1,s2,s4}. Oczywiście każde dwa stany z tego zbioru są ze sobą w relacji ρ1. Zauważmy, że f(s2,a)=f(s4,a)=s2, f(s2,b)=f(s4,b)=s4, więc s2ρ2s4 (bo s2ρ1s2 oraz s4ρ1s4). Ponieważ f(s1,b)=s5 i f(s2,b)=s4, a wiemy, że s5 nie jest w relacji ρ1 z s4, zatem stany s1 i s2 nie mogą być ze sobą w relacji ρ2, a to oznacza, że także stany s1 i s4 nie mogą być ze sobą w relacji ρ2.

W analogiczny sposób można sprawdzić, że relacja ρ2 nie podzieli zbioru {s0,s3,s5}. Ostatecznie, po pierwszym wykonaniu pętli algorytmu minimalizacji obliczyliśmy relację ρ2, która dzieli S na następujące podzbiory:

{s1},{s2,s4},{s0,s3,s5}.

W kolejnym kroku obliczamy ρ3. Zbiór {s1} oczywiście nie może być już podzielony na mniejsze podzbiory. Łatwo zauważyć, że ρ3 nie podzieli także zbioru {s2,s4}.

Rozważmy teraz zbiór {s0,s3,s5}. Mamy f(s3,a)=f(s5,a) oraz f(s3,b)=s3, f(s5,b)=s5 i wiadomo, że s3ρ2s5, zatem s3 i s5 będą ze sobą w relacji ρ3.

Ponieważ f(s0,a)=s2 i f(s3,a)=s1, ale s2 i s1 nie są ze sobą w relacji ρ2, zatem nie mogą być także ze sobą w relacji ρ3. Relacja ρ3 dzieli więc zbiór {s0,s3,s5} na zbiory {s0} oraz {s3,s5}.

Podział zbioru S przez relację ρ3 wygląda więc

następująco:

{s0},{s1},{s2,s4},{s3,s5}.

Relacja ρ4 nie podzieli już ani zbioru {s2,s4}, ani zbioru {s3,s5}, więc uzyskujemy równość ρ4=ρ3 i ponieważ ciąg relacji się ustabilizował, algorytm kończy działanie.

Podsumowując, mamy:

ρ1
{s1,s2,s4},{s0,s3,s5},
ρ2
{s1},{s2,s4},{s0,s3,s5},
ρ3
{s1},{s2,s4},{s0},{s3,s5}.ρ3=ρ4 i równoważny minimalny automat 𝒜L=(S,f*,s0,T) ma 4 stany.

q0={s0},q1={s1},q2={s2,s4},q3={s3,s5}, T={q1,q2}. Jak łatwo zauważyć jest to automat z przykładu 3.1 zamieszczonego w wykładzie 4 (patrz przykład 3.1. wykład 4). }}

<flash>file=ja-lekcja05-w-rys3.swf|width=250|height=250</flash>

<div.thumbcaption>Rysunek 3

Jednym z najczęściej stosowanych algorytmów automatu minimalnego jest algorytm, który buduje "tabelkę" na podstawie której określa się automat minimalny. Poprawność tego algorytmu również uzasadnia twierdzenie 1.1 (patrz twierdzenie 1.1.).

W algorytmie tym wyznaczać będziemy tzw. stany rozróżnialne. Algorytm działa w czasie O(|A|n2), gdzie |A| jest mocą alfabetu, a n -- liczbą stanów automatu wejściowego, czyli podlegajacego minimalizacji. Złożoność pamięciowa jest również O(|A|n2). Prezentowany algorytm nosi nazwę algorytmu Hopcrofta-Ullmana. Znana w literaturze jest pewna zmodyfikowana wersja tego algorytmu. Jest to algorytm Aho-Sethiego-Ullmana, który ma tę samą złożoność czasową, ale lepszą złożoność pamięciową - O(|A|n). Natomiast w ramach ćwiczeń prezentujemy jeszcze jeden algorytm, znany jako algorytm minimalizacji Hopcrofta. Czas działania tego algorytmu wynosi O(nlogn).

Niech będzie relacją zdefiniowaną przez funkcję przejść automatu w następujący sposób:

pqwA* (f(p,w)Tf(q,w)T).

Definicja 1.2.

Stany p i q są równoważne, jeśli pq.

Jeśli stany nie są równoważne, to będziemy mówić, że są rozróżnialne.

Zadaniem algorytmu jest wyznaczenie stanów równoważnych, celem ich utożsamienia ze sobą. Algorytm musi zdecydować, dla każdej pary stanów (p,q), czy są one rozróżnialne. Jeśli pod koniec działania algorytmu okaże się, że nie stwierdziliśmy rozróżnialności tych stanów, to znaczy, że są one równoważne; następuje ich utożsamienie, czyli "połączenie" ich w jeden stan. Gdy takiego połączenia dokonamy dla wszystkich par stanów, wobec których nie stwierdziliśmy ich rozróżnialności, powstanie automat o minimalnej liczbie stanów.

W praktyce algorytm nie wyznacza stanów równoważnych, ale stany rozróżnialne, gdyż jest to po prostu łatwiejsze. Po wyznaczeniu wszystkich par stanów rozróżnialnych pozostałe pary stanowić będą stany równoważne.

W algorytmie wykorzystywać będziemy tablicę list [p,q], po jednej liście dla każdej pary stanów. Funkcja initialize ([p,q]) inicjalizuje listę pustą, funkcja zdejmij ([p,q]) zdejmuje jeden z elementów, natomiast funkcja włóż ([p,q],x) wkłada element x na listę [p,q]. Funkcja empty ([p,q]) zwraca wartość true gdy lista jest pusta, oraz false w przeciwnym wypadku. Zwróćmy uwagę, że elementami każdej z list [p,q] są pary stanów (s,t)S×S.

Algorytm Minimalizuj3 - algorytm minimalizacji wykorzystujący relację



  1  Wejście: 𝒜=(S,A,f,s0,T) - automat
  2  Wyjście: 𝒜=(S,A,f,s0,T) - automat minimalny taki, że L(𝒜)=L(𝒜).
  3  for each pT do
  4    for each qST do
  5      zaznaczone[p,q]1;
  6      initialize([p,q])
  7    end for
  8  end for
  9  for each (p,q)(T×T)((ST)×(ST)) do
 10    zaznaczone[p,q]0;
 11  end for
 12  for each (p,q)(T×T)((ST)×(ST)) do
 13    flagfalse
 14    for each  aA do
 15      if  zaznaczone[f(p,a),f(q,a)]=1 then
 16        flagtrue;
 17      end if
 18    end for
 19    if flag=true then
 20      OZNACZ(p,q);  para (f(p,a),f(q,a)) była oznaczona dla pewnego a;
 21    else
 22      for each  aA do
 23        if f(p,a)f(q,a) then
 24          włóż([p,q],(f(p,a),f(q,a)));
 25        end if
 26      end for
 27    end if
 28  end for
 29  Parser nie mógł rozpoznać (nieznana funkcja „\slash”): {\displaystyle \displaystyle S' \leftarrow S \slash_\equiv}
;  relacja  jest dopełnieniem tabeli zaznaczone
 30  for each Parser nie mógł rozpoznać (nieznana funkcja „\slash”): {\displaystyle \displaystyle [s]_{ \equiv} \in S \slash_\equiv}
 do
 31    for each aA do
 32      f([s],a)[f(s,a)];
 33    end for
 34  end for
 35  s0[s0];
 36  T{[t]: tT};
 37  return 𝒜=(S,A,f,s0,T);


Występujaca w algorytmie procedura Oznacz opisana jest poniżej.

Algorytm



 1  procedure OZNACZ(p,qS)
 2  if zaznaczone[p,q]=1 
 3    return
 4  end if
 5  zaznaczone[p,q]1
 6  while not empty([p,q]) do
 7    OZNACZ(zdejmij([p,q]))
 8  end while
 9  end procedure



Działanie algorytmu łatwo przedstawić na tabelce, która złożona jest z kwadratów - pól, odpowiadających parom stanów automatu. Fakt znalezienia przez algorytm pary stanów rozróżnialnych zaznaczamy symbolem "x" w polu tabelki odpowiadającym tej parze, co wykorzystamy w przykładzie.

Przykład 1.4.

Zminimalizujemy automat przedstawiony na rysunku 4, używając algorytmu Minimalizuj3.

<flash>file=ja-lekcja05-w-rys4.swf|width=250|height=250</flash>

<div.thumbcaption>Rysunek 4

Proces działania algorytmu i konstrukcji tabelki przedstawiony jest na poniższej animacji 1.

TUTAJ ANIMACJA. ja-lekcja5-w-anim1.

Wypełniona tabelka po zakończeniu działania algorytmu przedstawiona jest na rysunku 5.

<flash>file=ja-lekcja05-w-rys5.swf|width=250|height=250</flash>

<div.thumbcaption>Rysunek 5

Z tabelki odczytujemy, że stanami równoważnymi są stany s1,s5, stany s2,s8 oraz stany s4,s6. Automat minimalny przedstawiony jest na rysunku 6.

<flash>file=ja-lekcja05-w-rys6.swf|width=250|height=250</flash>

<div.thumbcaption>Rysunek 6