tjut 2855

#include "cstdio"
#include "cstring"
#define LL long long
#define mod m
#define K 2
LL n,m;
struct Matrix
{
    LL mat[K][K];
    Matrix() {memset(mat,0,sizeof(mat));}
    Matrix(LL *val)
    {
        int idx=0;
        for(int i=0;i<K;i++)
            for(int j=0;j<K;j++)
              mat[i][j]=val[idx++];
    }
};
Matrix operator * (Matrix a,Matrix b)
{
    Matrix ret;
    for(int i=0;i<K;i++)
        for(int j=0;j<K;j++)
    {
        ret.mat[i][j]=0;
        for(int k=0;k<K;k++)
            ret.mat[i][j]+=((a.mat[i][k]*b.mat[k][j])%mod);
    }
    return ret;
}
Matrix operator ^ (Matrix a,LL n)
{
    Matrix ret,base=a;
    for(int i=0;i<K;i++) ret.mat[i][i]=1;
    while(n)
    {
        if(n&1) ret=ret*base;
        base=base*base;
        n>>=1;
    }
    return ret;
}
int main()
{
    //freopen("in.txt","r",stdin);
    int T;
    scanf("%d",&T);
    while(T--)
    {
        scanf("%I64d%I64d",&n,&m);
        if(n==0) printf("0\n");
        else
        {
            LL obj=n*2;
            LL bval[]={1,0,0,0};
            LL pval[]={1,1,1,0};
            Matrix Base(bval),Pow(pval),ans=Pow^(obj-1);
            ans=Base*ans;
            printf("%I64d\n",ans.mat[0][0]%mod);
        }
    }
}

评论
添加红包

请填写红包祝福语或标题

红包个数最小为10个

红包金额最低5元

当前余额3.43前往充值 >
需支付:10.00
成就一亿技术人!
领取后你会自动成为博主和红包主的粉丝 规则
hope_wisdom
发出的红包
实付
使用余额支付
点击重新获取
扫码支付
钱包余额 0

抵扣说明:

1.余额是钱包充值的虚拟货币,按照1:1的比例进行支付金额的抵扣。
2.余额无法直接购买下载,可以购买VIP、付费专栏及课程。

余额充值