このワークシートは[url=https://www.geogebra.org/m/twxxx3yq]Math by Code[/url]の一部です。[br][br]これまで、FPシリーズとして、[br]OOPのデザインパタンの純粋関数型への移植から始めた。[br]そして、データを載せるコンテナとしてのリスト、辞書を扱ったね。[br]その次に道具としてのモジュール、operator,functools,itertoolsと進んできた。[br][br]今回は趣向を変えて、プログラミングの発想・記法である[b]再帰[/b]を扱おう。

再帰は発想と記法であり、ツールではない。[br][br]再帰recursionは自己参照で、反復itelateと対照的な意味で扱われる。[br]itertoolsは生真面目な顔をして堂々とさぼるのが上手なyield,genetatorを部品とするレイジー関数たちがバケツリレーで役立つ世界を作ってくれていた。[br][br]再帰は、道具というよりは、もっと広く、深い。[br][br][color=#0000ff][b][size=150][size=200]再帰は、発想であり、記述の特徴だ。[br][/size][/size][/b][/color][br]たとえば、reduceというドミノ倒し関数がfunctoolsにあったが、[br]累乗、階乗、累加という発想自体は再帰的と言える。[br][br]reduceでは初期値とラムダ式があったように、[br]端っこと変化方法のペアがないと[br]無限に続いてしまう。[br][br]再帰関数は終了条件baseと再帰式recursionがある。[br]やはり、端っこと変化方法がペアになっている。[br][br]しかし、再帰関数のよくある例は階乗、フィボナッチ数列くらいだ。[br]苦労して学んでも、例が乏しい本が多い。[br][br]再帰はふつうにやると、[br]スタックの山ができて、いつも１から計算するため計算の山ができる。[br]それを避けるためのわざは、[br]パッケージFunctoolsのlru_cacheやcashというデコレタを頭につけるメモ化があった。[br]他にもforループ回しに直す方法があるが、スタックの山がパンクしたりはしないが、メモリは食うので、スピードに課題がある。[br][br]そこに現れる救世主は[b]TCO（末尾呼び出し最適化）[/b]だ。[br]しかし、残念ながらPythonという言語は裏でTCOに変換してくれないので、自力でかかないといけない。[br][br]ということで、これからやるのは、[br][br][b]２．発想としての再帰[br][br]３．TCOにしてスタックがパンクしない方法[br][/b][br]このあたりを探ってみよう。[br][br]

発想としての再帰[br][br][b][size=150]＜累乗を再帰的に計算する＞[br][/size][/b][br]２の10乗の計算は2の5乗の2乗です。[br]２の5乗は２の4乗の2倍です。[br]２の4乗は２の2乗の2乗です。[br]このように、指数が偶数なら半分にして計算したものの2乗[br]指数が奇数なら1個取って偶数乗にする。[br]こうすれば、計算が速くなることは経験でわかります。[br]これをコードにしてみよう。[br][br]def fastexp(a: int, n: int) -> int:[br] if n == 0:[br] return 1[br] elif n % 2 == 1:[br] return a * fastexp(a, n - 1)[br] else:[br] t = fastexp(a, n // 2)[br] return t * t[br]fastexp(2, 100)#1267650600228229401496703205376[br]計算は一瞬で終わります。[br][br][b][size=150]＜整数の和を再帰的に定義する＞[/size][/b][br][br]自然数とはペアノの公理から１から始まる後続者です。[br]aの前の数をP(a),aの後の数をS(a)とかくと、[br]add(a,b)=｛a≠0ならadd(P(a),S(a))、a=0ならb｝と定義します。[br]たし算という計算を自然数のシフトとして再定義できるわけですね。[br][br]def add(a: int, b: int) -> int:[br] if a == 0:[br] return b[br] else:[br] return add(a - 1, b + 1)[br]add(999,1001)#2000[br]再帰的な定義の守備範囲は広そうですね。[br][br][b][size=150]＜再帰構造を再帰で調べる＞[/size][/b][br][br]たとえば、2分木TreeNodeがあり、[b]木の生成は（値value、左left、右right）[/b]とかき、[br]入れ子ができるします。これを中置方式で、左、値、右の順に[br]走査するときは、中置再帰走査[b]inorder_traversal[/b]して結果を出力します。[br]木の構造に合わせて、[b]左走査＋値＋右走査という順に分解して掘り下げ[/b]ればいいね。[br]再帰構造には再帰関数、当然といえば当然でした。[br][br]class TreeNode:[br] def __init__(self, value, left=None, right=None):[br] self.value = value[br] self.left = left[br] self.right = right[br]def inorder_traversal(node):[br] if node is None:[br] return [] # Base case: return an empty list if the node is empty.[br] return inorder_traversal(node.left) + [node.value] + inorder_traversal(node.right)[br] # 木のデータ:[br]#　　10[br]# / ＼[br]# 5　15[br]# / ＼[br]# 12 20[br]root = TreeNode(10,[br] TreeNode(5),[br] TreeNode(15, TreeNode(12), TreeNode(20)))[br]result = inorder_traversal(root)[br]print("走査結果:", result)#走査結果: [5, 10, 12, 15, 20][br][br][b][size=150]＜ソートの王者も再帰的＞[/size][/b][br]速いソートいえばクイックソートアルゴリズムですね。[br]スタートとなるリストlstがあります。[br]適当な軸となる要素pivotをリストからとります。[br]それより小さい要素群Lと残りGにわけます。[br]つまり、lst=L+pivot+Gでやめたら、第一段階の並べ替えができます。[br]しかし、こうかくと次につながりません。[br]次は、LとGそれぞれスタートとしたいので、[br]quicksort(L)+pivot+quicksort(G)を返し、ピボットだけか空のリストが終了条件です。[br]def quicksort(lst):[br] if len(lst) <= 1:[br] return lst # Base case: 対象サイズが１以下[br] else:[br] pivot = lst[0][br] less = [x for x in lst[1:] if x < pivot][br] greater_or_equal = [x for x in lst[1:] if x >= pivot][br] return quicksort(less) + [pivot] + quicksort(greater_or_equal)[br]#例[br]unsorted_list = [3, 9, 30, 90, 7, 10, 1][br]sorted_list = quicksort(unsorted_list)[br]print("クイックソート:", sorted_list)#クイックソート: [1, 3, 7, 9, 10, 30, 90][br][br][b][size=150]＜振り返り＞[br][/size][/b][br]いくらでも、例が作れそうです。このくらいにしておきます。[br][br]再帰的とは、対象を分割することで、処理対象のサイズをどんどん小さくする。[br]そして、小さくしたものをくっつけて、巨大な全体を再構成するという発想だといえるね。[br]シダの葉の細かな仕組みをつかめば、全体をつかんだことになる。[br][br]本来フラクタルでないイテラブルな対象でさえ、[br]フラクタルにとらえる発想ともいえる。[br][br][b][color=#0000ff][size=200]分割統治。[br][br]全体を部分と同型に分解することで、[br]部分を制覇すれば全体が管理できる。[br][/size][/color][/b][br][br]次はTCO、お尻から再帰する方法を見てみよう。[br][br][br]

TCO末尾呼び出し最適化TailCallOptimazition[br][br][b][size=150]＜TCOを目指す書き方＞[br][/size][/b][br]TCO末尾呼び出し最適化TailCallOptimazition[br]名前がすごいけど、普通の再帰関数の書き方だと、[br]計算式の山をつくっておき、結局終了条件という端っことから、[br]前から順に積んでいく。[br]この[b]大から小に積んで、小から大に計算する[/b]というUターンをやめて、[br][b]大から小に順に計算してしまう[/b]というストレートなやり方がTCOを「目指す」書き方だ。[br][br]階乗計算で記述の発想と表現と効率を比較してみよう。[br][br]def fact(n: int) -> int:[br] if n == 0:[br] return 1[br] else:[br] return [b]n*fact(n-1)[/b][br]これは、終了条件と再帰式がある由緒正しい再帰だ。[br][b]発想も思想もよいけど、スタックを消費する。効率が悪い[/b]。[br][br]def facti(n: int) -> int:[br] if n == 0:[br] return 1[br] f = 1[br] for i in range(2, n+1):[br] f = f * i[br] return f[br]これは、スタックがあふれないが再帰とは関係なforループ文だ。[br][b]安全をとって、思想を捨てている[/b]。[br][br]def fact_tail(n, accumulator= 1):[br] if n==0:[br] print(f"n={n},accmulator={accumulator}")[br] [b]return accumulator[/b][br] else:[br] print(f"n={n},accmulator={accumulator}")[br] [b]return fact_tail(n-1, n*accumulator)[/b][br]fact_tail(10)[br][OUT][br]n=10,accmulator=1[br]n=9,accmulator=10[br]n=8,accmulator=90[br]n=7,accmulator=720[br]n=6,accmulator=5040[br]n=5,accmulator=30240[br]n=4,accmulator=151200[br]n=3,accmulator=604800[br]n=2,accmulator=1814400[br]n=1,accmulator=3628800[br]n=0,accmulator=3628800[br]3628800[br]これは、分割統治の再帰の発想で、お尻からの累積値をつけたものだ。[br]効率もよいね。たった、11個の行で10の階乗のステップを確認できた。[br][br]しかし、「目指す」書き方ではあるけど、Pythonは内部でスタックに積んでいて[br]やはりスタックの山はパンクしてします。[br]なぜか？[br][b]import sys[br]print(sys.getrecursionlimit())[br][OUT][br]3000になっていました。[br][br][/b]だから、print文は無駄な表示をするのでコメントアウトします。[br][b]def fact_tail(n, accumulator= 1):[br] if n==0:[br] #print(f"n={n},accmulator={accumulator}")[br] return accumulator[br] else:[br] #print(f"n={n},accmulator={accumulator}")[br] return fact_tail(n-1, n*accumulator)[br]fact_tail(10)[br][/b]print(fact_tail(2000))なら、スタックの限界内だから、一瞬で表示します。[br]しかし、[br]print(fact_tail(4000))ではスタックの山が高くなりすぎて崩壊します。[br]残念すぎる。[br][br]ではどうする？[br][br][b]＜本当のTCO、呼び出し最適化とは＞[br][/b][br]先ほどのTCOを目指すやり方を本当のTCOにします。[br]つまり、スタックの山を積み上げない方法です。[br][br][color=#0000ff]import sys[br]# Pythonの文字列表現の限界を10万桁まで拡張する[br]sys.set_int_max_str_digits(100000)[br]def factorial_raw(n, accumulator=1):[br] if n == 0:[br] return accumulator[br] else:[br] return lambda: factorial_raw(n-1, n*accumulator)[br]def trampoline(thunk):[br] while callable(thunk):[br] thunk = thunk()[br] return thunk[br][br]result = trampoline(factorial_raw(4000))[br]print(f"{len(str(result))} 桁")[br]print(result)[br][OUT][br][br][/color][b]12674 桁[br]18288019515140650133147431755739190....................略[br][br][/b]スタックに積んだらすぐに下すことを繰り返すので、[br]山は積みあがりません。[br]トランポリンという関数[b]はthunkという手形関数[/b]を読み込みます。[br][b]thunkが呼び出しできる、つまり関数であるうちは読み込み実行します[/b]。[br]しかし、[br]そうでないとき、[br][b]値になったらwhileをぬけて値になったthunkを返します[/b]。[br][br]だから、トランポリンにfactorial_rawをのっけると、[br]n==0という端っこにくるまで、[b]引数のない形だけのラムダ式[/b]のため、[br][b]（n-1, n*accumulator）[/b]ここだけが進むのです。[br][br]トランポリンがカウンター代わりに値じゃないから次だ。[br]と進めてくれるからです。[br][br]トランポリンには床と空中しかありません。[br][b]０階と１階しかないスタック[/b]とも言えます。[br][br]だから、直進どころか、[b]横滑りして最後の値になってから飛ぶのやめる[/b]わけです。[br][br]いやあ、だれだか知りませんが、世の中には発想の天才がいるもんですね。

[b][size=150]＜振り返り＞[br][/size][/b][br]「再帰」愛のMITの人々やハッカーと[br]「Python」愛のグイド・ヴァンロッサムの[br][b][color=#0000ff]壮絶は戦い[/color][/b]を見ることができた感じがしました。[br][br]geogebraのiteration,iterationListはイテレタとして[br]リカージョン風のことをやるというものです。[br][br]イテレタ自体はSequence(式、変数、初期値、終わり、ステップ）という数列があります。[br][br]また、zipがあるので、Zip(式、変数１、リスト１、変数２、リスト２）でリスト演算をイテレタなしでできます。[br][br]そんな道具箱をもつgeobgebraで、再帰的な発想を深めてみよう。[br][br][b][size=150][color=#9900ff][u]課題：geogebraでfibonacci以外の再帰関数を作ってみよう。[br][/u][/color][/size][/b][br]タイトル「イテレタで再帰もどき」[br]n=slider(1,170,1)[br]「iterationListの数列[i]ak+1[/i] が[i]akとkに依存して[/i]決められる。」と[url=https://geogebra.github.io/docs/manual/ja/commands/IterationList/]マニュアル[/url]ある。[br]まずは、マニュアル例を普通の階乗計算に直す。[br][br]f(k,a)=(k+1) a[br]factorial = iterationList(f,{1,1},n)[br]とすると、1,(1+1)1, (2+1)2,......となり、自然数列の積、つまり、階乗数列が出る。[br]#{1,2,6,24,120,....}[br]a=factorial(n)[br][b]text1="factorial(" + n +")=" + a +""[br][/b][br]次に、[br]g(k)=k+1[br]h(k,a)=g(k) a [br]とすることで、kの式の自由度を上げる。[br]fact2 = iterationList(h,{1,1},n)[br]はfactorialと同じ動きをすることが確認できる。[br]i(k,a)=g(k)+ a[br]とすることで、自然数列の和、つまり、三角数列ができるかも。[br]triArray=IterationList(i,{1,1},n) #{1,3,6,10,....}[br]b=triArray(n)[br]できました。[br]text2=[b]"triangle(" + n +")=" + b +""[br][br]今度はaの式を関数にしてコラッツ数列を表示してみよう。[br][/b]start = Slider(1, 100, 1)[br]m = Slider(1, 50, 1)[br]Mod(a,2)がバグるので、floor[a/2:]== a/2として、条件分岐させます。[br]collaz(a) = If(floor[a/2:]== a/2, a / 2, 3 * a + 1)[br]collatzArray = IterationList(collaz, start, m)[br]text3="collaz(" + start+")="+collazArray[br]１に到達できるかはわかりませんが、スライダーmを動かして見よう。[br]1にとどかないときはmの設定を少し大きくするなどしてみよう。[br]お気づきの通り、終了条件のないコラッツ数列なので、１の次は４となり、[br]1,4,2,の繰り返しをしだすこともわかります。[br][br][b]終了条件がない「再帰もどきのイテレタ」の弱み[/b]ですね。[br]