n階常系數(shù)線性齊次遞推關(guān)系的解法有很多，特征方程法、生成函數(shù)法，但是對(duì)于編程最實(shí)用的是矩陣解法。
我們定義所要求的f(n) = Ak-1 * f(n - 1) + Ak-2 * f(n - 2) + ... + A0 * f(n-k)，其中f(0)...f(k-1)的初值已經(jīng)給好。
構(gòu)造k * k的矩陣M：
其中A =(Ak-1 Ak-2 ... A1)，I是單位矩陣。
然后構(gòu)造一個(gè)k * 1列向量b：
這樣，M * b之后b0的值就是f(k)，以此類推，M ^ n * b之后b0的值就是f(k-1+n)，算法復(fù)雜度O(k ^ 3 * logn)。