#include <NTL/matrix.h>
#include <NTL/vec_vec_GF2.h>
typedef Mat<GF2> mat_GF2;
void conv(mat_GF2& X, const vec_vec_GF2& A);
mat_GF2 to_mat_GF2(const vec_vec_GF2& A);
void add(mat_GF2& X, const mat_GF2& A, const mat_GF2& B);
void sub(mat_GF2& X, const mat_GF2& A, const mat_GF2& B);
void negate(mat_GF2& X, const mat_GF2& A);
void mul(mat_GF2& X, const mat_GF2& A, const mat_GF2& B);
void mul(vec_GF2& x, const mat_GF2& A, const vec_GF2& b);
void mul(vec_GF2& x, const vec_GF2& a, const mat_GF2& B);
void mul(mat_GF2& X, const mat_GF2& A, GF2 b);
void mul(mat_GF2& X, const mat_GF2& A, long b);
void mul(mat_GF2& X, GF2 a, const mat_GF2& B);
void mul(mat_GF2& X, long a, const mat_GF2& B);
void determinant(GF2& d, const mat_GF2& A);
GF2 determinant(const mat_GF2& A);
void transpose(mat_GF2& X, const mat_GF2& A);
mat_GF2 transpose(const mat_GF2& A);
void solve(GF2& d, vec_GF2& x, const mat_GF2& A, const vec_GF2& b);
void solve(GF2& d, const mat_GF2& A, vec_GF2& x, const vec_GF2& b);
void inv(GF2& d, mat_GF2& X, const mat_GF2& A);
void sqr(mat_GF2& X, const mat_GF2& A);
mat_GF2 sqr(const mat_GF2& A);
void inv(mat_GF2& X, const mat_GF2& A);
mat_GF2 inv(const mat_GF2& A);
void power(mat_GF2& X, const mat_GF2& A, const ZZ& e);
mat_GF2 power(const mat_GF2& A, const ZZ& e);
void power(mat_GF2& X, const mat_GF2& A, long e);
mat_GF2 power(const mat_GF2& A, long e);
void ident(mat_GF2& X, long n);
mat_GF2 ident_mat_GF2(long n);
long IsIdent(const mat_GF2& A, long n);
void diag(mat_GF2& X, long n, GF2 d);
mat_GF2 diag(long n, GF2 d);
long IsDiag(const mat_GF2& A, long n, long d);
void random(mat_GF2& x, long n, long m);
mat_GF2 random_mat_GF2(long n, long m);
long gauss(mat_GF2& M);
long gauss(mat_GF2& M, long w);
void image(mat_GF2& X, const mat_GF2& A);
void kernel(mat_GF2& X, const mat_GF2& A);
void clear(mat_GF2& X);
long IsZero(const mat_GF2& A);
mat_GF2 operator+(const mat_GF2& a, const mat_GF2& b);
mat_GF2 operator-(const mat_GF2& a, const mat_GF2& b);
mat_GF2 operator*(const mat_GF2& a, const mat_GF2& b);
mat_GF2 operator-(const mat_GF2& a);
mat_GF2 operator*(const mat_GF2& a, GF2 b);
mat_GF2 operator*(const mat_GF2& a, long b);
mat_GF2 operator*(GF2 a, const mat_GF2& b);
mat_GF2 operator*(long a, const mat_GF2& b);
vec_GF2 operator*(const mat_GF2& a, const vec_GF2& b);
vec_GF2 operator*(const vec_GF2& a, const mat_GF2& b);
mat_GF2& operator+=(mat_GF2& x, const mat_GF2& a);
mat_GF2& operator-=(mat_GF2& x, const mat_GF2& a);
mat_GF2& operator*=(mat_GF2& x, const mat_GF2& a);
mat_GF2& operator*=(mat_GF2& x, GF2 a);
mat_GF2& operator*=(mat_GF2& x, long a);
vec_GF2& operator*=(vec_GF2& x, const mat_GF2& a);