このワークシートは[url=https://www.geogebra.org/m/twxxx3yq]Math by Code[/url]の一部です。[br][br][br]
[b]探検というコトバ[/b]からどんなようすをイメージしますか？[br][br]幅広く舗装されたまっすぐの太い通りをひたすら直進する。[br]これを探検といいませんよね。[br][br]探検という感じがするのは、道が曲がりくねっていて、[br]先が見えない、樹木やぬかるみが邪魔をする。[br]行き止まりがあるかもしれない。[br]まさに迷路の探索だ。[br][br]でも、探索というのは冒険家の専売特許でないです。[br][br]みなさんがふだん使っているPCやスマホでも探索はある。[br]そう、[b]データ探索[/b]だね。[br][br]今回は身近なデータ探索につながる３つのパタン、[br][b]イテレタ、インタプリタ、ビジター[/b][br]題して、[b]探索３兄弟[/b]をやろう。
＜イテレタ＞[br][br]イテレタiteratorは[b]反復iterするもの反復子[/b]のことです。[br][br]GOF本では、[br][b][color=#0000ff][size=150]イテレタ：[br]集約オブジェクトの内部が公開されてなくても、その要素に順にアクセスすることを助ける。[br][/size][/color][/b][br]たとえば、本棚の並ぶ本の値段はカバーの後ろにかいてありますが、[br]自分の持ってる本で、一番定価が高い本を探したいとしよう。[br]本棚に本が何冊あるかは別として、[br]一番上の段からみる、１つの段は左から順にみる。[br]一番下の段の一番右の本が最後の本になる。[br][br]本棚の本たちをbooksとしよう。[br]そこで、[br]python なら、[br]for book in books:[br] print(book.price)[br][br]この順番に見る働きをクラスに格上げしたものが[b]イテレタモジュール[/b]です。[br]イテレタを使うには、反復できるオブジェクトを用意します。[br][b]イテラブル（順に進める）なアグリゲート（集約オブジェクト）[/b]です。[br]順不同でよい集まり例えば、マップ、辞書、セット（集合）はイテラブルではありません。[br]イテラブルはnext（次）をつかって、順番に調べるのですね。[br]つまり、[b][color=#0000ff]最初から最後まで順番がつく集まりを順に調べるクラスがイテラブル[/color][/b]です。[br][br]なんて、[b]地味[/b]なんでしょうね。[br][br]なんで、こんな[b]単純なパタン[/b]に名前をつけたのでしょうか。[br][br]それは、順番をつけて全件調べるという[b]地味な調査が一番確実[/b]でよく使われるからです。[br]さらに大切なのは、[b][color=#0000ff][size=150]インターフェース、共通キーワードを決めておくと、[br]調べる対象が変化しても、調べる対象の実装が混ざっていても平気になるから[/size][/color][/b]です。[br]イテラブルな調査対象は、[b]リストと配列とタプルくらい[/b]ですが、[br]それが変更になっても、まざっていも大丈夫。[br]ユーザーは同じキーワードで調査をイテラブルさんに対してお願いできちゃうのですから。[br][br]それでは、実装してみよう。[br]昔、「[b]世界のナベアツ[/b]」というのがはやりました。[br][b]３の倍数か３のつく数のときに変顔をする[/b]という芸人です。[br]ふつうに考えると、１から順に整数Nを表示していくのですが、[br]整数Nが「３の倍数」（A条件）か「３のつく整数」（B条件）のときに、「変顔」を表示するので、[br]NがAまたはBであることを[b]if文で作り込むのが普通[/b]のプログラミングでしょうね。[br][br]あえてこれを、イテレタパタン、[b]順次検索共通仕様[/b]でやってみましょう。[br]ナベアツイテレタというインターフェースを作ります。[br][b]調査の最大数100[/b]のコンストラクターを持ちます。[br][b]nabe(x)というメソッド[/b]を持ちます。[br]この子クラスが３の倍数リストイテレタ[b]multi3[/b]と３のつく整数リストイテレタ[b]hasDidit3[/b]です。[br]ともにnabe(x)というメソッドをもち、[b]リストに属するときはtrueを返します[/b]。[br]だから、[br]M= multi3()[br]D=hasdDigit3()[br]でインスタンスを作れば、[br]ユーザーは[br]for i in range(1,100+1):[br] if M.nabe(i) or D.nabe(i):[br] print("変顔”)[br] else:[br] print(i)[br]とするだけです。[br]multi3とhasDigit3の中身を知らなくてもよいのです。[br]たとえば、multi3は0,1,2,3,4,5,.....のインデックス順に{true, false, false, true, false,false,....}というリストを[br]作っているかもしれません。それとも単純に３の倍数のリストにしているかもしれません。[br]一方のhasDigit3は[br]A={ 30+x for x in range(10)} , B={ 3+ 10 x for x in range(10)},C={33}[br] D=list(（A｜B） - C)のように、集合演算で作っているかもしれません。[br][br] このようにリストの作り方は違いますが、[br]nabe(n)の実装をしておけば、ユーザーは同じ命令が使えます。[br]さあ、実装してみよう。
[size=150]＜[b]順次検索共通仕様[/b]＞[/size][br]#nabeatsu.py[br][color=#0000ff]# ==========================================[br]# 【インターフェース】ナベアツイテレータの共通契約[br]# ==========================================[br][/color]class NabeatsuIterator:[br] """内部のリスト構造がどうであれ、共通のキーワード『nabe(x)』を保証する"""[br] def __init__(self, max_num=100):[br] self.max_num = max_num[br] self.target_list = self._create_list() # 子クラスごとに異なるリストが作られる[br][br] def _create_list(self):[br] """【型枠】リストの具体的な作り方は子クラスに委ねる"""[br] raise NotImplementedError()[br][br] def nabe(self, x):[br] """ユーザーが使う共通キーワード。リストに属していればTrueを返す"""[br] return [b]x in self.target_list[/b][br][br][br][color=#0000ff]# ==========================================[br]# 【子クラス1】3の倍数リストイテレータ（数理的なアプローチ）[br]# ==========================================[br][/color]class Multi3Iterator(NabeatsuIterator):[br] def _create_list(self):[br] # 1からmax_numまでの間で、3の倍数だけを格納したリストを作る[br][b] return [i for i in range(1, self.max_num + 1) if i % 3 == 0][br][/b][br][color=#0000ff]# ==========================================[br]# 【子クラス2】3のつく整数リストイテレータ（集合を駆使したアプローチ）[br]# ==========================================[br][/color]class HasDigit3Iterator(NabeatsuIterator):[br] def _create_list(self):[br] # 10の位が3の集合[br] A = {30 + x for x in range(10)}[br] # 1の位が3の集合[br] B = {3 + 10 * x for x in range(10)}[br] # 重複する33の集合[br] C = {33}[br] [br] # 集合の和（OR）から重複を引き、1からmax_numに収まるものだけをリスト化[br][b] union_set = (A | B) - C[br] return [i for i in range(1, self.max_num + 1) if i in union_set][br][/b][br][br][color=#0000ff]# ==========================================[br]# 【実行・ユーザー側のコード】中身を知らなくても、同じ命令で全件調査！[br]# ==========================================[br][/color][br][b]# 1. 異なる内部実装を持つ、2つのイテレータインスタンスを生成[br]M = Multi3Iterator(100)[br]D = HasDigit3Iterator(100)[br][/b][br]print("ナベアツさん登場！ \n")[br][br]# 2. ユーザーは共通の「nabe(i)」だけで、1から100まで確実に走査する[br]for i in range(1, 100 + 1):[br][b][color=#0000ff] # Mの中身が倍数リストか、Dの中身が集合演算か、ユーザーは一切知らなくてよい[br][/color] if M.nabe(i) or D.nabe(i):[br] print("（変顔）", end=" ")[br] else:[br] print(f"{i}", end=" ")[br][/b] [br] # 10個ごとに改行して見やすく整える[br] if i % 10 == 0:[br] print()[br][OUT][br]1 2 （変顔） 4 5 （変顔） 7 8 （変顔） 10 [br]11 （変顔） （変顔） 14 （変顔） 16 17 （変顔） 19 20 [br]（変顔） 22 （変顔） （変顔） 25 26 （変顔） 28 29 （変顔） [br]（変顔） （変顔） （変顔） （変顔） （変顔） （変顔） （変顔） （変顔） （変顔） 40 [br]41 （変顔） （変顔） 44 （変顔） 46 47 （変顔） 49 50 [br]（変顔） 52 （変顔） （変顔） 55 56 （変顔） 58 59 （変顔） [br]61 62 （変顔） 64 65 （変顔） 67 68 （変顔） 70 [br]71 （変顔） （変顔） 74 （変顔） 76 77 （変顔） 79 80 [br]（変顔） 82 （変顔） （変顔） 85 86 （変顔） 88 89 （変顔） [br]91 92 （変顔） 94 95 （変顔） 97 98 （変顔） 100
[b][size=150]＜インタプリタは通訳一式＞[/size][/b][br][br]インタプリタというのは、その名の通り通訳ですね。[br]通訳するのは今はAIもやりますが、人間のイメージが強いですね。[br]同時通訳は聞いてるだけで大変そう。忙しそうだし。[br]でも、インタプリタは楽ですよ。[br]なにしろ、文法から何からはっきりしているからです。[br]たとえば、RPNという計算記法があります。[br]逆ポーランド記法の略で、[br]中身は演算を真ん中にかくふつうの記法（中置）A＋Bを「AB+の順にかく（後置）記法」のことです。[br]後置記法RPNの何がよいかというと、[br]中置だと、かっこのある式 (1 + 2) × 3[br]これを後置にすると、 1 2 + 3 ×[br]で、いいのです。[br]前とか中というコトバは数式ができあがってから人間が名付けた言い方です。[br]これをPCがわかるように伝えないと、プログラミングには使えません。[br]それがインタプリタです。[br]そう、PCにわかるように[b]文法と単語まで教えます[/b]。それが、[b]BNF[/b]情報です。[br]１けたの計算をするRPNをBNF（バッカス・ナウア式）で表します。[br]BNFの基本記号は３つです。[br][code]A::=[/code] B: 「AはBで定義される」[br][code]A|B[/code] : 「AまたはB」[br]＜A＞：Aは構成要素[br]RPNのルールは３つです。[br]＜digit＞ ::= 0 | 1 | 2 | ... | 9[br]＜operator＞ ::= + | - | * | /[br]＜expression＞ ::= ＜digit＞ | ＜expression＞ ＜expression＞ ＜operator＞[br]数式の定義が自己参照しているからこそ再帰的に読み取ることができるのですね。[br][br]さて、データ構造といえばスタックかキューかハッシュの出番かな？[br]左から順に、つまり入れた順に使うわけではないから、スタックでしょうか。[br]stackにSさんと名前をつけよう。[br]RPN＝「4 [u]3 2[u]＋[/u][/u] 1 ー ＊」を空白区切りでSさんにわたそう。[br]Sさんは数値ならどんどん上に積んでいく。[br]Sに左から4、3、2の順を渡すと、上か順に[u]ｘ、ｙ、残りをysと名付けるなら、[/u]x=２, y=３, ys=４となる。[br][color=#0000ff]もとの式の左右の順番とｘ，ｙの順番が逆になってしまうので計算順に注意しよう。[br][/color]次は（＋）が出たから、上からx=２,y=３を取り出して、y+x=３+２=５をSさんに戻す。[br]残った４の上に５が積まれるので、上から順にx=5,y=4となる。[br]次の1はSさんに渡す。Sさんは上から順にx=1,y=5,ys=4となる。[br]次は（ー）が出たから、上からx=1,y=5を取り出して、y-x=5-1=4を戻す。[br]Sさんは上から順にx=4,y=4となる。[br]次も（＊）が出たから、上からx=4, y=4を取り出し、y*x=4*4=16を戻す。[br]x=16が入っている。[br]RPNが空になったとき、Sさんのスタックが１つの数値だけになりそれが答えになるね。[br][br][b][size=150]＜RPNインタプリタをHaskellで作ってみよう。＞[/size][/b][br]RPN=文字列を受取り、和を返す関数calcDiff( RPN)を作ってみよう。[br]文字列RPNをリストLに分解するwords関数を作っておき、使おう。[br]Lを左から順に読み取って処理は、左畳み込みfoldl関数にfoldStack関数にLと空のリストを渡す。[br]最後に畳み込んだあとのリストの先頭をhead関数で読み込む。[br]この３つの関数をパイプラインでつないだ、合成関数がcalcDiffです。[br]ただし、foldStack関数だけは未定義なので、where句で末尾にケース分けして反応を定義します。[br]foldStack関数が使うスタック（だたのリストです）Sは最初は空です。S＝[ ][br]Lからの読み取りが数値numなら, Sの中身xsの先頭にnumを追加する。Sの先頭から順にx,y,ysとなる。[br]Lからの読み取りが”-”なら,Sの先頭からx,yなので、y-xを計算してSの先頭に追加する。[br]最後にSの先頭(head)を読み取る。[br]これをHaskellでかくとこうなります。[br][color=#0000ff]calcDiff = head . foldl foldStack [] . words [br] where [br] foldStack (x:y:ys) "-" =(y-x):ys[br] foldStack xs num = read num :xs[br][/color]ghci>calcDiff " 10 - 2 "[br]8[br]ghci>calcDiff " 2 - 10 "[br]-8[br]Haskellなら、読み取った演算記号が-以外の場合の式をwhere句に、好きなだけ追加するだけで[br]四則計算ができるように改造できます。[br]たとえば、sum や ln など、haskellが読める関数がそのままかけます。[br][color=#0000ff][size=150]calcRPN = head . foldl foldStack [] . words [br] where [br] foldStack (x:y:ys) "+" =(y+x):ys[br] foldStack (x:y:ys) "-" =(y-x):ys[br] foldStack (x:y:ys) "*" =(y*x):ys[br] foldStack (x:y:ys) "/" =(y/x):ys[br] foldStack (x:y:ys) "^" =(y**x):ys[br] foldStack (x:xs) "ln" = log x:xs[br] foldStack xs "sum" =[sum xs][br] foldStack xs num = read num :xs[br][/size][/color][b][br][color=#9900ff][u]課題：juliaでRPNインタプリタをつくろう。[/u][/color][/b][br]juliaにもfoldl,foldr関数があるのですが、簡単な演算のくり返しでは作動しますが、[br]一般の関数は難しいようです。仕方ないのでwhileで代用してます。[br]words関数とfoldStack関数を作っておき、単体で作動テストをしてから、[br]それらを組み込んだcalcRPN関数を動かしてます。[br][color=#0000ff]function words(x) # 1スペースで区切る[br] return split(x," ")[br]end[br]function foldStack(lst)[br] ope = ["*","+","-","/"][br] opeS = Set(ope)[br] Stack =[][br] item = popfirst!(lst)[br] # y x ope を抜き出して Stack の上からx ,y の順につまれた数で[br] # z= y ope x を計算します。[br] while item ∉ opeS[br] pushfirst!(Stack,item)[br] item = popfirst!(lst)[br] end[br] if item in ope[br] x = parse.(Int,popfirst!(Stack))[br] y = parse.(Int,popfirst!(Stack))[br] z = 0[br] if item == ope[1][br] z = y * x[br] elseif item == ope[2][br] z = y + x[br] elseif item == ope[3][br] z = y - x[br] elseif item == ope[4][br] z = y / x[br] end[br] end[br] #取り出した３要素y,x,opeを１要素zにしてリストlstに戻し、[br] #計算しなかったStackの残りも、リストlstにもどします。[br] pushfirst!(Stack,z)[br] while length(Stack)>0[br] p = string(popfirst!(Stack))[br] pushfirst!(lst,p)[br] end[br] return lst[br]end[br][/color][color=#0000ff]function calcRPN(RPN)[br] res = RPN |> words[br] while length(res)>2 [br] res = foldStack(res)[br] end[br] return parse.(Int,res)[br]end[br]"4 3 2 + 1 - *" |> calcRPN[br][/color]#==========================[br][OUT] [br]1-element Vector{Int64}:[br] 16[br]
[b][color=#9900ff][u][size=150]課題：これを「Python」で実装してみよう。[br][/size][/u][/color][/b][br]from functools import reduce[br][br][b]def calcDiff(RPN: str) -> int:[br][/b] """【Interpreter】[br] 逆ポーランド記法（RPN）の文字列を通訳し、数式の意味を解釈して計算結果を返す。[br] words ➔ foldl(foldStack) ➔ head のパイプライン合成。[br] """[br] [br] # 1. 【words関数】文字列をスペースで区切ってリスト（トークン列）に分解する[br][b] def words(text: str) -> list:[br][/b] return text.split()[br][br] # 2. 【foldStack関数】左畳み込み(foldl)の各ステップで、スタックに値を積む、または計算する[br] # Pythonの reduce(foldl) に合わせて、第1引数をスタック(S)、第2引数を読み取った文字(item)とする[br] def foldStack(S: list, item: str) -> list:[br] operators = ["+", "-", "*", "/"][br] [br] if item in operators:[br] # 演算子を読み込んだ場合：スタックの先頭からx, yを取り出す（Sの先頭から順にx, yになる）[br] x = S.pop(0)[br] y = S.pop(0)[br] [br] # 演算子に応じた計算を行い、結果zをスタックの先頭(xs)に戻す[br] if item == "+": z = y + x[br] elif item == "-": z = y - x[br] elif item == "*": z = y * x[br] elif item == "/": z = int(y / x) # 割り算は整数丸め[br] [br] S.insert(0, z)[br] else:[br] # 数値を読み込んだ場合：整数に変換してスタックの先頭(xs)に追加する[br] S.insert(0, int(item))[br] [br] return S[br][br] # 3. 【head関数】畳み込んだあとのスタックの先頭要素を読み込む[br][b] def head(S: list) -> int:[br][/b] return S[0][br][br] # ========================================================[br] # 【パイプライン合成の実行】[br] # 提示された数式「 1 2 + 3 * 」の文脈を、左から右へと畳み込んで解釈していく[br] # ========================================================[br] tokens = words(RPN) # ① words: 文字列をリストLに分解[br] final_stack = reduce(foldStack, tokens, []) # ② foldl: 空リスト[]からスタートしてトークンを畳み込む[br] return head(final_stack) # ③ head: 先頭の計算結果を回収[br][br][br]# ==========================================[br]# 【作動テスト】Juliaのサンプルと同じ数式で検証[br]# ==========================================[br]print("インタプリタ\n")[br]rpn_expr1 = "4 3 2 + 1 - *"[br]result1 = calcDiff(rpn_expr1)[br]print(f"RPN表記: {rpn_expr1}")[br]print(f"[OUT] ➔ {result1}")[br]
[b][size=150]＜振り返り＞[br][/size][/b][br]ビジターを3番目に扱う予定だったけれど、[br]イテレタ、インタプリタにかなりの行数をとってしまた。[br]発想は単純だけれども、[br]実装するとなると場所をとってしまう。[br]やはり、[br]地味に確実に動くものを作るということは、[br]いざやってみると細かな設定とぬけのない記述が必要になり、[br]順次的な処理では、クラス数は少ないのに[br]１つ１つのクラスの処理が複雑に流れがちだからか。[br][br]ということで、ビジターを使った探検は次回のお楽しみにしよう。[br][br]上にあげた、geogebraで作ったRPNインタプリタの作り方は[br][url=https://www.geogebra.org/m/twxxx3yq#material/fydswe7g]calcRPN関数で電卓を作ろう[/url]のページの後半を参照してください。