矩阵快速幂（模板）

利用快速幂的思想 根据矩阵的结合律 可以递归二分求解 

struct Mat

{

    int mat[N][N];

};

int n;

Mat operator * (Mat a,Mat b)

{

    Mat c;

    memset(c.mat,0,sizeof(c.mat));

    int i,j,k;

    for(k =0 ; k < n ; k++)

    {

        for(i = 0 ; i < n ;i++)

        {

            if(a.mat[i][k]==0) continue;//优化

            for(j = 0 ;j < n ;j++)

            {

                if(b.mat[k][j]==0) continue;//优化

                c.mat[i][j] = (c.mat[i][j]+(a.mat[i][k]*b.mat[k][j])%mod)%mod;

            }

        }

    }

    return c;

}

Mat operator ^(Mat a,int k)

{

    Mat c;

    int i,j;

    for(i =0 ; i < n ;i++)

        for(j = 0; j < n ;j++)

        c.mat[i][j] = (i==j);

    for(; k ;k >>= 1)

    {

        if(k&1) c = c*a;

        a = a*a;

    }

    return c;

}