poj 3070 矩阵快速幂简单题

最新推荐文章于 2024-10-30 23:49:49 发布

流風回雪_YZK

最新推荐文章于 2024-10-30 23:49:49 发布

阅读量1k

点赞数 1

CC 4.0 BY-SA版权

分类专栏：数学

本文链接：https://blog.youkuaiyun.com/u011498819/article/details/39322875

数学专栏收录该内容

27 篇文章

订阅专栏

本文通过一个具体的示例讲解了如何使用矩阵快速幂算法解决特定类型的问题。代码中定义了一个结构体来表示矩阵，并实现了矩阵乘法和快速幂运算。通过对输入的整数n进行处理，演示了如何高效地计算特定矩阵的幂次方。

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

基本运用，基本是模板题。

求fi【n】. (1,1) *( 1 )

( 1,0) ( 0)

#include<iostream>
#include<cstring>
using namespace std;
struct juz
{
    int bat[3][3];
    int x,y;     //行 列
};
juz mutp(juz a,juz b)
{
    juz c;
    c.x=a.x;c.y=b.y;
    memset(c.bat,0,sizeof(c.bat));
    for(int k=0;k<a.y;k++)
          for(int i=0;i<a.x;i++)
          if(a.bat[i][k])
          {
              for(int j=0;j<b.y;j++)
              {
                  c.bat[i][j]+=(a.bat[i][k]*b.bat[k][j])%10000;
              }
          }
    return c;
}
juz quickf(juz a,int k)
{
    juz c=a;
    for(int i=0;i<a.x;i++)
      for(int j=0;j<a.x;j++)
          c.bat[i][j]=(i==j);
    while(k>=1)
    {
        if(k%2)
            c=mutp(c,a);
        k=k/2; a=mutp(a,a);
    }
    return c;
}
int main()
{
    int n;
    while(cin>>n&&n!=-1)
    {
        if(n==0)
         {
            cout<<0<<endl;continue;
         }
        juz a,b,c;
        a.x=2;a.y=2; b.x=2;b.y=1;
        a.bat[0][0]=1;a.bat[0][1]=1;a.bat[1][0]=1;a.bat[1][1]=0;
        b.bat[0][0]=1;b.bat[1][0]=0;
        c=quickf(a,n-1);
        c=mutp(c,b);
        cout<<c.bat[0][0]<<endl;
    }
    return 0;
}