Polynomials¶
Univariate polynomials over finite fields are supported with the Poly
class.
Create a polynomial¶
Create a polynomial by specifying its coefficients in degree-descending order and the finite field its over.
In [1]: GF = galois.GF(2**8)
In [2]: galois.Poly([1, 0, 0, 55, 23], field=GF)
Out[2]: Poly(x^4 + 55x + 23, GF(2^8))
Or pass a FieldArray
of coefficients without explicitly specifying the finite field.
In [3]: coeffs = GF([1, 0, 0, 55, 23]); coeffs
Out[3]: GF([ 1, 0, 0, 55, 23], order=2^8)
In [4]: galois.Poly(coeffs)
Out[4]: Poly(x^4 + 55x + 23, GF(2^8))
Use set_printoptions()
to display the polynomial coefficients in degree-ascending order.
In [5]: galois.set_printoptions(coeffs="asc")
In [6]: galois.Poly(coeffs)
Out[6]: Poly(23 + 55x + x^4, GF(2^8))
Element representation¶
As with FieldArray
instances, the finite field element representation of the polynomial coefficients may be changed
by setting the repr
keyword argument of GF()
or using the repr()
classmethod.
In [7]: GF = galois.GF(3**5)
# Display f(x) using the default integer representation
In [8]: f = galois.Poly([13, 0, 4, 2], field=GF); print(f)
13x^3 + 4x + 2
# Display f(x) using the polynomial representation
In [9]: GF.repr("poly"); print(f)
(α^2 + α + 1)x^3 + (α + 1)x + 2
# Display f(x) using the power representation
In [10]: GF.repr("power"); print(f)
(α^10)x^3 + (α^69)x + α^121
In [11]: GF.repr("int");
See Element Representation for more details.
Alternate constructors¶
There are several additional ways to create a polynomial. These alternate constructors are included as classmethods in Poly
.
By convention, alternate constructors use PascalCase
while other classmethods use snake_case
.
Create a polynomial by specifying its non-zero degrees and coefficients using Degrees()
.
In [12]: galois.Poly.Degrees([1000, 1], coeffs=[1, 179], field=GF)
Out[12]: Poly(x^1000 + 179x, GF(3^5))
Create a polynomial from its integer representation using Int()
. Additionally, one may create a polynomial from
a binary, octal, or hexadecimal string of its integer representation.
In [13]: galois.Poly.Int(268, field=GF)
Out[13]: Poly(x + 25, GF(3^5))
In [14]: galois.Poly.Int(int("0b1011", 2))
Out[14]: Poly(x^3 + x + 1, GF(2))
In [15]: galois.Poly.Int(int("0o5034", 8), field=galois.GF(2**3))
Out[15]: Poly(5x^3 + 3x + 4, GF(2^3))
In [16]: galois.Poly.Int(int("0xf700a275", 16), field=galois.GF(2**8))
Out[16]: Poly(247x^3 + 162x + 117, GF(2^8))
Create a polynomial from its string representation using Str()
.
In [17]: galois.Poly.Str("x^5 + 143", field=GF)
Out[17]: Poly(x^5 + 143, GF(3^5))
Create a polynomial from its roots using Roots()
.
In [18]: f = galois.Poly.Roots([137, 22, 51], field=GF); f
Out[18]: Poly(x^3 + 180x^2 + 19x + 58, GF(3^5))
In [19]: f.roots()
Out[19]: GF([ 22, 51, 137], order=3^5)
The Zero()
, One()
, and Identity()
classmethods create common,
simple polynomials. They are included for convenience.
In [20]: galois.Poly.Zero(GF)
Out[20]: Poly(0, GF(3^5))
In [21]: galois.Poly.One(GF)
Out[21]: Poly(1, GF(3^5))
In [22]: galois.Poly.Identity(GF)
Out[22]: Poly(x, GF(3^5))
Random polynomials of a given degree are easily created with Random()
.
In [23]: galois.Poly.Random(4, field=GF)
Out[23]: Poly(68x^4 + 37x^3 + 227x^2 + 209x + 45, GF(3^5))
Methods¶
Polynomial objects have several methods that modify or perform operations on the polynomial. Below are some examples.
Compute the derivative of a polynomial using derivative()
.
In [24]: GF = galois.GF(7)
In [25]: f = galois.Poly([1, 0, 5, 2, 3], field=GF); f
Out[25]: Poly(x^4 + 5x^2 + 2x + 3, GF(7))
In [26]: f.derivative()
Out[26]: Poly(4x^3 + 3x + 2, GF(7))
Compute the roots of a polynomial using roots()
.
In [27]: f.roots()
Out[27]: GF([5, 6], order=7)
Properties¶
Polynomial objects have several instance properties. Below are some examples.
Find the non-zero degrees and coefficients of the polynomial using nonzero_degrees
and nonzero_coeffs
.
In [28]: GF = galois.GF(7)
In [29]: f = galois.Poly([1, 0, 3], field=GF); f
Out[29]: Poly(x^2 + 3, GF(7))
In [30]: f.nonzero_degrees
Out[30]: array([2, 0])
In [31]: f.nonzero_coeffs
Out[31]: GF([1, 3], order=7)
Find the integer equivalent of the polynomial using int()
, see __int__()
. Additionally, one may
convert a polynomial into the binary, octal, or hexadecimal string of its integer representation.
In [32]: int(f)
Out[32]: 52
In [33]: g = galois.Poly([1, 0, 1, 1]); g
Out[33]: Poly(x^3 + x + 1, GF(2))
In [34]: bin(g)
Out[34]: '0b1011'
In [35]: g = galois.Poly([5, 0, 3, 4], field=galois.GF(2**3)); g
Out[35]: Poly(5x^3 + 3x + 4, GF(2^3))
In [36]: oct(g)
Out[36]: '0o5034'
In [37]: g = galois.Poly([0xf7, 0x00, 0xa2, 0x75], field=galois.GF(2**8)); g
Out[37]: Poly(247x^3 + 162x + 117, GF(2^8))
In [38]: hex(g)
Out[38]: '0xf700a275'
Get the string representation of the polynomial using str()
.
In [39]: str(f)
Out[39]: 'x^2 + 3'
Special polynomials¶
The galois
library also includes several functions to find certain special polynomials. Below are some examples.
Find one or all irreducible polynomials with irreducible_poly()
and irreducible_polys()
.
In [40]: galois.irreducible_poly(3, 3)
Out[40]: Poly(x^3 + 2x + 1, GF(3))
In [41]: list(galois.irreducible_polys(3, 3))
Out[41]:
[Poly(x^3 + 2x + 1, GF(3)),
Poly(x^3 + 2x + 2, GF(3)),
Poly(x^3 + x^2 + 2, GF(3)),
Poly(x^3 + x^2 + x + 2, GF(3)),
Poly(x^3 + x^2 + 2x + 1, GF(3)),
Poly(x^3 + 2x^2 + 1, GF(3)),
Poly(x^3 + 2x^2 + x + 1, GF(3)),
Poly(x^3 + 2x^2 + 2x + 2, GF(3))]
Find one or all primitive polynomials with primitive_poly()
and primitive_polys()
.
In [42]: galois.primitive_poly(3, 3)
Out[42]: Poly(x^3 + 2x + 1, GF(3))
In [43]: list(galois.primitive_polys(3, 3))
Out[43]:
[Poly(x^3 + 2x + 1, GF(3)),
Poly(x^3 + x^2 + 2x + 1, GF(3)),
Poly(x^3 + 2x^2 + 1, GF(3)),
Poly(x^3 + 2x^2 + x + 1, GF(3))]
Find the Conway polynomial using conway_poly()
.
In [44]: galois.conway_poly(3, 3)
Out[44]: Poly(x^3 + 2x + 1, GF(3))