矩阵乘法+快速幂求斐波那契

#include<iostream>
#include<algorithm>
#include<cmath>
using namespace std;
int n,m;
int a[2][2]={1,1,1,0};
int b[2][2]={1,0,0,1};
void jzcf(int a[2][2],int b[2][2])
{
    int i,j,k;
    int c[2][2];
    for(i=0;i<2;i++)
    {
        for(j=0;j<2;j++)
        {
            c[i][j]=0;
            for(k=0;k<2;k++)
            {
                c[i][j]=(c[i][j]+a[i][k]*b[k][j])%m;
            }
        }
    }
    for(i=0;i<2;i++)
    {
        for(j=0;j<2;j++)
        {
            b[i][j]=c[i][j];
        }
    }
}
void ksm(int n)
{
    if(n==0)
    {
        cout<<0<<endl;
        return;
    }
    else if(n==1)
    {
        cout<<1<<endl;
        return;
    }
    n--;
    while(n)
    {
        if(n%2==1)
        {
            jzcf(a,b);
        }
        jzcf(a,a);
//        cout<<a[0][0]<<endl;
        n/=2;
    }
    cout<<b[0][0]<<endl;
}
int main()
{
    int t;
    while(~scanf("%d",&t))
    {
        while(t--)
        {
            
            cin>>n>>m;
            ksm(n);
//            a[2][2]={1,1,1,0};
//            b[2][2]={1,0,0,1};
            a[0][0]=a[0][1]=a[1][0]=b[0][0]=b[1][1]=1;
            a[1][1]=b[0][1]=b[1][0]=0;
        }
    }
    return 0;
}