このページは[url=https://www.geogebra.org/m/ddm5798j]マス旅[/url]の一部です。[br][br]状態遷移図、オートマトン。[br]これらは、いろんな場面で聞くコトバですね。[br]イベントや条件判断によって、場合分けして状態を移り進む[br]仕組みや図というイメージです。（正確な定義はあとに回します。）[br][br]よくあるたとえが「自動販売機の反応」、「パスカルの三角形を剰余で色分けする」などですね。[br][url=https://www.geogebra.org/m/twxxx3yq#material/bykfdpmc]シェルピンスキー[/url]でセルオートマトンをやりました。忘れている人はざっと目を通しておきましょう。[br][br]さて、前回はクリーネ代数から正規表現が始まったことをやり、[br]合成数判定で、正規表現を使うと[br]まるで「[b]再帰的と言えるような検索[/b]」ができる体験をしましたね。[br]「[b]正規表現は、ツリーのトラバースや状態遷移図につながる[/b]」と言われてもピンとくるでしょう。[br][br]そこのところを、さらに調べてみよう。

[br][b][size=150]＜FA＞[br][/size][/b]オートマトンAは入力を1文字ずつ読みながら状態を移り、最終的に「受理」か「拒否」を判定するものです。[br]有限オートマトンFAは「状態を覚えながら1文字ずつ読む判定機」で、５つ組（状態集合Q、入力記号Σ、遷移関数δ、初期状態q0、受理状態F）を持ちます。[br]q0,Fはもちろん、Qの要素だから、特に最初の3要素がベースになりますね。[br]状態遷移図の読み方は省略します。[br]くわしくはこちら（[url=https://www.geogebra.org/m/ddm5798j#material/q6sj8zvg]状態発生をモデル化[/url]）の「状態遷移図」[br][br]また、FAは決定性DFAと非決定性NFAに2分されます。[br][br][b][size=150]＜DFA＞[br][/size][/b][br]DFAのD（Deterministic決定性）は何が決まっているということでこの名がついているかです。[br]それぞれの状態Qでそれぞれの入力記号Σに対して、どの状態Qに移るかがはっきり一意的に決まるというのが決定性があるということです。[br]つまり、δ：Q×Σ→Qが関数と言えるものです。[br]たとえば、[br]５つの組（Q,Σ、δ、q0、F）のうち、δ以外を決めてみる。[br]Σ＝{0,1},Q={q0:直前が１,q1:直前が0、q2:=[01]を見つけた=F}[br]q0で1を読むと、この1が「直前が1」の意味となるから、q0状態のままとなり、自己ループにはまる。[br]q0で0を読むと、この0が「直前が0」の意味となるから、q1状態に移る。q1にいて0を読んでも「直前が0」の意味となるから、q1から抜けられないが、q1にいて1を読むと[01]を見つけたから、q2に移動し、受理状態となる。q2に到達したあとは、次に0が来ても1が来ても[01]を含んでいることには変わりないので、q2に留まり続けます[br]このように、関数δ(qi,j)=qkで、i,jからkが１つに決まるので、これはDFAと言えるね。[br][br][b][size=150]＜NFA＞[br][/size][/b][br]NFAのN(Non)は一意的な状態への移動が決まるの否定です。[br]だから、１つの状態の否定だから、複数の状態や0個の状態も[br]ありという、あいまいなマシンです。[br]つまり、空集合も入れた状態Qのあらゆる部分へ移動できます。[br]言い換えると、Qの部分集合、P(Q)の中のどれか1つの集合に一斉に移る（状態のアンサンブル・並行世界）ということです。状態集合Qのサイズがnならば、2^n通りの可能性があるということですね。[br][br]ぜんぜん別世界に見えますが、当たるも八卦当たらぬも八卦、[br]いろいろやってみてうまくいったらラッキーという進み方ができるということですね。もちろん、行き止まりでアウトということもあります。[br]これはDFAには自己ループがあり受理できない場合があるので、NFAだからダメだということにはなりません。[br][br]前回の例でいうと、強欲マッチは１つの最長を決めますが、[br]無欲マッチはまったり粘り強く可能性を広げましたね。[br]この違いはDFAとNFAの違いに似てますね。

[br]NFAはゆるく無欲なので、書きやすそうですね。[br]DFAは明確なので受理判定が速そうですね。[br][br]そのいいとこどりができます。[br]人間がNFAで書いた表現をDFAに変換してからPCに実行してもらうのです。[br][br][size=150][b]＜べき集合構成＞[/b][/size][br][br][b]実際NFAはPC内部でDFAに変換されて実行されている[/b]そうです。[br][br]そこで活躍するのが「べき集合」です。[br]人間がNFAをプログラミングによってDFAに変換することもできるでしょう。[br][br]NFAの初期状態q0をそのままDFAの初期状態q0とします。[br]NFAの状態集合Qから作ったべき集合P(Q)＝Q'をDFAの状態集合とします。[br]新状態Q'のうちNFAの受理状態Fを含む状態すべてがDFAの受理状態集合F'となります。[br][b]NFAの「状態の組み合わせ（集合）」を、DFAの「1つの新しい状態」とみなす操作[/b]を[br]「[b]べき集合構成[/b]」と言います。[br][br][b][size=150]＜ムーアのアルゴリズム＞[br][/size][/b][br]実際は状態の数が膨大になるので、[br]受理群と非受理群に状態をグループ分けして、[br]文字を入れてそれぞの文字の行先が同じグループになるかどうかで、[br]グループを2分します。これを繰り返すと分割が起きなくなります。[br]この機能としての状態遷移図（遷移関数）の簡略化を図ることにより、[br][b]最小DFA[/b]が得られるそうです。[br][br][b][size=150]＜正規表現（有限オートマトン）の限界＞[br][/size][/b][br]a^n b^n（aをn個並べた後にbを同じだけ並べる）は正規表現ではできないのです。[br][br]「[b][color=#0000ff]正規表現でHTMLパースはするな[/color][/b]」[br][br]という格言につながりますね。[br][br]コンパイラの構文解析やHTMLのタグ対応にはプッシュダウン・オートマトン（スタック）が必要ですが、TCPプロトコル解析や単純なデータ切り出しなど、[b]入れ子構造がない場面[/b]では[br][b][size=200][color=#0000ff]正規表現（有限オートマトン）[br][/color][/size][/b]が実質的な標準として大活躍しています。

[br][color=#9900ff][b][size=150]課題：import re に頼らず、状態遷移表で「aが偶数個あるか判定するDFA」などを再現してみよう。文字集合Σ＝{a,b,c,d}とします。[br][/size][/b][/color][br][b][size=150]＜Python＞[br][/size][/b][br]#aが偶数個あると受理するDFA[br]#Σ={a,b,c,d}[br]#q_even(偶数個状態)=F(受理）,q_odd(奇数個状態),q0(開始)=q_even扱い[br]#qcur(現在),qnext(次)[br][br][b]def transit(cur, char):[/b][br] # アルファベットの定義[br] alphabet = {'a', 'b', 'c', 'd'}[br] [br] [b] # 定義外の文字が来たら、またはすでにエラー状態なら、エラー状態を維持[/b][br] if char not in alphabet or cur == 'q_err':[br] return 'q_err'[br] [br] [b]# aが来たら状態を反転[br][/b] if char == 'a':[br] if cur == 'q_even':[br] return 'q_odd'[br] if cur == 'q_odd':[br] return 'q_even'[br] [br] [b]# b, c, d の場合は現在の状態（q_even / q_odd）を維持[/b][br] return cur[br][br][b]def DFA(target):[br][/b] # 最初から「q_even（偶数個状態）」でスタートする[br] cur = 'q_even'[br] [br] for item in target:[br] cur = transit(cur, item)[br] [br] # 最終的に q_even の場合のみ受理（True）とする[br] return cur == 'q_even'[br][br][color=#0000ff]# --- テスト実行 ---[br]test_cases = [[br] "", # aは0個 (偶数) -> True[br] "bcd", # aは0個 (偶数) -> True[br] "ab", # aは1個 (奇数) -> False[br] "aba", # aは2個 (偶数) -> True[br] "aacaada", # aは4個 (偶数) -> True[br] "abcdaaaa", # aは5個 (奇数) -> False[br] "abx" # 定義外の文字 'x' があるため -> False[br]][br][br]for test in test_cases:[br] result = DFA(test)[br] print(f"文字列: '{test:<8}' -> 受理結果: {result}")[br][/color][size=150][color=#0000ff][b][br][/b][/color][b][color=#0000ff][size=100]では、if文によってδ関数を代用してましたが、[br]δ関数を辞書化によって、見やすくするとどうなるでしょうか。[br][/size][/color][br]＜辞書でδ関数を実装する＞[/b][/size][br][br]# --- 辞書によるDFAの定義 ---[br][br]# 状態集合 Q[br]states = {'q_even', 'q_odd', 'q_err'}[br][br]# アルファベット Σ[br]alphabet = {'a', 'b', 'c', 'd'}[br][br]# 初期状態 q0[br]start_state = 'q_even'[br][br]# 受理状態集合 F[br]accept_states = {'q_even'}[br][br][color=#0000ff][b]# 遷移関数 δ: (状態, 文字) -> 次の状態[/b][/color][br][b]# ※ 定義外の文字やエラー状態からの遷移は get() のデフォルト値で処理[br]delta = {[br] ('q_even', 'a'): 'q_odd',[br] ('q_even', 'b'): 'q_even',[br] ('q_even', 'c'): 'q_even',[br] ('q_even', 'd'): 'q_even',[br][br] ('q_odd', 'a'): 'q_even',[br] ('q_odd', 'b'): 'q_odd',[br] ('q_odd', 'c'): 'q_odd',[br] ('q_odd', 'd'): 'q_odd',[br]}[br][/b][br]def dfa_dict(target_string):[br] current_state = start_state[br] [br] for char in target_string:[br] # 辞書から次の状態を取得（未定義の文字が来たら 'q_err' へ移行）[br] current_state = delta.get((current_state, char), 'q_err')[br] [br] # エラー状態に入ったらその時点で処理を中断（トラップ状態）[br] if current_state == 'q_err':[br] return False[br][br] # 最終状態が受理状態集合 F に含まれるか判定[br] return current_state in accept_states[br][br][br]# --- テスト実行 ---[br]test_cases = [[br] "", # aは0個 (偶数) -> True[br] "bcd", # aは0個 (偶数) -> True[br] "ab", # aは1個 (奇数) -> False[br] "aba", # aは2個 (偶数) -> True[br] "aacaada", # aは4個 (偶数) -> True[br] "abcdaaaa", # aは5個 (奇数) -> False[br] "abx" # 定義外の文字 'x' があるため -> False[br]][br][br]for test in test_cases:[br] result = dfa_dict(test)[br] print(f"文字列: '{test:<8}' -> 受理結果: {result}")[br][br][b][size=150]＜振り返り＞[br][/size][/b][br]課題では正規表現を使わずにということですが、正規表現＝有限オートマトンだから、[br]オートマトンでかけるならば、正規表現でもコードが作れますね。[br][br]import re[br][br]# 正規表現パターンの定義[br]# ^ : 先頭[br]# [bcd]* : 最初の 'a' が出るまでの 'b','c','d' (0個以上)[br]# (?: ... )* : 非キャプチャグループの繰り返し[br]# a[bcd]* : 1個目の 'a' と、その後の 'b','c','d'[br]# a[bcd]* : 2個目の 'a' と、その後の 'b','c','d' (これで 'a' が偶数個ペアになる)[br]# $ : 末尾[br]pattern = re.compile(r"^[bcd]*(?:a[bcd]*a[bcd]*)*$")[br][br]def dfa_re(target_string):[br] # パターンに完全マッチするかどうかを判定[br] return bool(pattern.match(target_string))[br][br][br]# --- テスト実行 ---[br]test_cases = [[br] "", # aは0個 (偶数) -> True[br] "bcd", # aは0個 (偶数) -> True[br] "ab", # aは1個 (奇数) -> False[br] "aba", # aは2個 (偶数) -> True[br] "aacaada", # aは4個 (偶数) -> True[br] "abcdaaaa", # aは5個 (奇数) -> False[br] "abx" # 定義外の文字 'x' があるため -> False[br]][br][br]for test in test_cases:[br] result = dfa_re(test)[br] print(f"文字列: '{test:<8}' -> 受理結果: {result}")[br][br]

[b][size=150]＜geogebra＞[br][/size][/b][br][b][color=#0000ff][code]"a"[/code] → [code]-1[/code]（偶数・奇数の状態反転：1 ⇄ -1）[br][code]"b", "c", "d"[/code] → [code]1[/code]（状態維持：1 × 1 = 1, -1 × 1 = -1）[br][code]その他[/code] → [code]0[/code]（エラー状態への遷移：何をやっても 0）[br][/color][/b][br]このように文字を数値に置き換えることで、文字読み取りが代数演算に置き換わります。[br]さらに、演算結果を状態を表す点の位置として表示します。[br]この「代数的な演算」がそのまま「オートマトンの状態遷移」になりますよ。[br][br]「0をかけたら二度と抜け出せない（トラップ状態）」[br]というDFAのエラー特性を、数値の積 [code]Product(k)[/code] で表現します。[br][br]安定して動かすために、動いた結果をTakeでコマ切れにして作り、[br]それをProductによって、累積効果を先に計算しておきます。[br][br]スライダーによって読み取り位置を後付けで追いかけると、[br]スライダーで1文字ずつ入力を読み取っているような振る舞いが[br]表示でますね。[br][br][b]#数式ビューに次の内容を1行ずつコピペします。#は設定[br]S_e=(1,0) # 偶数状態（緑）[br]S_o=(-1,0) # 奇数状態（黄土）[br]Serr=(0,0) # エラー状態（赤）[br][/b][br]InputChar="abcdefa"[br]InputBox(InputChar)[br]Le=Length(InputChar)[br]digit=Sequence(If(Element(InputChar, k) ≟ "a", -1, If(Element(InputChar, k) ≟ "b" ∨ Element(InputChar, k) ≟ "c" ∨ Element(InputChar, k) ≟ "d", 1, 0)), k, 1, Le)[br][br]sub=Sequence(Take(digit,1,k),k,1,Le)[br]Ps=Zip(Product(k),k,sub)[br]P0=Join({{1},Ps})[br]Num=Slider(1,Le+1,1) #設定＞上級＞アニメーション＞モード＞増加（1回）[br]pos=Element(P0,Num)[br]dispChar="_"+InputChar[br][br][b]P=(pos,0) # 現在の状態を示す動点（青・サイズ最大）[br][/b][br]poschar=Element(dispChar,Num)[br]a=Num-1[br]text=""+a+"文字目"+poschar+""[br][br][b]#操作の仕方[br][/b][br]0文字目は開始前の状態を表します。[br]入力ボックス1の"abcdefaa"を適当に変えましょう。[br]Numを手動で１に戻します。[br]Numのプレイボタンを押すと、オートマトンに1文字ずつ[br]データが入り、Pが今の状態になります。[br]S_eがaが偶数個状態[br]S_oがaが奇数個状態[br]Serrが｛a,b,c,d｝以外の文字の入力によるエラー[br]Numがとまると判定がPの位置でわかります。[br]PがS_eで受理されますね。[br][br]タイトルは「a,b,c,dの文字列でaが偶数なら受理するオートマトン」[br][br]一度エラーになると、aが偶数個でもエラーになりますね。[br]これは、0を入れた積が0になることが効いてます。[br]