public static MatrixLong operator ^(MatrixLong a, long b)
{
if (a == null)
{
throw new ArgumentNullException("a");
}
if (a.Rows != a.Cols)
{
throw new ArgumentException("Can't power non-square matrix");
}
if (b < 0)
{
throw new ArgumentException("Power can't be negative");
}
MatrixLong res = new MatrixLong(a.Rows, a.Cols);
if (b == 0)
{
for (int i = 0; i < res.Rows; ++i)
{
res._Data[i, i] = 1;
}
}
else
{
res = a ^ (b >> 1);
res *= res;
if ((b & 1) != 0)
{
res *= a;
}
}
return(res);
}
public static MatrixLong operator *(MatrixLong a, MatrixLong b)
{
if (a == null)
{
throw new ArgumentNullException("a");
}
if (b == null)
{
throw new ArgumentNullException("b");
}
if (a.Cols != b.Rows)
{
throw new ArgumentException("Invalid matrix dimensions");
}
MatrixLong c = new MatrixLong(a.Rows, b.Cols);
for (int i = 0; i < a.Rows; i++)
{
for (int j = 0; j < b.Cols; ++j)
{
for (int k = 0; k < a.Cols; ++k)
{
c._Data[i, j] += NumTheoryUtils.Mul(a._Data[i, k], b._Data[k, j], MOD);
if (c._Data[i, j] >= MOD)
{
c._Data[i, j] -= MOD;
}
}
}
}
return(c);
}
