1365 Fib(N) mod Fib(K) [斐波那契相关]

最新推荐文章于 2021-05-19 16:50:31 发布

原创最新推荐文章于 2021-05-19 16:50:31 发布 · 290 阅读

0 ·

CC 4.0 BY-SA版权

数学-斐波那契相关专栏收录该内容

0 篇文章

订阅专栏

本文深入探讨了斐波那契数列的性质，包括递推公式、负系数项以及矩阵表示，详细解释了如何求解斐波那契数列中某项对另一项取模的问题，并提供了高效的算法实现。

$Fib(N)\ mod\ Fib(K)$

Fib(N)表示斐波那契数列的第N项（ $F (0) = 0, F (1) = 1$ ），给出N和K，求 $F i b (N) m o d F i b (K)$ 。由于结果太大，输出 $Mod\ 1000000007$ 的结果。

$\color{red}{正解部分}$

首先要知道

$F_i = F_kF_{i-k+1}+F_{k-1}F_{i-k}$ ,

在 $mod\ F_k$ 意义下继续化简 $\downarrow$

$F_i\\ = F_{k-1}F_{i-k}\\ = F_{k-1}(F_kF_{i-2k+1}+F_{k-1}F_{i-k-k})\\ = F_{k-1}^2F_{i-2k}\\ \ \dotsi \\ = F_{k-1}^{\lfloor \frac{i}{k} \rfloor}F_{i\%k}$

再要知道

$F_{k-1}^2+F_{k-1}F_{k}-F_k^2 = (-1)^k$

$证明 :$

首先要知道 $二阶行列式$ 的相关内容如下.

$二阶行列式$ : $\begin{pmatrix} a\ b \\ c\ d \end{pmatrix} =ad-bc$

再看斐波那契数列的转移矩阵如下,

$\begin{bmatrix} 1\ 1 \\ 1\ 0 \end{bmatrix}^k = \begin{bmatrix} F_{k+1}\ \ \ \ F_k \\ F_k\ \ \ \ \ F_{k-1} \end{bmatrix}$

矩阵相等, 行列式也相等,

$\therefore (-1)^k \\ = F_{k+1}F_{k-1} - F_k^2\\ = F_{k-1}^2+F_{k-1}F_k-F_k^2$

$得证$ .

所以在 $\mod F_k$ 意义下, $1)^k=F_{k-1}^2$ .

再看原式: $F_i = F_{k-1}^{\lfloor \frac{i}{k} \rfloor}F_{i\%k}$

设 $\lfloor \frac{i}{k} \rfloor$ , $y=i\%k$ ,

则 $F_i = F_{k-1}^xF_{y}=\begin{cases} x\&1 == 0, \ \ \ (-1)^{\frac{x}{2}k} F_y\\ x\&1 == 1, \ \ \ (-1)^{(\frac{x}{2}+1)k} F_{y-k} = (-1)^{(\frac{x}{2}+1)k+k-y-1} F_{k-y} \end{cases}$

其中 $F_{y-k}$ 可能为负数, 这里又要引入

$斐波那契的负系数项:F_{-i} = (-1)^{i-1}F_i$ .

然后就可以 $A C$ 辣 $!$

总结一下, 该题涉及的知识点有 $↓$

$F_i = F_kF_{i-k+1}+F_{k-1}F_{i-k}$
$F_{k-1}^2+F_{k-1}F_{k}-F_k^2 = (-1)^k$
$F_{-i} = (-1)^{i-1}F_i$

$\color{red}{实现部分}$

#include<bits/stdc++.h>
#define reg register
typedef long long ll;

const int mod = 1e9 + 7;

ll N;
ll K;

struct Matrix{ int C[4][4]; Matrix(){ memset(C, 0, sizeof C); } };

Matrix modify(const Matrix &a, const Matrix &b){
        Matrix s;
        for(reg int i = 1; i <= 2; i ++)
                for(reg int j = 1; j <= 2; j ++)
                        for(reg int k = 1; k <= 2; k ++){
                                int &t = s.C[i][j];
                                t = (1ll*t + (1ll*a.C[i][k]*b.C[k][j]%mod)) % mod;
                        }
        return s;
}

Matrix ma_Ksm(Matrix a, ll b){
        Matrix s;
        for(reg int i = 1; i <= 2; i ++) s.C[i][i] = 1;
        while(b){
                if(b & 1) s = modify(s, a);
                a = modify(a, a); b >>= 1;
        }
        return s;
}

ll get(ll k){
        Matrix res, I;
        res.C[1][1] = 1;
        I.C[1][1] = I.C[1][2] = I.C[2][1] = 1;
        I = ma_Ksm(I, k);
        res = modify(res, I);
        return res.C[1][2];
}

void Work(){
        scanf("%lld%lld", &N, &K);
        ll x = N/K, y = N%K;
        if(x & 1){
                x = (x+1)/2, y = K-y-1;
                x = (K & 1) * x;
                ll Ans = ((x+y&1)?-1:1) * get(y + 1);
                if(Ans < 0) Ans = (Ans + get(K) + mod) % mod;
                printf("%lld\n", Ans);
        }else{
                x >>= 1;
                if(!(K&1)) x = 0;
                ll Ans = ((x&1)?-1:1) * get(y);
                if(Ans < 0) Ans = (Ans + get(K) + mod) % mod;
                printf("%lld\n", Ans);
        }
}

int main(){ 
        freopen("fib.in", "r", stdin);
        freopen("fib.out", "w", stdout);
        int T; scanf("%d", &T);
        while(T --) Work();
        return 0;
}