rt

比如这样封装一个矩阵，那么我应该如何定义一个单位矩阵数组

struct Matrix
{
	ll M[N][N];
//	void clear(){memset(M,0,sizeof(M));}
//	void reset(){clear();for(int i=1;i<=L;i++)M[i][i]=1;}
	Matrix(){memset(M,0,sizeof(M));}
	Matrix(bool isT){if(isT)for(int i=1;i<=L;i++)M[i][i]=1;}
	Matrix friend operator *(const Matrix A,const Matrix B)
	{
		Matrix Ans;
		for(int i=1;i<=L;i++)
			for(int j=1;j<=L;j++)
				for(int k=1;k<=L;k++)
					Ans.M[i][j]=(Ans.M[i][j]+A.M[i][k]*B.M[k][j])%mod;
		return Ans;
	}
};