我有一个关于如何将'recursion'转换为'tail recursion'的问题。 这不是一个功课,只是当我试图从算法书中修改递归定理时弹出一个问题。 我熟悉使用递归(factorial和Fibonacci序列)的两个典型示例,并且还知道如何以递归方式和尾递归方式实现它们。 我的代码如下(我使用Perl只是为了简单,但可以很容易地转换为C / Java / C ++)
#this is the recursive function
sub recP {
my ($n) = @_;
if ($n==0 or $n==1 or $n==2) {
return 1;
} else {
return (recP($n-3)*recP($n-1))+1;
}
}
for (my $k=1;$k<10;$k++) {
print "*"x10,"\n";
print "recP($k)=", recP($k), "\n";
}
运行代码时,输出如下:
recP(1)=1
recP(2)=1
recP(3)=2
recP(4)=3
recP(5)=4
recP(6)=9
recP(7)=28
recP(8)=113
recP(9)=1018
递归函数在返回之前用不同的参数调用自己两次; 我尝试了几种方法将其转换为尾递归方式,但结果都是错误的。
任何人都可以查看代码并向我展示使其尾递归的正确方法吗?特别是我相信这个树递归的转换有一个例程(在返回之前多次调用递归函数),可以对此有所了解吗?所以我可以使用相同的逻辑来处理不同的问题。 提前谢谢。
答案 0 :(得分:18)
虽然您经常会看到以下内容作为将阶乘转换为尾部调用的示例:
int factorial(int n, int acc=1) {
if (n <= 1) return acc;
else return factorial(n-1, n*acc);
}
它不太正确,因为它需要乘法既是关联的又是可交换的。 (乘法是关联和可交换的,但上述内容不能作为不满足这些约束的其他操作的模型。)更好的解决方案可能是:
int factorial(int n, int k=1, int acc=1) {
if (n == 0) return acc;
else return factorial(n-1, k+1, acc*k);
}
这也可以作为斐波纳契变换的模型:
int fibonacci(int n, int a=1, int b=0) {
if (n == 0) return a;
else return fibonacci(n-1, a+b, a);
}
请注意,这些计算序列从头开始计算,而不是在调用堆栈中排队挂起的延续。因此,它们在结构上更像是迭代解决方案,而不是递归解决方案。但是,与迭代程序不同,它们从不修改任何变量;所有绑定都是不变的。这是一个有趣且有用的属性;在这些简单的情况下,它没有太大的区别,但是编写代码而不重新分配会使一些编译器优化变得更容易。
无论如何,最后一个确实提供了递归函数的模型;像斐波那契序列一样,我们需要保留相关的过去值,但我们需要其中三个而不是两个:
int mouse(int n, int a=1, int b=1, int c=1) {
if (n <=2 ) return a;
else return mouse(n-1, a*c+1, a, b);
}
<强>附加物强>
在评论中,提出了两个问题。我会尝试在这里回答它们(还有一个)。
首先,应该清楚(考虑到没有函数调用概念的底层机器架构)任何函数调用都可以被重新定义为goto(可能带有无限制的中间存储);此外,任何goto都可以表示为尾调用。因此,有可能(但不一定非常好)将任何递归重写为尾递归。
通常的机制是“延续传递风格”,这是一种奇特的说法,即每当你想调用一个函数时,你将当前函数的其余部分打包为一个新函数(“延续”),并将该延续传递给被调用的函数。由于每个函数都接收一个continuation作为参数,它必须通过调用它所接收的continuation来完成它创建的任何延续。
这可能足以使你的头旋转,所以我会用另一种方式:不是将参数和返回位置推入堆栈并调用函数(稍后将返回),而是推送参数和延续位置进入堆栈并转到一个函数,该函数稍后将转到连续位置。简而言之,您只需将堆栈作为显式参数,然后您就不需要返回。这种编程风格在事件驱动的代码中很常见(参见Python Twisted),编写(和读取)真的很痛苦。所以我强烈建议让编译器为你做这个转换,如果你能找到一个可以做到这一点。
@xxmouse 建议我从帽子中拉出递归方程式,并询问它是如何推导出来的。它只是原始的递归,但重新表述为单个元组的转换:
fn = fn-1*fn-3 + 1
=>
Fn = <Fn-11*Fn-13+1, Fn-11, Fn-12>
我不知道这是否更清楚,但这是我能做的最好的事情。请看一下斐波纳契案例,看一个稍微简单的案例。
@j_random_hacker 询问此转换的限制是什么。它适用于递归序列,其中每个元素可以由先前k
元素的某个公式表示,其中k
是常量。还有其他方法可以产生尾调用递归。例如:
// For didactic purposes only
bool is_odd(int n) { return n%2 == 1; }
int power(int x, int n, int acc=1) {
if (n == 0) return acc;
else if (is_odd(n)) return power(x, n-1, acc*x);
else return power(x*x, n/2, acc);
}
以上不与通常的非尾部调用递归相同,后者执行不同(但等效且等长)的乘法序列。
int squared(n) { return n * n; }
int power(int x, int n) {
if (n == 0) return 1;
else if (is_odd(n)) return x * power(x, n-1));
else return squared(power(x, n/2));
}
感谢Alexey Frunze进行以下测试: 输出(ideone):
mouse(0) = 1
mouse(1) = 1
mouse(2) = 1
mouse(3) = 2
mouse(4) = 3
mouse(5) = 4
mouse(6) = 9
mouse(7) = 28
mouse(8) = 113
mouse(9) = 1018
答案 1 :(得分:5)
使用谷歌,我找到了描述Tail Recursion的页面。基本上,您需要将功能分成至少两个其他功能:一个执行工作,保持&#34;积累&#34;当前值,另一个是您的workhouse函数的驱动程序。 C中的阶乘示例如下:
/* not tail recursive */
unsigned int
factorial1(unsigned int n)
{
if(n == 0)
return 1;
return n * factorial1(n-1);
}
/* tail recursive version */
unsigned int
factorial_driver(unsigned int n, unsigned int acc)
{
if(n == 0)
return acc;
/* notice that the multiplication happens in the function call */
return factorial_driver(n - 1, n * acc);
}
/* driver function for tail recursive factorial */
unsigned int
factorial2(unsigned int n)
{
return factorial_driver(n, 1);
}
答案 2 :(得分:3)
@Alexey Frunze的回答是可以的但不完全正确。确实可以将任何程序转换为一个程序,其中所有递归都是尾递归,通过将其转换为Continuation Passing Style。
我现在没有时间,但如果我有一些会议记录,我会尝试在CPS中重新实施您的计划。
答案 3 :(得分:1)
你可以这样做:
#include <stdio.h>
void fr(int n, int a[])
{
int tmp;
if (n == 0)
return;
tmp = a[0] * a[2] + 1;
a[2] = a[1];
a[1] = a[0];
a[0] = tmp;
fr(n - 1, a);
}
int f(int n)
{
int a[3] = { 1, 1, 1 };
if (n <= 2)
return 1;
fr(n - 2, a);
return a[0];
}
int main(void)
{
int k;
for (k = 0; k < 10; k++)
printf("f(%d) = %d\n", k, f(k));
return 0;
}
输出(ideone):
f(0) = 1
f(1) = 1
f(2) = 1
f(3) = 2
f(4) = 3
f(5) = 4
f(6) = 9
f(7) = 28
f(8) = 113
f(9) = 1018
编译器可能会将fr()
转换为以下内容:
void fr(int n, int a[])
{
int tmp;
label:
if (n == 0)
return;
tmp = a[0] * a[2] + 1;
a[2] = a[1];
a[1] = a[0];
a[0] = tmp;
n--;
goto label;
}
这将是尾部调用优化。
答案 4 :(得分:1)
问题是最后一个操作不是递归调用之一,而是在乘法中加1。你在C中的功能:
unsigned faa (int n) // Ordinary recursion
{
return n<3 ? 1 :
faa(n-3)*faa(n-1) + 1; // Call, call, multiply, add
}
如果更改请求值的顺序,可以将其中一个调用转换为循环:
unsigned foo (int n) // Similar to tail recursion
{ // (reverse order)
int i;
unsigned f;
for (i=3, f=1; i<=n; i++)
f = f*foo(i-3) + 1;
return f;
}
关键是考虑在原始函数中实际计算值的顺序,而不是它们被请求的顺序。
请注意,我假设您要删除一个递归调用。如果你想在函数末尾编写递归调用,期望编译器为你优化它,请参阅其他答案。
尽管如此,“The Right Thing(TM)”要做的是使用动态编程来避免多次计算相同的值:
unsigned fuu (int n) // Dynamic programming
{
int i;
unsigned A[4]={1,1,1,1};
for (i=3; i<=n; i++)
{
memmove (A+1, A, 3*sizeof(int));
A[0] = A[1]*A[3] + 1;
}
return A[0];
}
数组A包含序列的滑动窗口:A [0] == f(i),A [1] == f(i-1),A [2] == f(i-2)等等。
memmove
可能写成:
A[3] = A[2];
A[2] = A[1];
A[1] = A[0];