矩阵加法，减法，乘法c++

代码如下

struct matrix{
    int a[maxn][maxn];
    matrix() { memset(a,0,sizeof a); }//init
    matrix operator-(const matrix &T) const
    {
        matrix res;
        for(int i=1;i<=n;i++)
            for(int j=1;j<=n;j++)
                res.a[i][j]=(a[i][j]-T.a[i][j]+mod)%mod;
        return res;
    }
    matrix operator+(const matrix &T) const
    {
        matrix res;
        for(int i=1;i<=n;i++)
            for(int j=1;j<=n;j++)
                res.a[i][j]=(a[i][j]+T.a[i][j])%mod;
        return res;
    }
    matrix operator*(const matrix &T) const//n*p X p*m
    {
        matrix res;
        for(int i=1;i<=n;i++)
            for(int j=1;j<=m;j++)
                for(int k=1;k<=p;k++)
                    res.a[i][j]=(res.a[i][j]+a[i][k]*T.a[k][j]+mod*mod)%mod;
        return res;
    }
}