web-dev-qa-db-ja.com

高速フィボナッチ再帰

フィボナッチ再帰のアルゴリズムを思い出そうとしています。以下:

public int fibonacci(int n)  {
  if(n == 0)
    return 0;
  else if(n == 1)
    return 1;
  else
    return fibonacci(n - 1) + fibonacci(n - 2);
}

貪欲だから探しているのはnotです。これは指数関数的に大きくなります( Java再帰フィボナッチ数列 -初期引数が大きいほど、無駄な呼び出しが行われます)。

おそらく、以前のフィボナッチ値を呼び出すと、値を再度計算する代わりに取得する「循環引数シフト」のようなものがあります。

20
ducin

多分このように:

int fib(int term, int val = 1, int prev = 0)
{
 if(term == 0) return prev;
 return fib(term - 1, val+prev, val);
}

この関数は末尾再帰です。これは、非常に効率的に最適化し実行できることを意味します。実際には、単純なループに最適化されます。

40
duedl0r

この種の問題は線形回帰タイプであり、高速行列累乗法。この種のアプローチを簡潔に説明する blogpost を次に示します。

9
plesiv

memoization を使用すると、再帰フィボナッチの非常に高速なバージョンを実行できます(つまり、以前の結果を保存して再計算を回避します)。たとえば、これはPythonの概念実証です。ここでは、以前の結果を保存するために辞書が使用されます。

results = { 0:0, 1:1 }

def memofib(n):
    if n not in results:
        results[n] = memofib(n-1) + memofib(n-2)
    return results[n]

通常は「通常の」再帰バージョンをブロックする入力値に対してはすぐに戻ります。 intデータ型では大きな結果を保持するには不十分であり、任意の精度の整数を使用することをお勧めします。

まったく別のオプション-この反復バージョンを書き換えて...

def iterfib(n):
    a, b = 0, 1
    for i in xrange(n):
        a, b = b, a + b
    return a

...私のコードではloopと呼ばれる末尾再帰関数として:

def tailfib(n):
    return loop(n, 0, 1)

def loop(i, a, b):
    if i == 0:
        return a
    return loop(i-1, b, a+b)
7
Óscar López

フィボナッチ問題に関する興味深い記事を見つけました

ここにコードスニペット

# Returns F(n)
def fibonacci(n):
    if n < 0:
        raise ValueError("Negative arguments not implemented")
    return _fib(n)[0]


# Returns a Tuple (F(n), F(n+1))
def _fib(n):
    if n == 0:
        return (0, 1)
    else:
        a, b = _fib(n // 2)
        c = a * (2 * b - a)
        d = b * b + a * a
        if n % 2 == 0:
            return (c, d)
        else:
            return (d, c + d)

# added iterative version base on C# example
def iterFib(n):
    a = 0
    b = 1
    i=31
    while i>=0:
        d = a * (b * 2 - a)
        e = a * a + b * b
        a = d
        b = e
        if ((n >> i) & 1) != 0:
            c = a + b;
            a = b
            b = c
        i=i-1
    return a
3
Moch Lutfi

N番目のfib番号を取得してから、先行する番号を含む配列を作成するとします。

int a[n];
a[0] = 0;
a[1] =1;
a[i] = n[i-1]+n[n-2];
2
monkeyking

効率を高めるために再帰と遅延初期化キャッシュを使用するJavaScriptの例:

var cache = {};

function fibonacciOf (n) {
  if(n === 0) return 0;
  if(n === 1) return 1;
  var previous = cache[n-1] || fibonacciOf(n-1);
  cache[n-1] = previous;
  return previous + fibonacciOf(n-2);
};
1

duedl0rのアルゴリズムをSwiftに翻訳:

func fib(n: Int, previous: (Int, Int) = (0,1)) -> Int {
    guard n > 0 else { return 0 }
    if n == 1 { return previous.1 }
    return fib(n - 1, previous: (previous.1, previous.0 + previous.1))
}

働いた例:

fib(4)
= fib(4, (0,1) )
= fib(3, (1,1) )
= fib(2, (1,2) )
= fib(1, (2,3) )
= 3
0
Ken Ko

指数関数的な成長を止めるには、計算値を記憶する必要があります。

  1. 配列を使用して値を保存するだけです。
  2. すでに計算している場合は、配列を確認してください。
  3. 見つかった場合は、使用するか、計算して保存します。

これは、メモリを使用した再帰の高速化の実例です。

フィボナッチ数の計算

0
kta

高速なフィボナッチ計算のための良いアルゴリズムは(Pythonで)です:

def fib2(n):
    # return (fib(n), fib(n-1))
    if n ==  0: return (0,  1)
    if n == -1: return (1, -1)
    k, r = divmod(n, 2) # n=2k+r
    u_k, u_km1 = fib2(k)
    u_k_s, u_km1_s = u_k**2, u_km1**2  # Can be improved by parallel calls
    u_2kp1 = 4 * u_k_s - u_km1_s + (-2 if k%2 else 2)
    u_2km1 = u_k_s + u_km1_s
    u_2k   = u_2kp1 - u_2km1
    return (u_2kp1, u_2k) if r else (u_2k, u_2km1)

def fib(n):
    k, r = divmod(n, 2) # n=2k+r
    u_k, u_km1 = fib2(k)
    return (2*u_k+u_km1)*(2*u_k-u_km1)+(-2 if k%2 else 2) if r else u_k*(u_k+2*u_km1)

非常に高速な計算が必要な場合は、libgmpにリンクし、mpz_fib_ui()またはmpz_fib2_ui()関数を使用します。

0
HNoob