Starting with Numpy 1.25, only ndim-0 arrays are treated as scalars
instead of all arrays of size 1 were treated as scalars.

Python 3.11 introduced a length limitation in int -> str conversion as a
solution to CVE-2020-10735. This impacts galois to create FieldArray
with very large primes.

Example:

>>> import galois
>>> p = galois.mersenne_primes()[23]
>>> gf = galois.GF(p)

The solution does not impact unaffected versions of Python because
sys.set_int_max_str_digits() was introduced in Python 3.11.

This avoid `Literal[...]` bug, see jbms/sphinx-immaterial#275

Implements #506

Implements #491