Towards Verified Polynomial Factorisation

James H. Davenport

2024-09-15

Abstract:Computer algebra systems are really good at factoring polynomials, i.e. writing f as a product of irreducible factors. It is relatively easy to verify that we have a factorisation, but verifying that these factors are irreducible is a much harder problem. This paper reports work-in-progress to do such verification in Lean.

Symbolic Computation

What problem does this paper attempt to address?

Integer Polynomial Factorization by Recombination of Real Factors: Re-evaluating an Old Technique in Modern Era

Shahriar Iravanian

2024-10-21

Abstract:Polynomial factorization over $ZZ$ is of great historical and practical importance. Currently, the standard technique is to factor the polynomial over finite fields first and then to lift to integers. Factorization over finite fields can be done in polynomial time using Berlekamp or Cantor-Zassenhaus algorithms. Lifting from the finite field to $ZZ$ requires a combinatorial algorithm. The van Hoeij algorithm casts the combinatorial problem as a knapsack-equivalent problem, which is then solved using lattice-reduction (the LLL algorithm) in polynomial time, which is implemented in many computer algebra systems (CAS). In this paper, we revisit the old idea of starting with factorization over $RR$ instead of a finite field, followed by recombination of the resulting linear and quadratic factors. We transform the problem into an integer subset sum problem, which is then solved using the Horowizt-Sinha algorithm. This algorithm can factor a random integer polynomial of degree $d$ in a time complexity of $O(2^(d slash 4))$. While the resulting algorithm is exponential, consistent with the integer subset sum problem being in NP, it has a few advantages. First, it is simple and easy to implement. Second, it is almost embarrassingly parallelizable. We demonstrate this by implementing the algorithm in a Graphic Processing Unit (GPU). The resulting code can factor a degree 100 polynomial is a few tenths of a second, comparable to some standard CAS. This shows that it is possible to use current hardware, especially massively parallel systems like GPU, to the benefit of symbolic algebra.

Symbolic Computation,Computational Complexity
Derandomizing Multivariate Polynomial Factoring for Low Degree Factors

Pranjal Dutta,Amit Sinhababu,Thomas Thierauf

2024-11-26

Abstract:For a polynomial $f$ from a class $\mathcal{C}$ of polynomials, we show that the problem to compute all the constant degree irreducible factors of $f$ reduces in polynomial time to polynomial identity tests (PIT) for class $\mathcal{C}$ and divisibility tests of $f$ by constant degree polynomials. We apply the result to several classes $\mathcal{C}$ and obtain the constant degree factors in \begin{enumerate} \item polynomial time, for $\mathcal{C}$ being polynomials that have only constant degree factors, \item quasipolynomial time, for $\mathcal{C}$ being sparse polynomials, \item subexponential time, for $\mathcal{C}$ being polynomials that have constant-depth circuits. \end{enumerate} Result 2 and 3 were already shown by Kumar, Ramanathan, and Saptharishi with a different proof and their time complexities necessarily depend on black-box PITs for a related bigger class $\mathcal{C}'$. Our complexities vary on whether the input is given as a blackbox or whitebox. We also show that the problem to compute the sparse factors of polynomial from a class $\mathcal{C}$ reduces in polynomial time to PIT for class $\mathcal{C}$, divisibility tests of $f$ by sparse polynomials, and irreducibility preserving bivariate projections for sparse polynomials. For $\mathcal{C}$ being sparse polynomials, it follows that it suffices to derandomize irreducibility preserving bivariate projections for sparse polynomials in order to compute all the sparse irreducible factors efficiently. When we consider factors of sparse polynomials that are sums of univariate polynomials, a subclass of sparse polynomials, we obtain a polynomial time algorithm. This was already shown by Volkovich with a different proof.

Computational Complexity
Polynomial Factorization Over Henselian Fields

Maria Alberich-Carramiñana,Jordi Guàrdia,Enric Nart,Adrien Poteaux,Joaquim Roé,Martin Weimann

DOI: https://doi.org/10.1007/s10208-024-09646-x

IF: 3.4388

2024-02-22

Foundations of Computational Mathematics

Abstract:We present an algorithm that, given an irreducible polynomial g over a general valued field ( K , v ), finds the factorization of g over the Henselianization of K under certain conditions. The analysis leading to the algorithm follows the footsteps of Ore, Mac Lane, Okutsu, Montes, Vaquié and Herrera–Olalla–Mahboub–Spivakovsky, whose work we review in our context. The correctness is based on a key new result (Theorem 4.10), exhibiting relations between generalized Newton polygons and factorization in the context of an arbitrary valuation. This allows us to develop a polynomial factorization algorithm and an irreducibility test that go beyond the classical discrete, rank-one case. These foundational results may find applications for various computational tasks involved in arithmetic of function fields, desingularization of hypersurfaces, multivariate Puiseux series or valuation theory.

mathematics, applied,computer science, theory & methods
First steps towards Computational Polynomials in Lean

James Harold Davenport

2024-09-15

Abstract:The proof assistant Lean has support for abstract polynomials, but this is not necessarily the same as support for computations with polynomials. Lean is also a functional programming language, so it should be possible to implement computational polynomials in Lean. It turns out not to be as easy as the naive author thought.

Symbolic Computation
Exact Bivariate Polynomial Factorization over ℚ by Approximation of Roots

Yong Feng,Wenyuan Wu,Jingzhong Zhang,Jingwei Chen

DOI: https://doi.org/10.1007/s11424-014-2170-5

2015-01-01

Journal of Systems Science and Complexity

Abstract:Factorization of polynomials is one of the foundations of symbolic computation. Its applications arise in numerous branches of mathematics and other sciences. However, the present advanced programming languages such as C++ and J++, do not support symbolic computation directly. Hence, it leads to difficulties in applying factorization in engineering fields. In this paper, the authors present an algorithm which use numerical method to obtain exact factors of a bivariate polynomial with rational coefficients. The proposed method can be directly implemented in efficient programming language such C++ together with the GNU Multiple-Precision Library. In addition, the numerical computation part often only requires double precision and is easily parallelizable.
Factoring sparse polynomials fast

Alexander Demin,Joris van der Hoeven

2023-12-29

Abstract:Consider a sparse polynomial in several variables given explicitly as a sum of non-zero terms with coefficients in an effective field. In this paper, we present several algorithms for factoring such polynomials and related tasks (such as gcd computation, square-free factorization, content-free factorization, and root extraction). Our methods are all based on sparse interpolation, but follow two main lines of attack: iteration on the number of variables and more direct reductions to the univariate or bivariate case. We present detailed probabilistic complexity bounds in terms of the complexity of sparse interpolation and evaluation.

Symbolic Computation,Commutative Algebra
Exact Bivariate Polynomial Factorization Over A"E By Approximation Of Roots

Feng Yong,Wu Wenyuan,Zhang Jingzhong,Chen Jingwei

DOI: https://doi.org/10.1007/s11424-014-2170-5

2015-01-01

Abstract:Factorization of polynomials is one of the foundations of symbolic computation. Its applications arise in numerous branches of mathematics and other sciences. However, the present advanced programming languages such as C++ and J++, do not support symbolic computation directly. Hence, it leads to difficulties in applying factorization in engineering fields. In this paper, the authors present an algorithm which use numerical method to obtain exact factors of a bivariate polynomial with rational coefficients. The proposed method can be directly implemented in efficient programming language such C++ together with the GNU Multiple-Precision Library. In addition, the numerical computation part often only requires double precision and is easily parallelizable.
Exact polynomial factorization by approximate high degree algebraic numbers

Jing-wei Chen,Yong Feng,Xiao-lin Qin,Jing-zhong Zhang

DOI: https://doi.org/10.1145/1577190.1577199

2009-01-01

Abstract:For factoring polynomials in two variables with rational coefficients, an algorithm using transcendental evaluation was presented by Hulst and Lenstra. In their algorithm, transcendence measure was computed. However, a constant c is necessary to compute the transcendence measure. The size of c involved the transcendence measure can influence the efficiency of the algorithm greatly. In this paper, we overcome the problem arising in Hulst and Lenstra's algorithm and propose a new polynomial time algorithm for factoring bivariate polynomials with rational coefficients. Using an approximate algebraic number of high degree instead of a variable of a bivariate polynomial, we can get a univariate one. A factor of the resulting univariate polynomial can then be obtained by a numerical root finder and the purely numerical LLL algorithm. The high degree of the algebraic number guarantees that this factor corresponds to a factor of the original bivariate polynomial. We prove that our algorithm saves a (log2(mn))2+ε factor in bit-complexity comparing with the algorithm presented by Hulst and Lenstra, where (n, m) represents the bi-degree of the polynomial to be factored. We also demonstrate on many significant experiments that our algorithm is practical. Moreover our algorithm can be generalized to polynomials with variables more than two.
Some factorization results on polynomials having integer coefficients

Jitender Singh,Rishu Garg

DOI: https://doi.org/10.1080/00927872.2023.2301530

2024-01-15

Communications in Algebra

Abstract:In this article, we prove some factorization results for several classes of polynomials having integer coefficients, which in particular yield several classes of irreducible polynomials. Such classes of polynomials are devised by imposing some sufficiency conditions on their coefficients along with some conditions on the prime factorization of the constant term or the leading coefficient of the underlying polynomials in conjunction with the information about their root location.

mathematics
Multivariate Polynomial Factorization by Interpolation Method

Jingzhong Zhang,Yong Feng

DOI: https://doi.org/10.48550/arxiv.math/0701670

2007-01-01

Abstract:Factorization of polynomials arises in numerous areas in symbolic computation. It is an important capability in many symbolic and algebraic computation. There are two type of factorization of polynomials. One is convention polynomial factorization, and the other approximate polynomial factorization. Conventional factorization algorithms use symbolic methods to get exact factors of a polynomial while approximate factorization algorithms use numerical methods to get approximate factors of a polynomial. Symbolic computation often confront intermediate expression swell problem, which lower the efficiency of factorization. The numerical computation is famous for its high efficiency, but it only gives approximate results. In this paper, we present an algorithm which use approximate method to get exact factors of a multivariate polynomial. Compared with other methods, this method has the numerical computation advantage of high efficiency for some class of polynomials with factors of lower degree. The experimental results show that the method is more efficient than {\it factor} in Maple 9.5 for polynomials with more variables and higher degree.
Formalizing Factorization on Euclidean Domains and Abstract Euclidean Algorithms

Thaynara Arielly de Lima,Andréia Borges Avelar,André Luiz Galdino,Mauricio Ayala-Rincón

DOI: https://doi.org/10.4204/EPTCS.402.5

2024-04-23

Abstract:This paper discusses the extension of the Prototype Verification System (PVS) sub-theory for rings, part of the PVS algebra theory, with theorems related to the division algorithm for Euclidean rings and Unique Factorization Domains that are general structures where an analog of the Fundamental Theorem of Arithmetic holds. First, we formalize the general abstract notions of divisibility, prime, and irreducible elements in commutative rings, essential to deal with unique factorization domains. Then, we formalize the landmark theorem, establishing that every principal ideal domain is a unique factorization domain. Finally, we specify the theory of Euclidean domains and formally verify that the rings of integers, the Gaussian integers, and arbitrary fields are Euclidean domains. To highlight the benefits of such a general abstract discipline of formalization, we specify a Euclidean gcd algorithm for Euclidean domains and formalize its correctness. Also, we show how this correctness is inherited under adequate parameterizations for the structures of integers and Gaussian integers.

Logic in Computer Science
Exact Bivariate Polynomial Factorization in Q by Approximation of Roots

Yong Feng,Wenyuan Wu,Jingzhong Zhang

DOI: https://doi.org/10.48550/arXiv.1008.3596

2010-08-21

Abstract:Factorization of polynomials is one of the foundations of symbolic computation. Its applications arise in numerous branches of mathematics and other sciences. However, the present advanced programming languages such as C++ and J++, do not support symbolic computation directly. Hence, it leads to difficulties in applying factorization in engineering fields. In this paper, we present an algorithm which use numerical method to obtain exact factors of a bivariate polynomial with rational coefficients. Our method can be directly implemented in efficient programming language such C++ together with the GNU Multiple-Precision Library. In addition, the numerical computation part often only requires double precision and is easily parallelizable.

Algebraic Geometry,Symbolic Computation,Numerical Analysis
The numerical factorization of polynomials

Wenyuan Wu,Zhonggang Zeng

DOI: https://doi.org/10.1007/s10208-015-9289-1

2021-03-09

Abstract:Polynomial factorization in conventional sense is an ill-posed problem due to its discontinuity with respect to coefficient perturbations, making it a challenge for numerical computation using empirical data. As a regularization, this paper formulates the notion of numerical factorization based on the geometry of polynomial spaces and the stratification of factorization manifolds. Furthermore, this paper establishes the existence, uniqueness, Lipschitz continuity, condition number, and convergence of the numerical factorization to the underlying exact factorization, leading to a robust and efficient algorithm with a MATLAB implementation capable of accurate polynomial factorizations using floating point arithmetic even if the coefficients are perturbed.

Numerical Analysis
Efficiently factoring polynomials modulo $p^4$

Ashish Dwivedi,Rajat Mittal,Nitin Saxena

DOI: https://doi.org/10.48550/arXiv.1901.06628

2019-01-20

Abstract:Polynomial factoring has famous practical algorithms over fields-- finite, rational \& $p$-adic. However, modulo prime powers it gets hard as there is non-unique factorization and a combinatorial blowup ensues. For example, $x^2+p \bmod p^2$ is irreducible, but $x^2+px \bmod p^2$ has exponentially many factors! We present the first randomized poly(deg $f, \log p$) time algorithm to factor a given univariate integral $f(x)$ modulo $p^k$, for a prime $p$ and $k \leq 4$. Thus, we solve the open question of factoring modulo $p^3$ posed in (Sircana, ISSAC'17). Our method reduces the general problem of factoring $f(x) \bmod p^k$ to that of {\em root finding} in a related polynomial $E(y) \bmod\langle p^k, \varphi(x)^\ell \rangle$ for some irreducible $\varphi \bmod p$. We could efficiently solve the latter for $k\le4$, by incrementally transforming $E(y)$. Moreover, we discover an efficient and strong generalization of Hensel lifting to lift factors of $f(x) \bmod p$ to those $\bmod\ p^4$ (if possible). This was previously unknown, as the case of repeated factors of $f(x) \bmod p$ forbids classical Hensel lifting.

Computational Complexity,Data Structures and Algorithms,Symbolic Computation,Number Theory
Cyclotomic Factors and LRS-Degeneracy

John Abbott,Nico Mexis

2024-03-14

Abstract:We present three new, practical algorithms for polynomials in $\mathbb{Z}[x]$: one to test if a polynomial is cyclotomic, one to determine which cyclotomic polynomials are factors, and one to determine whether the given polynomial is LRS-degenerate. A polynomial is ``LRS-degenerate'' iff it has two distinct roots $\alpha, \beta$ such that $\beta = \zeta \alpha$ for some root of unity $\zeta$. All three algorithms are based on ``intelligent brute force''. The first two produce the indexes of the cyclotomic polynomials; the third produces a list of degeneracy orders. The algorithms are implemented in CoCoALib.

Commutative Algebra,Number Theory
On Factor Left Prime Factorization Problems for Multivariate Polynomial Matrices

Dong Lu,Dingkang Wang,Fanghui Xiao

DOI: https://doi.org/10.48550/arXiv.2010.07007

2020-10-14

Abstract:This paper is concerned with factor left prime factorization problems for multivariate polynomial matrices without full row rank. We propose a necessary and sufficient condition for the existence of factor left prime factorizations of a class of multivariate polynomial matrices, and then design an algorithm to compute all factor left prime factorizations if they exist. We implement the algorithm on the computer algebra system Maple, and two examples are given to illustrate the effectiveness of the algorithm. The results presented in this paper are also true for the existence of factor right prime factorizations of multivariate polynomial matrices without full column rank.

Symbolic Computation
An algorithm for factoring non-monic quadratic polynomials Or: How I learned to stop using the quadratic formula and love undoing FOIL

Corey Thomas Bruns

DOI: https://doi.org/10.48550/arXiv.0912.1296

2009-12-08

Abstract:We give an algorithm for factoring quadratic polynomials over any UFD, Z in particular. We prove the correctness of this algorithm and give examples over Z and Z[i].

Rings and Algebras,History and Overview
Towards Optimal Bounds for Implicit Factorization Problem

Yao Lu,Liqiang Peng,Rui Zhang,Lei Hu,Dongdai Lin

DOI: https://doi.org/10.1007/978-3-319-31301-6_26

2016-01-01

Abstract:We propose a new algorithm to solve the Implicit Factorization Problem, which was introduced by May and Ritzenhofen at PKC'09. In 2011, Sarkar and Maitra IEEE TIT 576: 4002---4013, 2011 improved May and Ritzenhofen's results by making use of the technique for solving multivariate approximate common divisors problem. In this paper, based on the observation that the desired root of the equations that derived by Sarkar and Maitra contains large prime factors, which are already determined by some known integers, we develop new techniques to acquire better bounds. We show that our attack is the best among all known attacks, and give experimental results to verify the correctness. Additionally, for the first time, we can experimentally handle the implicit factorization for the case of balanced RSA moduli.
Topology and Factorization of Polynomials

Hani Shaker

DOI: https://doi.org/10.48550/arXiv.0704.3363

2008-04-02

Abstract:For any polynomial $P \in \mathbb{C}[X_1,X_2,...,X_n]$, we describe a $\mathbb{C}$-vector space $F(P)$ of solutions of a linear system of equations coming from some algebraic partial differential equations such that the dimension of $F(P)$ is the number of irreducible factors of $P$. Moreover, the knowledge of $F(P)$ gives a complete factorization of the polynomial $P$ by taking gcd's. This generalizes previous results by Ruppert and Gao in the case $n=2$.

Algebraic Geometry,Algebraic Topology
SAT and Lattice Reduction for Integer Factorization

Yameen Ajani,Curtis Bright

2024-06-29

Abstract:The difficulty of factoring large integers into primes is the basis for cryptosystems such as RSA. Due to the widespread popularity of RSA, there have been many proposed attacks on the factorization problem such as side-channel attacks where some bits of the prime factors are available. When enough bits of the prime factors are known, two methods that are effective at solving the factorization problem are satisfiability (SAT) solvers and Coppersmith's method. The SAT approach reduces the factorization problem to a Boolean satisfiability problem, while Coppersmith's approach uses lattice basis reduction. Both methods have their advantages, but they also have their limitations: Coppersmith's method does not apply when the known bit positions are randomized, while SAT-based methods can take advantage of known bits in arbitrary locations, but have no knowledge of the algebraic structure exploited by Coppersmith's method. In this paper we describe a new hybrid SAT and computer algebra approach to efficiently solve random leaked-bit factorization problems. Specifically, Coppersmith's method is invoked by a SAT solver to determine whether a partial bit assignment can be extended to a complete assignment. Our hybrid implementation solves random leaked-bit factorization problems significantly faster than either a pure SAT or pure computer algebra approach.

Cryptography and Security,Discrete Mathematics,Logic in Computer Science,Symbolic Computation

Towards Verified Polynomial Factorisation

Integer Polynomial Factorization by Recombination of Real Factors: Re-evaluating an Old Technique in Modern Era

Derandomizing Multivariate Polynomial Factoring for Low Degree Factors

Polynomial Factorization Over Henselian Fields

First steps towards Computational Polynomials in Lean

Exact Bivariate Polynomial Factorization over ℚ by Approximation of Roots

Factoring sparse polynomials fast

Exact Bivariate Polynomial Factorization Over A"E By Approximation Of Roots

Exact polynomial factorization by approximate high degree algebraic numbers

Some factorization results on polynomials having integer coefficients

Multivariate Polynomial Factorization by Interpolation Method

Formalizing Factorization on Euclidean Domains and Abstract Euclidean Algorithms

Exact Bivariate Polynomial Factorization in Q by Approximation of Roots

The numerical factorization of polynomials

Efficiently factoring polynomials modulo $p^4$

Cyclotomic Factors and LRS-Degeneracy

On Factor Left Prime Factorization Problems for Multivariate Polynomial Matrices

An algorithm for factoring non-monic quadratic polynomials Or: How I learned to stop using the quadratic formula and love undoing FOIL

Towards Optimal Bounds for Implicit Factorization Problem

Topology and Factorization of Polynomials

SAT and Lattice Reduction for Integer Factorization