51nod 1358 浮波那契

最新推荐文章于 2022-02-11 13:58:53 发布

YYyyCCCcccBb

最新推荐文章于 2022-02-11 13:58:53 发布

阅读量206

点赞数

分类专栏： 51nod基础题矩阵

本文链接：https://blog.csdn.net/qq_33951440/article/details/78436558

版权

51nod基础题同时被 2 个专栏收录

100 篇文章 0 订阅

订阅专栏

矩阵

16 篇文章 0 订阅

订阅专栏

1358 浮波那契

基准时间限制：1 秒空间限制：131072 KB 分值: 40 难度：4级算法题

关注

TengBieBie已经学习了很多关于斐波那切数列的性质，所以他感到一些些厌烦。现在他遇到了一个新的数列，这个数列叫做Float-Bonacci。这里有一个关于Float-Bonacci的定义。

对于一个具体的n，TengBieBie想要快速计算FB(n).

但是TengBieBie对FB的了解非常少，所以他向你求助。

你的任务是计算FB(n).FB(n)可能非常大，请输出FB(n)%1,000,000,007 (1e9+7)即可。

Input

输入共一行，在一行中给出一个整数n (1<=n<=1,000,000,000)。

Output

对于每一个n，在一行中输出FB(n)%1,000,000,007 (1e9+7)。

Input示例

Output示例

思路：

和别的题特殊的位置就在于这个3.4。只要处理好3.4 就可以解题了。

可以将G[n] = FB[n/10]

这样就可以展开来算。FB[X]= G[X*10]

这样G[n][1] = G[n-1][10]+G[n-1][34]

#include<bits/stdc++.h>
using namespace std;
typedef long long ll;
const int N=36;
ll mod=1e9+7;
struct node{
    ll a[N][N];
    void init()
    {
        for(int i=1;i<N;i++)
        {
            for(int j=1;j<N;j++)
                a[i][j]=0;
        }
    }
}  ;
int n,m;

node mul(node x,node y)
{
    node res;
    res.init();
    for(int i=1;i<N;i++)
    {
        for(int j=1;j<N;j++)
        {
            for(int k=1;k<N;k++)
            {
                res.a[i][j]+=x.a[i][k] * y.a[k][j];
                res.a[i][j]%=mod;
            }
        }
    }
    return res;
}
node _quick( node base, ll m)
{
    node res;
    res.init();
    for(int i=1;i<N;i++)
    {
        res.a[i][i]=1;
    }
    while(m){
        if(m&1)
            res=mul(base,res);
        base=mul(base,base);
        m/=2;
    }
    return res;
}
int main()
{
    int n;
    scanf("%d",&n);
    if(n<=4)
    {
        printf("1\n");
        return 0;
    }
    node base;
    base.init();
    for(int i=2;i<=34;i++)
        base.a[i][i-1]=1;
    base.a[1][10]=base.a[1][34]=1;
    base=_quick(base,(ll)(n-4)*10);
    ll res=0;
    for(int i=1;i<35;i++)
        res=(res+base.a[1][i])%mod;
    printf("%lld\n",res);
}