Fibonacci 前N项平方和矩阵快速幂

最新推荐文章于 2024-12-06 00:48:37 发布

原创最新推荐文章于 2024-12-06 00:48:37 发布 · 1.8k 阅读

3 ·

CC 4.0 BY-SA版权

本文介绍了一种使用矩阵快速幂的方法来高效计算斐波那契数列的前n项平方和。通过将斐波那契数列的递推公式转化为矩阵形式，实现了对大数的有效处理，并提供了完整的AC代码实现。

由普通快速幂同理需要一个单位变量S=1

但是矩阵快速幂是矩阵相乘所以需要构造一个单位矩阵

由于本人没有学习过线性代数对矩阵方面不是很敏感所以摸索Fn项是矩阵的哪一个位置很久

这题其他知识相信大家在网上都找的到

但这题真正考得是如何把 Fn的前n项平方和转化

由图可知前N项平方和可以转化成 Fn*(Fn+Fn-1)

AC代码如下

#include<iostream>
#include<algorithm>
#include<string.h>
#include<stdio.h>
#include <queue>
#include<cstdlib>
#define ll long long
#define MOD 1000000007
using namespace std;
struct nobe
{
    ll a[2][2];
};
ll n;
ll sum;
nobe mut(nobe x,nobe y)
{
    nobe res;
    memset(res.a,0,sizeof(res.a));
    for(ll i=0;i<2;i++)
        for(ll j=0;j<2;j++)
        for(ll k=0;k<2;k++)
            res.a[i][j]=(res.a[i][j]+x.a[i][k]*y.a[k][j])%MOD;
        return res;
}
void quick(ll n)
{
    nobe c,res;
    c.a[0][0]=1;
    c.a[0][1]=1;
    c.a[1][0]=1;
    c.a[1][1]=0;
    memset(res.a,0,sizeof(res.a));
    for(ll i=0;i<2;i++)
        res.a[i][i]=1;
    while(n)
    {
        if(n&1)
        res=mut(res,c);
        c=mut(c,c);
        n=n>>1;
    }
       printf("%lld\n",((res.a[0][1]%MOD)*(res.a[0][1]%MOD+res.a[1][1]%MOD))%MOD);

}
int main()
{
    ll i,t;
    while(scanf("%lld",&n)!=EOF)
    {
           quick(n);
    }


    return 0;
}