Fibonacci hoj 矩阵快速幂裸题

最新推荐文章于 2022-01-22 13:02:11 发布

ehi11

最新推荐文章于 2022-01-22 13:02:11 发布

阅读量737

点赞数 1

CC 4.0 BY-SA版权

分类专栏：数论 poj hoj ACM 矩阵快速幂文章标签： include

本文链接：https://blog.youkuaiyun.com/ehi11/article/details/7975444

ACM 同时被 3 个专栏收录

180 篇文章

订阅专栏

hoj

136 篇文章

订阅专栏

poj

48 篇文章

订阅专栏

本文介绍了一种高效的矩阵快速幂求解方法，并通过实例展示了其在特定场景下的应用，涉及数学运算与编程实现。

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

#include <stdio.h>
#include <cstring>
const int mod=10000;
int A[3][3],ans[3][3];
void mul(int a[][3],int b[][3])
{
    int tmp[3][3];
    for(int i=1; i<=2; i++)
        for(int j=1; j<=2; j++)
        {
            tmp[i][j]=0;
            for(int k=1; k<=2; k++)
            {
                tmp[i][j]+=a[i][k]*b[k][j];
            }
            if(tmp[i][j]>=mod) tmp[i][j]%=mod;
        }
    memcpy(a,tmp,sizeof(tmp));
}
void pow(int a[][3],int b[][3],int n)
{
    int tmp[3][3];
    memcpy(tmp,b,sizeof(int)*9);
    a[1][1]=a[2][2]=1;
    a[1][2]=a[2][1]=0;
    while(n)
    {
        if(n&1)
        mul(a,tmp);
        mul(tmp,tmp);
        n>>=1;
    }
}
int main()
{
    int n;
    while(scanf("%d",&n)==1)
    {
        if(n==-1) break;
        A[1][1]=A[1][2]=A[2][1]=1;
        A[2][2]=0;
        pow(ans,A,n);
        printf("%d\n",ans[1][2]);
    }
    return 0;
}