#include <NTL/matrix.h>
#include <NTL/vec_vec_ZZ_p.h>
typedef Mat<ZZ_p> mat_ZZ_p;
void add(mat_ZZ_p& X, const mat_ZZ_p& A, const mat_ZZ_p& B);
void sub(mat_ZZ_p& X, const mat_ZZ_p& A, const mat_ZZ_p& B);
void negate(mat_ZZ_p& X, const mat_ZZ_p& A);
void mul(mat_ZZ_p& X, const mat_ZZ_p& A, const mat_ZZ_p& B);
void mul(vec_ZZ_p& x, const mat_ZZ_p& A, const vec_ZZ_p& b);
void mul(vec_ZZ_p& x, const vec_ZZ_p& a, const mat_ZZ_p& B);
void mul(mat_ZZ_p& X, const mat_ZZ_p& A, const ZZ_p& b);
void mul(mat_ZZ_p& X, const mat_ZZ_p& A, long b);
void mul(mat_ZZ_p& X, const ZZ_p& a, const mat_ZZ_p& B);
void mul(mat_ZZ_p& X, long a, const mat_ZZ_p& B);
void determinant(ZZ_p& d, const mat_ZZ_p& A);
ZZ_p determinant(const mat_ZZ_p& a);
void transpose(mat_ZZ_p& X, const mat_ZZ_p& A);
mat_ZZ_p transpose(const mat_ZZ_p& A);
void solve(ZZ_p& d, vec_ZZ_p& x, const mat_ZZ_p& A, const vec_ZZ_p& b);
void solve(zz_p& d, const mat_zz_p& A, vec_zz_p& x, const vec_zz_p& b);
void inv(ZZ_p& d, mat_ZZ_p& X, const mat_ZZ_p& A);
void sqr(mat_ZZ_p& X, const mat_ZZ_p& A);
mat_ZZ_p sqr(const mat_ZZ_p& A);
void inv(mat_ZZ_p& X, const mat_ZZ_p& A);
mat_ZZ_p inv(const mat_ZZ_p& A);
void power(mat_ZZ_p& X, const mat_ZZ_p& A, const ZZ& e);
mat_ZZ_p power(const mat_ZZ_p& A, const ZZ& e);
void power(mat_ZZ_p& X, const mat_ZZ_p& A, long e);
mat_ZZ_p power(const mat_ZZ_p& A, long e);
void ident(mat_ZZ_p& X, long n);
mat_ZZ_p ident_mat_ZZ_p(long n);
long IsIdent(const mat_ZZ_p& A, long n);
void diag(mat_ZZ_p& X, long n, const ZZ_p& d);
mat_ZZ_p diag(long n, const ZZ_p& d);
long IsDiag(const mat_ZZ_p& A, long n, const ZZ_p& d);
void random(mat_ZZ_p& x, long n, long m);
mat_ZZ_p random_mat_ZZ_p(long n, long m);
long gauss(mat_ZZ_p& M);
long gauss(mat_ZZ_p& M, long w);
void image(mat_ZZ_p& X, const mat_ZZ_p& A);
void kernel(mat_ZZ_p& X, const mat_ZZ_p& A);
void clear(mat_ZZ_p& a);
long IsZero(const mat_ZZ_p& a);
mat_ZZ_p operator+(const mat_ZZ_p& a, const mat_ZZ_p& b);
mat_ZZ_p operator-(const mat_ZZ_p& a, const mat_ZZ_p& b);
mat_ZZ_p operator*(const mat_ZZ_p& a, const mat_ZZ_p& b);
mat_ZZ_p operator-(const mat_ZZ_p& a);
mat_ZZ_p operator*(const mat_ZZ_p& a, const ZZ_p& b);
mat_ZZ_p operator*(const mat_ZZ_p& a, long b);
mat_ZZ_p operator*(const ZZ_p& a, const mat_ZZ_p& b);
mat_ZZ_p operator*(long a, const mat_ZZ_p& b);
vec_ZZ_p operator*(const mat_ZZ_p& a, const vec_ZZ_p& b);
vec_ZZ_p operator*(const vec_ZZ_p& a, const mat_ZZ_p& b);
mat_ZZ_p& operator+=(mat_ZZ_p& x, const mat_ZZ_p& a);
mat_ZZ_p& operator-=(mat_ZZ_p& x, const mat_ZZ_p& a);
mat_ZZ_p& operator*=(mat_ZZ_p& x, const mat_ZZ_p& a);
mat_ZZ_p& operator*=(mat_ZZ_p& x, const ZZ_p& a);
mat_ZZ_p& operator*=(mat_ZZ_p& x, long a);
vec_ZZ_p& operator*=(vec_ZZ_p& x, const mat_ZZ_p& a);