1. 程式人生 > >矩陣快速冪求斐波那契(模板)

矩陣快速冪求斐波那契(模板)

#include<bits/stdc++.h>
using namespace std;
typedef long long ll;
const int MOD=10000;
structmat
{
    lla[2][2];
};
matmat_mul(matx,maty)
{
    matres;
    memset(res.a,0,sizeof(res.a));
    for(int i=0;i<2;i++)
        for(int j=0;j<2;j++)
        for(int 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 mat_pow(int n) { matc,res; c.a[0][0]=c.a[0][1]=c.a[1][0]=1; c.a[1][1]=0; memset(res.a,0,sizeof(res.a)); for(int i=0;i<2;i++) res.a[i][i]=1; while(n) { if(n&1) res=mat_mul(res,c); c=mat_mul(c,c); n=
n>>1; } printf("%I64d\n",res.a[0][1]); } int main() { int n; while(~scanf("%d",&n)&&n!=-1) { mat_pow(n); } return 0; }