Haskell Advent Calendar 2023 の12日目の記事です。
『関数プログラミング入門』(オーム社)などを読んでいると foldr や foldl の融合則というのが出てきます。たとえば foldr の融合則は次の形です。
foldr の融合則
この法則を使うとたとえばリストを2回分走査する計算を1回分にまとめることができたりして便利なわけです。
が、この法則の見た目がややこしいですし、条件部分 h (f x y) = g x (h y)
も意味がよく分からないと思っていました。
これがもっと一般的な、より短い形から導けるということを知ってびっくりしたので、学んだことを復習する意味で共有したいと思います。
一般の融合則
それがこの形です
圏論の言葉にはなりましたが、式自体はすごく簡潔になりました。むりやり可換図式にするとこんな感じでしょうか(本当は三角形にしたかった)。
0id↓⏐0!A!BA↓⏐fB
この融合則で圏 C を具体的に用意してあげると foldr の融合則になるわけです。私は一般の融合則を後から知ったので、foldr の融合則がこんな当たり前っぽい式から導けるんだということにむちゃくちゃ驚きました。
F-代数
そのための準備としてまず F-代数というものを用意します。F-代数 とは、圏 C と関手 F:C→C について X←FX の形をした射のことです。
また、F-代数の間の F-準同型 というものを考えます。F-代数 X←FX から Y←FY への F-準同型 とは、射 h:X→Y であって以下を可換にするものです。
Xh↓⏐YFX↓⏐FhFY
このとき、F-代数を対象とし、F-準同型を射とする圏を考えることができます。この圏に始対象があるとき、それを F-始代数 といいます。
リスト型 [A] を F-始代数 として定義する
いま、対象 A を固定して、対象 X と射 f に対して関手 F を次のように定義します。
FX=1+A×X
Ff=id1+idA×f
ただし 1 は終対象、× と + は積と余積です。このとき、この F についての代数に始代数が存在することが示されます。都合上それを [A]←1+A×[A] と書くことにすると、始代数から他の代数への可換図式は次のように書けます。
[A](∣e,f∣)↓⏐X[nil,cons][e,f]1+A×[A]↓⏐id1+idA×(∣e,f∣)1+A×X
ここで [A] から X への射、つまり始対象からの射は [e,f] に応じて1つに決まりますので、それを ([e,f]) という記法で表しています。
ここまで始代数をまるでリスト型やその上の関数かのように [A] や nil や cons と書いていますが、そう書いてよいことは後で分かります。いま、この図の表す等式を書き下してみます。
(∣e,f∣)∘[nil,cons]=[e,f]∘(id1+idA×(∣e,f∣))
これは次の2つの式に分解されます。
- (∣e,f∣)∘nil=e∘id1
- (∣e,f∣)∘cons=f∘(idA×(∣e,f∣))
ここで、もとになっている圏 C をたとえば集合と写像の圏とします。すると nil, cons, e, f は それぞれ次の形の写像になります。ただし 1 は元を一つだけ持つ集合(一点集合)を表します。
nil:1→[A]
cons:A×[A]→[A]
e:1→X
f:A×X→X
ここで上の1番めの式より
(∣e,f∣)∘nil=e
一点集合 1 の元を () として両辺に () を適用すると、
((∣e,f∣)∘nil)()=e()
左辺を計算して nil() を nil、e() を e と書くと
(∣e,f∣)nil=e
となります。
また2番目の式の両辺に (x, xs) を適用して、
((∣e,f∣)∘cons)(x,xs)=(f∘(idA×(∣e,f∣)))(x,xs)
左辺と右辺をそれぞれ計算して
左辺
=((∣e,f∣)∘cons)(x,xs)
=(∣e,f∣)(cons(x,xs))
右辺
=f((idA×(∣e,f∣))(x,xs)
=f(x,(∣e,f∣)xs)
左辺 = 右辺の形に戻して
(∣e,f∣)(cons(x,xs))=f(x,(∣e,f∣)xs)
まとめると
となりますが、これを foldr の定義
と見比べると、 (∣e,f∣) は foldr f e
と同じ意味だと分かります。同様に、nil は []
と、cons は :
と、[A] は [a]
と同じ意味ということになります。
foldr の融合則を導く
さきほどのF-代数の圏に融合則を適用します。次の可換図式で、
[A](∣e,f∣)↓⏐Xh↓⏐Y[nil,cons][e,f][d,g]1+A×[A]↓⏐id1+idA×(∣e,f∣)1+A×X↓⏐id1+idA×h1+A×Y
一番上の [A]←1+A×[A] は始対象ですから、融合則は h∘![e,f]=![d,g] の形になります。![e,f] はつまり (∣e,f∣) のことなので、
h∘(∣e,f∣)=(∣d,g∣)
となります。また、下側の四角形の可換性を書き下すと次のようになります。
h∘[e,f]=[d,g]∘(id1+idA×h)
これは次の2つの式に分解されます。
- h∘e=d∘id1
- h∘f=g∘(idA×h)
1番の式から d=h∘e であることが分かるので、融合則は結局
h∘(∣e,f∣)=(∣h∘e,g∣)
となります。前と同様に集合と写像の圏だと思って両辺に xs を適用すると
左辺
=(h∘(∣e,f∣))xs
=h((∣e,f∣)xs)
右辺
=(∣h∘e,g∣)xs
左辺 = 右辺の形に戻して
h((∣e,f∣)xs)=(∣h∘e,g∣)xs
2番の式でも同様に両辺に (x, y) を適用すると
左辺
=(h∘f)(x,y)
=h(f(x,y))
右辺
=(g∘(idA×h))(x,y)
=g((idA×h)(x,y))
=g(x,hy)
左辺 = 右辺の形に戻して
h(f(x,y))=g(x,hy)
まとめると
となりますが、これは (∣e,f∣) が foldr f e
だったことを思い出すと、冒頭の foldr の融合則
と同じ意味であることが分かります。
まとめ
翻って考えると、foldr の融合則の条件 h (f x y) = g x (h y)
は h が準同型という意味だったんですね。実のところまだよく分かってないのですが、前よりはイメージがつかめた気がします。
参考文献
圏論勉強会:第7回 @ワークスアプリケーションズ(中村晃一)
http://nineties.github.io/category-seminar/7.html
「代数的データ型」以降のページで、上の内容のほとんどがより丁寧に書かれています。本記事でやったのはリストの catamophism の融合則が foldr の融合則と一致することを具体的に確かめたことだけです。
『関数プログラミング入門』(オーム社)
『Algorithm Design With Haskell』(Cambridge University Press)
foldr の融合則についてはこの本で知りました。
Discussion