#include <NTL/matrix.h>
#include <NTL/vec_vec_ZZ.h>
typedef Mat<ZZ> mat_ZZ;
void add(mat_ZZ& X, const mat_ZZ& A, const mat_ZZ& B);
void sub(mat_ZZ& X, const mat_ZZ& A, const mat_ZZ& B);
void negate(mat_ZZ& X, const mat_ZZ& A);
void mul(mat_ZZ& X, const mat_ZZ& A, const mat_ZZ& B);
void mul(vec_ZZ& x, const mat_ZZ& A, const vec_ZZ& b);
void mul(vec_ZZ& x, const vec_ZZ& a, const mat_ZZ& B);
void mul(mat_ZZ& X, const mat_ZZ& A, const ZZ& b);
void mul(mat_ZZ& X, const mat_ZZ& A, long b);
void mul(mat_ZZ& X, const ZZ& a, const mat_ZZ& B);
void mul(mat_ZZ& X, long a, const mat_ZZ& B);
void determinant(ZZ& d, const mat_ZZ& A, long deterministic=0);
ZZ determinant(const mat_ZZ& a, long deterministic=0);
void solve(ZZ& d, vec_ZZ& x,
const mat_ZZ& A, const vec_ZZ& b,
long deterministic=0)
void solve1(ZZ& d, vec_ZZ& x, const mat_ZZ& A, const vec_ZZ& b);
void inv(ZZ& d, mat_ZZ& X, const mat_ZZ& A, long deterministic=0);
void sqr(mat_ZZ& X, const mat_ZZ& A);
mat_ZZ sqr(const mat_ZZ& A);
void inv(mat_ZZ& X, const mat_ZZ& A);
mat_ZZ inv(const mat_ZZ& A);
void power(mat_ZZ& X, const mat_ZZ& A, const ZZ& e);
mat_ZZ power(const mat_ZZ& A, const ZZ& e);
void power(mat_ZZ& X, const mat_ZZ& A, long e);
mat_ZZ power(const mat_ZZ& A, long e);
void ident(mat_ZZ& X, long n);
mat_ZZ ident_mat_ZZ(long n);
long IsIdent(const mat_ZZ& A, long n);
void diag(mat_ZZ& X, long n, const ZZ& d);
mat_ZZ diag(long n, const ZZ& d);
long IsDiag(const mat_ZZ& A, long n, const ZZ& d);
void transpose(mat_ZZ& X, const mat_ZZ& A);
mat_ZZ transpose(const mat_ZZ& A);
long CRT(mat_ZZ& a, ZZ& prod, const mat_zz_p& A);
void clear(mat_ZZ& a);
long IsZero(const mat_ZZ& a);
mat_ZZ operator+(const mat_ZZ& a, const mat_ZZ& b);
mat_ZZ operator-(const mat_ZZ& a, const mat_ZZ& b);
mat_ZZ operator*(const mat_ZZ& a, const mat_ZZ& b);
mat_ZZ operator-(const mat_ZZ& a);
mat_ZZ operator*(const mat_ZZ& a, const ZZ& b);
mat_ZZ operator*(const mat_ZZ& a, long b);
mat_ZZ operator*(const ZZ& a, const mat_ZZ& b);
mat_ZZ operator*(long a, const mat_ZZ& b);
vec_ZZ operator*(const mat_ZZ& a, const vec_ZZ& b);
vec_ZZ operator*(const vec_ZZ& a, const mat_ZZ& b);
mat_ZZ& operator+=(mat_ZZ& x, const mat_ZZ& a);
mat_ZZ& operator-=(mat_ZZ& x, const mat_ZZ& a);
mat_ZZ& operator*=(mat_ZZ& x, const mat_ZZ& a);
mat_ZZ& operator*=(mat_ZZ& x, const ZZ& a);
mat_ZZ& operator*=(mat_ZZ& x, long a);
vec_ZZ& operator*=(vec_ZZ& x, const mat_ZZ& a);