【51NOD1242】斐波那契数列的第N项(矩阵快速幂)

最新推荐文章于 2020-11-28 22:54:22 发布

oranges_c

最新推荐文章于 2020-11-28 22:54:22 发布

阅读量403

点赞数

文章标签：其他-矩阵快速幂

本文链接：https://blog.csdn.net/oranges_c/article/details/53894578

版权

其他矩阵快速幂专栏收录该内容

3 篇文章 0 订阅

订阅专栏

本文记录了一位初学者在解决51NOD1242题目的过程中，学习并应用矩阵快速幂解决斐波那契数列第N项问题的经验。在实现过程中需要注意避免因多组输入导致的时间超限错误。

摘要由CSDN通过智能技术生成

记录一个菜逼的成长。。

题目链接
ps:只是记下模板。

这题有个坑，用多组输入就TLE.

//#pragma comment(linker, "/STACK:1024000000,1024000000")
#include <cstdio>
#include <iostream>
#include <cstring>
#include <string>
#include <algorithm>
#include <cstdlib>
#include <vector>
#include <set>
#include <map>
#include <queue>
#include <stack>
#include <list>
#include <deque>
#include <cctype>
#include <bitset>
#include <cmath>
#include <cassert>
using namespace std;
#define ALL(v) (v).begin(),(v).end()
#define cl(a,b) memset(a,b,sizeof(a))
#define bp __builtin_popcount
#define pb push_back
#define mp make_pair
#define fin freopen("D://in.txt","r",stdin)
#define fout freopen("D://out.txt","w",stdout)
#define lson t<<1,l,mid
#define rson t<<1|1,mid+1,r
#define seglen (node[t].r-node[t].l+1)
#define pi 3.1415926
#define exp  2.718281828459
#define lowbit(x) (x)&(-x)
typedef long long LL;
typedef unsigned long long ULL;
typedef pair<int,int> PII;
typedef pair<LL,LL> PLL;
typedef vector<PII> VPII;
const int INF = 0x3f3f3f3f;
//const int MOD = 1e9 + 7;
template <typename T>
inline void read(T &x){
    T ans=0;
    char last=' ',ch=getchar();
    while(ch<'0' || ch>'9')last=ch,ch=getchar();
    while(ch>='0' && ch<='9')ans=ans*10+ch-'0',ch=getchar();
    if(last=='-')ans=-ans;
    x = ans;
}
const int MOD = 1000000009;
typedef vector<LL> vec;
typedef vector<vec> mat;
LL n;
mat mul(mat &A,mat &B)
{
    mat C(A.size(),vec(B[0].size()));
    for( int i = 0; i < A.size(); i++ ){
        for( int k = 0; k < B.size();k++ ){
            for( int j = 0; j < B[0].size();j++ ){
                C[i][j] = (C[i][j] + A[i][k] * B[k][j]);
                C[i][j] %= MOD;
            }
        }
    }
    return C;
}
mat pow(mat A,LL n)
{
    mat B(A.size(),vec(A.size()));
    for( int i = 0; i < A.size(); i++ ){
        B[i][i] = 1;
    }
    while(n > 0){
        if(n & 1)B = mul(B,A);
        A = mul(A,A);
        n >>= 1;
    }
    return B;
}
void solve()
{
    mat A(2,vec(2));
    A[0][0] = A[0][1] = A[1][0] = 1;
    A[1][1] = 0;
    A = pow(A,n);
    printf("%lld\n",A[1][0]);
}
int main()
{
    scanf("%lld",&n);
        solve();
    return 0;
}