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