マトリックス高速べき乗テンプレート

typedef vector vec;
typedef vector mat;

mat mul(mat &A,mat& B)
{
    //    A m*n B n*h   C=A*B  m*h
    mat C(A.size(),vec(B[0].size()));
    for(int i=0;i0){
        if(n&1) B=mul(B,A);
        A=mul(A,A);
        n>>=1;
    }
    return B;
}