SDNU1085.爬楼梯再加强版（矩阵快速幂+矩阵关系推导）

最新推荐文章于 2023-11-12 20:52:12 发布

圣帝天龙

最新推荐文章于 2023-11-12 20:52:12 发布

阅读量549

点赞数 1

CC 4.0 BY-SA版权

分类专栏：数学基础数论文章标签： ACM

本文链接：https://blog.csdn.net/qq_41658955/article/details/87289934

博客详细介绍了使用矩阵快速幂和矩阵关系进行高级爬楼梯问题的推导过程，分享了作者的心得体会。

摘要生成于 C知道，由 DeepSeek-R1 满血版支持，前往体验 >

#include<bits/stdc++.h>
using namespace std;
const int maxn=3;
const int mod=1e9+7;
struct mat
{
    long long int m[maxn][maxn];
} ans;
mat mul(mat A,mat B)
{
    mat tmp;
    for(int i=0; i<maxn; ++i)
        for(int j=0; j<maxn; ++j)
            tmp.m[i][j]=0;
    for(int i=0; i<maxn; ++i)
        for(int j=0; j<maxn; ++j)
            for(int k=0; k<maxn; ++k)
                tmp.m[i][j]=(tmp.m[i][j]+(A.m[i][k])*(B.m[k][j])%mod)%mod;
    return tmp;
}
mat quickpow(mat res,long long int N)
{
    for(int i=0; i<maxn; ++i)
        for(int j=0; j<maxn; ++j)
        {
            if(i==j)
                ans.m[i][j]=1;
            else
                ans.m[i][j]=0;
        }
    while(N)
    {
        if(N&1)
            ans=mul(ans,res);
        res=mul(res,res);
        N>>=1;
    }
    return ans;
}
int main()
{
    mat res;
    long long int n;
    while(~scanf("%lld",&n))
    {
        for(int i=0;