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

#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]);
    }
}