2011-01-12 3 views
6

У меня есть рекурсия для решения.Как эффективно вычислить рекуррентные отношения в математике?

f(m,n)=Sum[f[m - 1, n - 1 - i] + f[m - 3, n - 5 - i], {i, 2, n - 2*m + 2}] + f[m - 1, n - 3] + f[m - 3, n - 7] 
f(0,n)=1, f(1,n)=n 

Однако следующий код ММА очень неэффективен

f[m_, n_] := Module[{}, 
    If[m < 0, Return[0];]; 
    If[m == 0, Return[1];]; 
    If[m == 1, Return[n];]; 
    Return[Sum[f[m - 1, n - 1 - i] + f[m - 3, n - 5 - i], {i, 2, n - 2*m + 2}] + f[m - 1, n - 3] + f[m - 3, n - 7]];] 

Он принимает нестерпимо долго, чтобы вычислить F [40,20]. Может ли кто-нибудь предложить эффективный способ сделать это? Большое спасибо!

+3

Это не «решает» рекурсию. То, о чем вы просите, это «реализация функции двух переменных, определяемых рекурсией». Решая рекурсию, потребуется найти прямую формулу для m и n, не включающую рекурсию. – ogerard

ответ

12

Стандартный трюк заключается в сохранении промежуточных значений. Следующие принимает 0,000025 секунд

f[m_, n_] := 0 /; m < 0; 
f[0, n_] := 1; 
f[1, n_] := n; 
f[m_, n_] := (f[m, n] = 
    Sum[f[m - 1, n - 1 - i] + f[m - 3, n - 5 - i], {i, 2, 
     n - 2*m + 2}] + f[m - 1, n - 3] + f[m - 3, n - 7]); 
AbsoluteTiming[f[40, 20]] 

 Смежные вопросы

  • Нет связанных вопросов^_^