C++ acm比赛专用模板【大数矩阵快速幂斐波拉契取模模板】

最近解题遇到了 也就整理了一下 关于C++版的大数矩阵快速幂斐波拉契取模模板

#include 
#include 
#include 
using namespace std;
const int mod=10000;
typedef struct
{
    long long m[2][2];
}matrix;
matrix I={1,0,0,1};
matrix P={0,1,1,1};
matrix mul(matrix a,matrix b)
{
    int i,j,k;
    matrix c;
    for(i=0;i<2;i++)
      for(j=0;j<2;j++)
      {
          c.m[i][j]=0;
          for(k=0;k<2;k++)
            c.m[i][j]+=(a.m[i][k]*b.m[k][j])%mod;
          c.m[i][j]%=mod;
      }
    return c;
}
matrix quick_mod(int n)
{
    matrix a=P,b=I;
    while(n>0)
    {
        if(n&1)
          b=mul(b,a);
        n=n>>1;
        a=mul(a,a);
    }
    return b;
}
int main()
{
    int t;
    scanf("%d",&t);
    while(t--){
        int n;
        scanf("%d",&n);
        matrix temp=quick_mod(n-2);
        printf("%d\n",(temp.m[1][0]+temp.m[1][1])%mod);
    }
    return 0;
}

学如逆水行舟，不进则退