使用Lambda表达式编写递归函数
2009-08-31 01:40 by 老赵, 12633 visits其实这从来不是一个很简单的事情，虽然有些朋友认为这很简单。
伪递归
我的意思是，某些看上去像递归的做法，事实上并非是递归，所以我把它叫做是“伪”递归。
例如，我们想要使用Lambda表达式编写一个计算递归的fac函数，一开始我们总会设法这样做：
Func<int, int> fac = x => x <= 1 ? 1 : x * fac(x - 1);
不过此时编译器会无情地告诉我们，fac还没有定义。于是您可能会想，这个简单，分两行写咯。于是有朋友就会给出这样的代码：
Func<int, int> fac = null; fac = x => x <= 1 ? 1 : x * fac(x - 1);
这样看起来也很“递归”，执行起来似乎也没有问题。但是，其实这并没有使用Lambda表达式构造一个递归函数，为什么呢？因为我们使用Lambda表达式构造的其实只是一个普通的匿名方法，它是这样的：
x => x <= 1 ? 1 : x * fac(x - 1);
既然是“匿名方法”，这个构造的东西是没有名字的——因此用Lambda表达式写递归“从来不是一个很简单的事情”。那么这个Lambda表达式里的fac是什么呢？是一个“委托”。因此，这只是个“调用了一个委托”的Lambda表达式。“委托对象”和“匿名方法”是有区别的，前者是一个实际的对象，而后者只是个“定义方式”，只是“委托对象”可以成为“匿名方法”的载体而已。这个Lambda表达式构造的“委托对象”在调用时，它会去寻找fac这个引用所指向的委托对象。请注意，这里是根据“引用”去找“对象”，这意味着Lambda表达式构造的委托对象在调用时，fac可能已经不再指向当初的委托对象了。例如：
Func<int, int> fac = null; fac = x => x <= 1 ? 1 : x * fac(x - 1); Console.WriteLine(fac(5)); // 120; Func<int, int> facAlias = fac; fac = x => x; Console.WriteLine(facAlias(5)); // 20
第一次打印出的120是正确的结果。不过facAlias从fac那里“接过”了使用Lambda表达式构造的委托对象之后，我们让fac引用指向了新的匿名方法x => x。于是facAlias在调用时：
facAlias(5) <— facAlias是x => x <= 1 ? 1 : x * fac(x – 1) = 5 <= 1 ? 1 : 5 * fac(5 - 1) = 5 * fac(4) <— 注意此时fac是x => x = 5 * 4 = 20
自然就不对了。
因此，使用Lambda表达式构造一个递归函数不是一件容易的事情。
把自己传给自己吧
可能已经有朋友知道“标准”的做法是什么样的，不过我这里还想谈一下我当时遇到这个问题时想到的一个做法。比较笨（非常符合我的特点），但是可以解决问题。
我的想法是，既然使用“Lambda表达式来构造一个递归函数”的难点是因为“我们正在构造的东西是没有名字的”，因此“我们无法调用自身”。那么，如果我们换种写法，把我们正在调用的匿名函数作为参数传给自己，那么不就可以在匿名函数的方法体中，通过调用参数来调用自身了吗？于是，原本我们构造fac方法的Lambda表达式：
x => x <= 1 ? 1 : x * fac(x - 1);
就需要变成：
(f, x) => x <= 1 ? 1 : x * f(f, x - 1);
请注意，这里的f参数是一个函数，它也是我们正在使用Lambda表达式定义的匿名委托（就是“(f, x) => ...”这一长串）。为了递归调用，它还必须把自身作为第一个参数传入下一层的调用中去，所以最后不是fac(x - 1)而是f(f, x - 1)。我们可以把这个匿名函数放到一个叫做selfFac的变量中去：
var selfFac = (f, x) => x <= 1 ? 1 : x * f(f, x - 1);
在第一次调用selfFac时，我们必须把它自身传递进去。于是我们可以这样来获得阶乘的结果：
Console.WriteLine(selfFac(selfFac, 5)); // 120;
但是这段代码没法编译通过，因为编译器不知道selfFac应该是什么类型的委托对象。不过根据selfFac(selfFac, 5)的调用方式，我们可以推断出，这个委托类型会接受两个参数，第一个是它自身的类型，第二个是个整型，而返回的也是个整型。于是，我们可以得出委托的签名了：
delegate int SelfFactorial(SelfFactorial selfFac, int x);
哎，但是这一点都不通用啊。没关系，我们可以写的通用一些：
delegate TResult SelfApplicable<T, TResult>(SelfApplicable<T, TResult> self, T arg);
这样，我们便可以定义selfFac，甚至于selfFib（菲波纳契数列）：
SelfApplicable<int, int> selfFib = (f, x) => x <= 1 ? 1 : f(f, x - 1) + f(f, x - 2); Console.WriteLine(selfFib(selfFib, 5)); // 8
但是，这还不是我们所需要的递归函数啊。没错，我们需要的是传入一个x就可以得到结果的函数，这种每次还需要把自己传进去的东西算什么？不过这个倒也容易，在selfXxx的基础上再前进一步就可以了：
SelfApplicable<int, int> selfFac = (f, x) => x <= 1 ? 1 : x * f(f, x - 1); Func<int, int> fac = x => selfFac(selfFac, x); SelfApplicable<int, int> selfFib = (f, x) => x <= 1 ? 1 : f(f, x - 1) + f(f, x - 2); Func<int, int> fib = x => selfFib(selfFib, x);
为此，我们甚至可以总结出一个辅助方法：
static Func<T, TResult> Make<T, TResult>(SelfApplicable<T, TResult> self) { return x => self(self, x); }
于是乎：
var fac = Make<int, int>((f, x) => x <= 1 ? 1 : x * f(f, x - 1)); var fib = Make<int, int>((f, x) => x <= 1 ? 1 : f(f, x - 1) + f(f, x - 2));
这样我们便使用Lambda表达式定义了递归函数。当然，需要两个参数的递归函数定义方式也比较类似。首先是SelfApplicable委托类型和对应的辅助方法：
// 委托类型 delegate TResult SelfApplicable<T1, T2, TResult>(SelfApplicable<T1, T2, TResult> self, T1 arg1, T2 arg2); // 辅助方法 static Func<T1, T2, TResult> Make<T1, T2, TResult>(SelfApplicable<T1, T2, TResult> self) { return (x, y) => self(self, x, y); }
于是使用“辗转相除法”计算最大公约数的gcd函数便是：
var gcd = Make<int, int, int>((f, x, y) => y == 0 ? x : f(f, y, x % y)); Console.WriteLine(gcd(20, 36)); // 4
这也是我目前凭“个人能力”能够走出的最远距离了。
不动点组合子
但是装配脑袋很早给了我们更好的解决方法：
static Func<T, TResult> Fix<T, TResult>(Func<Func<T, TResult>, Func<T, TResult>> f) { return x => f(Fix(f))(x); } static Func<T1, T2, TResult> Fix<T1, T2, TResult>(Func<Func<T1, T2, TResult>, Func<T1, T2, TResult>> f) { return (x, y) => f(Fix(f))(x, y); }
Fix求出的是函数f的不动点，它就是我们所需要的递归函数：
var fac = Fix<int, int>(f => x => x <= 1 ? 1 : x * f(x - 1)); var fib = Fix<int, int>(f => x => x <= 1 ? 1 : f(x - 1) + f(x - 2)); var gcd = Fix<int, int, int>(f => (x, y) => y == 0 ? x : f(y, x % y));
用脑袋的话来说，Fix方法应该被视为是内置方法。您比较Fix方法内部和之前的Make方法内部的写法，就能够意识到两种做法之间的差距了。
由于我的脑袋不如装配脑袋的脑袋装配的那么好，即使看来一些推导过程之后还是无法做到100%的理解，我还需要阅读更多的内容。希望在以后的某一天，我可以把这部分内容融会贯通地理解下来，并且可以详细地解释给大家听。在这之前，我还是听脑袋的话，把Fix强行记在脑袋里吧。
最后，希望大家多多参与一些如“函数式链表快速排序”这样的趣味编程——不过，千万不要学脑袋这样做：
var qsort = Fix<IEnumerable<int>, IEnumerable<int>>(f => l => l.Any() ? f(l.Skip(1).Where(e => e < l.First())).Concat(Enumerable.Repeat(l.First(), 1)).Concat(f(l.Skip(1).Where(e => e >= l.First()))) : Enumerable.Empty<int>());
当然，偶尔玩玩是有益无害的。
占个沙发