#pragma once //table-driven galois field modulo 2 //do not use with GF(2^17) or larger namespace nall { template struct GaloisField { using type = GaloisField; GaloisField(u32 x = 0) : x(x) {} operator field() const { return x; } auto operator^(field y) const -> type { return x ^ y; } auto operator+(field y) const -> type { return x ^ y; } auto operator-(field y) const -> type { return x ^ y; } auto operator*(field y) const -> type { return x && y ? exp(log(x) + log(y)) : 0; } auto operator/(field y) const -> type { return x && y ? exp(log(x) + Elements - log(y)) : 0; } auto& operator =(field y) { return x = y, *this; } auto& operator^=(field y) { return x = operator^(y), *this; } auto& operator+=(field y) { return x = operator^(y), *this; } auto& operator-=(field y) { return x = operator^(y), *this; } auto& operator*=(field y) { return x = operator*(y), *this; } auto& operator/=(field y) { return x = operator/(y), *this; } auto pow(field y) const -> type { return exp(log(x) * y); } auto inv() const -> type { return exp(Elements - log(x)); } // 1/x static auto log(u32 x) -> u32 { enum : u32 { Size = bit::round(Elements), Mask = Size - 1 }; static array log = [] { u32 shift = 0, polynomial = Polynomial; while(polynomial >>= 1) shift++; shift--; array log; field x = 1; for(u32 n : range(Elements)) { log[x] = n; x = x << 1 ^ (x >> shift ? Polynomial : 0); } log[0] = 0; //-inf (undefined) return log; }(); return log[x & Mask]; } static auto exp(u32 x) -> u32 { static array exp = [] { u32 shift = 0, polynomial = Polynomial; while(polynomial >>= 1) shift++; shift--; array exp; field x = 1; for(u32 n : range(Elements)) { exp[n] = x; x = x << 1 ^ (x >> shift ? Polynomial : 0); } return exp; }(); return exp[x % Elements]; } field x; }; }