使用 Lambda 表达式编写递归三:实现 Y 组合子
2013-04-10 17:31 鹤冲天 阅读(7402) 评论(5) 编辑 收藏 举报本系列文章目录:
也许你我都难以理解,为什么有人对她痴迷疯狂,铭记在心中不说,还要刻在身上:
她让人绞尽脑汁,也琢磨不定!她让人心力憔悴,又百般回味!
她,看似平淡,却深藏玄机!她,貌不惊人,却天下无敌!
她是谁?她就是 Y 组合子:Y = λf.(λx.f (x x)) (λx.f (x x)),不动点组合子中最著名的一个。
小小开场抒情后,开始本文的内容,使用 c# 实现 Y 组合子。
本文继续使用上一篇文章中的类型假定,假定递归函数:
- 参数为 int;
- 返回值为 long。
实现 Y 组合子算法
Y 组合子 η-展开
Y 组合算子的定义:
1
|
Y = λf.(λx.f(x x)) (λx.f(x x)) |
根据前文所述,对于传值调用,应使用 η-展开后的组合子:
1
|
Y = λf.(λx.f(λv.((x x)v))) (λx.f(λv.((x x)v))) |
这个也称为 Z 组合子,是将 (x x) η-展开为 λv. ((x x) n)) 。
不过我更倾向于另外一种 η-展开形式:即将 Y 组合子中的 f (x x) η-展开为 λn. (f (x x) n)) ,最终得出:
1
|
Y = λf.(λx.λn.(f(x x)n))) (λx.λn.(f(x x)n))) |
(不知道这个组合子有没有名字,没有的话就姑且称为鹤冲天组合子吧,呵呵!)
进行一次 α-变换 变换(将参数 f 改成 g),得出:
1
|
Y = λg.(λx.λn.(g (x x) n))) (λx.λn.(g(x x)n))) |
简化 Y 组合子
仔细观察展开后的 Y,会发现:
1
|
Y = λg.(λx.λn.(g(x x)n))) (λx.λn.(g(x x)n))) |
右侧两个高亮部分是相同的,定义一个中间的变量 h:
1
|
h = λx.λn.(g(x x)n) |
则 Y 可简化为:
1
2
|
Y = λg.h h Y = λg.h(h) |
Y 可表示为 Y = g => h(h)。
h 类型确定及实现
先来确定 h 的类型。
Y = λg. h(h) 变换一步得 Y g = h(h),再变换一次 Y(g) = h(h),可以看出 h 的返回值即是 Y 的返回值。
前文中已确定 Y 返回值的类型为 Func<int, long>,所以 h 的返回值类型为 Func<int, long>。
再来确定 h 的参数类型,h 调用自身 h(h),因此 h 参数的类型就是 h 的类型,h 参数和 h 是同一类型,都是未确定的。
通过以下这个奇妙的委托,可以来表示自身调用有逻辑:
1
|
delegate TResult SelfApplicable<TResult>(SelfApplicable<TResult> self);
|
借助 SelfApplicable<TResult> ,h 类型可表示为 SelfApplicable<Func<int, long>> 。
根据 h 的定义
1
2
|
h = λx.λn.(g(x x)n) h = λx.λn.(g(x(x))(n)) |
由本系列第一篇文章中总结出的小规律,可写出 lambda表达式如下:
1
|
SelfApplicable<Func<int, long>> h = x => n => g(x(x))(n); |
说明:如果前面 Y 没有展开的话,你得出将是:
1
|
|
这样最终得出的 Y 可以编译通过,但运行时会出现死循环,最终导致 StackOverflowException。
Y 组合子实现
综合以上两小节的结论:
1
2
|
Y = λg.h(h) SelfApplicable<Func<int, long>> h = x => n => g(x(x))(n); |
可写出 Y组合子的实现代码:
1 2 3 4 |
Func<Func<Func<int, long>, Func<int, long>>, Func<int, long>> Y = g => { SelfApplicable<Func<int, long>> h = x => n => g(x(x))(n); return h(h); }; |
与之等同的方法为:
1 2 3 4 |
public static Func<int, long> Y(Func<Func<int, long>, Func<int, long>> g) { SelfApplicable<Func<int, long>> h = x => n => g(x(x))(n); return h(h); } |
Y 组合子已实现,下面我们写出几个常用音频函数。
单步函数
根据上文确定出的 g 类型,可以写出:
阶乘的单步函数:
1
|
Func<Func<int, long>, Func<int, long>> g = f => n => n == 0 ? 1 : n * f(n - 1); |
斐波那契数列求值单步函数:
1
|
Func<Func<int, long>, Func<int, long>> g = f => n => n < 2 ? n : f(n - 1) + f(n - 2); |
进行阶乘计算
综合以上代码:
1 2 3 4 5 6 7 8 9 10 |
delegate TResult SelfApplicable<TResult>(SelfApplicable<TResult> self); static void Main(string[] args) { Func<Func<Func<int, long>, Func<int, long>>, Func<int, long>> Y = g => { SelfApplicable<Func<int, long>> h = x => n => g(x(x))(n); return h(h); }; Func<int, long> fact = Y(f => n => n == 0 ? 1 : n * f(n - 1)); long result = fact(5); // 结果为 120; } |
对 Y 组合子进行封装
考虑到 Y 的复杂度及重用,可以把相关代码封装成如下:
1 2 3 4 5 6 7 |
public static class YCombitator { delegate TResult SelfApplicable<TResult>(SelfApplicable<TResult> self); public static Func<TInput, TResult> Fix<TInput, TResult>(Func<Func<TInput, TResult>, Func<TInput, TResult>> g) { SelfApplicable<Func<TInput, TResult>> h = x => n => g(x(x))(n); return h(h); } } |
使用时就方便多了:
1 2 3 4 5 |
var factorial = YCombitator.Fix<int, int>(f => n => n == 0 ? 1 : n * f(n - 1)); var result1 = factorial(5); // 120 var fibonacci = YCombitator.Fix<int, int>(f => n => n < 2 ? n : f(n - 1) + f(n - 2)); var result2 = fibonacci(5); // 5 |
总结
通过本文及前面两篇文章的努力,终于使用 c# 实现 Y 组合子,达到了使用 lambda 构造递归函数的目的。
不过,本系列还没有结束:在msdn一篇文章中,我看到这样一种写法:
1 2 3 4 5 6 7 8 9 10 |
public class Program{ delegate T SelfApplicable<T>(SelfApplicable<T> self); static void Main(string[] args) { SelfApplicable<Func<Func<Func<int,int>, Func<int,int>>, Func<int,int>>> Y = y => f => x => f(y(y)(f))(x); Func<Func<Func<int, int>, Func<int, int>>, Func<int, int>> Fix = Y(Y); Func<Func<int,int>, Func<int,int>> F = fac => x => x == 0 ? 1 : x * fac(x-1); Func<int,int> factorial = Fix(F); var result = factorial(5); // 120 } } |
看上去比本文中的实现要简单方便,怎么得出的?我会在下一篇文章中告诉大家,敬请期待!
附:相关资源:
-------------------
思想火花,照亮世界