BZOJ_5118_Fib数列2_矩阵乘法+欧拉定理

Posted fcwww

tags:

篇首语:本文由小常识网(cha138.com)小编为大家整理,主要介绍了BZOJ_5118_Fib数列2_矩阵乘法+欧拉定理相关的知识,希望对你有一定的参考价值。

BZOJ_5118_Fib数列2_矩阵乘法+欧拉定理

Description

Fib定义为Fib(0)=0,Fib(1)=1,对于n≥2,Fib(n)=Fib(n-1)+Fib(n-2)
现给出N,求Fib(2^n).

Input

本题有多组数据。第一行一个整数T,表示数据组数。
接下来T行每行一个整数N,含义如题目所示。
n≤10^15, T≤5

Output

输出共T行,每行一个整数为所求答案。
由于答案可能过大,请将答案mod 1125899839733759后输出

Sample Input

2
2
31

Sample Output

3
343812777493853

根据欧拉定理,有$a^{n}modp=a^{nmod\varphi(p)}modp$。
然后矩阵乘法即可。需要用到快速乘
 
代码:
#include <stdio.h>
#include <string.h>
#include <algorithm>
using namespace std;
typedef long long ll;
const ll p=1125899839733759ll;
ll n;
ll qc(ll x,ll y,ll mod) {
    ll re=0;
    for(;y;y>>=1ll,x=(x+x)%mod) if(y&1ll) re=(re+x)%mod;
    return re;
}
ll qp(ll x,ll y,ll mod) {
    ll re=1;
    for(;y;y>>=1ll,x=qc(x,x,mod)) if(y&1ll) re=qc(re,x,mod);
    return re;
}
struct Mat {
    ll v[2][2];
    Mat() {memset(v,0,sizeof(v));}
    Mat operator * (const Mat &x) const {
        Mat re; int i,j,k;
        for(i=0;i<2;i++) {
            for(j=0;j<2;j++) {
                for(k=0;k<2;k++) {
                    (re.v[i][j]+=qc(v[i][k],x.v[k][j],p))%=p;
                }
            }
        }
        return re;
    }
};
Mat pow(Mat &x,ll y) {
    Mat I;
    I.v[0][0]=I.v[1][1]=1;
    while(y) {
        if(y&1ll) I=I*x;
        x=x*x;
        y>>=1ll;
    }
    return I;
}
int main() {
    int T;
    scanf("%d",&T);
    while(T--) {
        scanf("%lld",&n);
        n=qp(2,n,p-1);
        Mat x;
        x.v[0][1]=x.v[1][0]=x.v[1][1]=1;
        Mat T=pow(x,n);
        printf("%lld\n",T.v[1][0]);
    }
}

 

以上是关于BZOJ_5118_Fib数列2_矩阵乘法+欧拉定理的主要内容,如果未能解决你的问题,请参考以下文章

bzoj5118Fib数列2 费马小定理+矩阵乘法

bzoj5118: Fib数列2(费马小定理+矩阵快速幂)

[bzoj4818][Sdoi2017]序列计数_矩阵乘法_欧拉筛

矩阵乘法加速fib数列

「CodePlus 2017 12 月赛」可做题2(矩阵快速幂+exgcd)

BZOJ_5015_[Snoi2017]礼物_矩阵乘法