*模板--矩阵

const int maxn=65;
struct Matrix
{
    int n,m[maxn][maxn];
    void init(int sz)
    {
        memst(m,0,sizeof(m));
        n=sz;
    }
    Matrix(int sz=0){init(sz);}
    void setE()
    {
        for(int i=0;i<n;i++) m[i][i]=1;
    }
    Matrix operator* (const Matrix &a)
    {
        Matrix ans(n);
        for(int k=0;k<n;k++)
            for(int i=0;i<n;i++)
                for(int j=0;j<n;j++)
                    ans.m[i][j]=ans.m[i][j]+m[i][k]*a.m[k][j];
        return ans;
    }
};