Until recently, the only known method of finding the roots of polynomials over prime power rings, other than fields, was brute force. One reason for this is the lack of a division algorithm, obstructing the use of greatest common divisors. Fix a prime integer p and f in (Z/pnZ)[x] any nonzero polynomial of degree d whose coefficients are not all divisible by p. For the case n=2, we prove a new efficient algorithm to count the roots of f in Z/p2Z within time (d+size(f)+log p),2+o(1), based on a formula conjectured by Cheng, Gao, Rojas, and Wan.
more »
« less
The Weierstrass–Durand–Kerner root finder is not generally convergent
Finding roots of univariate polynomials is one of the fundamental tasks of numerics, and there is still a wide gap between root finders that are well understood in theory and those that perform well in practice. We investigate the root-finding method of Weierstrass, also known as the Durand–Kerner-method: this is a root finder that tries to approximate all roots of a given polynomial in parallel. This method has been introduced 130 years ago and has since then a good reputation for finding all roots in practice except in obvious cases of symmetry. Nonetheless, very little is known about its global dynamics and convergence properties. We show that the Weierstrass method, like the well-known Newton method, is not generally convergent: there are open sets of polynomials p p of every degree d ≥ 3 d \ge 3 such that the dynamics of the Weierstrass method applied to p p exhibits attracting periodic orbits. Specifically, all polynomials sufficiently close to Z 3 + Z + 175 Z^3 + Z + 175 have attracting cycles of period 4 4 . Here, period 4 4 is minimal: we show that for cubic polynomials, there are no periodic orbits of length 2 2 or 3 3 that attract open sets of starting points. We also establish another convergence problem for the Weierstrass method: for almost every polynomial of degree d ≥ 3 d\ge 3 there are orbits that are defined for all iterates but converge to ∞ \infty ; this is a problem that does not occur for Newton’s method. Our results are obtained by first interpreting the original problem coming from numerical mathematics in terms of higher-dimensional complex dynamics, then phrasing the question in algebraic terms in such a way that we could finally answer it by applying methods from computer algebra. The main innovation here is the translation into an algebraic question, which is amenable to (exact) computational methods close to the limits of current computer algebra systems.
more »
« less
- Award ID(s):
- 1928930
- PAR ID:
- 10420546
- Date Published:
- Journal Name:
- Mathematics of Computation
- Volume:
- 92
- Issue:
- 340
- ISSN:
- 0025-5718
- Page Range / eLocation ID:
- 839 to 866
- Format(s):
- Medium: X
- Sponsoring Org:
- National Science Foundation
More Like this
-
-
We present a computationally-efficient and numerically-robust algorithm for finding real roots of polynomials. It begins with determining the intervals where the given polynomial is monotonic. Then, it performs a robust variant of Newton iterations to find the real root within each interval, providing fast and guaranteed convergence and satisfying the given error bound, as permitted by the numerical precision used. For cubic polynomials, the algorithm is more accurate and faster than both the analytical solution and directly applying Newton iterations. It trivially extends to polynomials with arbitrary degrees, but it is limited to finding the real roots only and has quadratic worst-case complexity in terms of the polynomial's degree. We show that our method outperforms alternative polynomial solutions we tested up to degree 20. We also present an example rendering application with a known efficient numerical solution and show that our method provides faster, more accurate, and more robust solutions by solving polynomials of degree 10.more » « less
-
Abstract Fix an integer$$d\ge 2$$ . The parameters$$c_0\in \overline{\mathbb {Q}}$$ for which the unicritical polynomial$$f_{d,c}(z)=z^d+c\in \mathbb {C}[z]$$ has finite postcritical orbit, also known asMisiurewiczparameters, play a significant role in complex dynamics. Recent work of Buff, Epstein, and Koch proved the first known cases of a long-standing dynamical conjecture of Milnor using their arithmetic properties, about which relatively little is otherwise known. Continuing our work from a companion paper, we address further arithmetic properties of Misiurewicz parameters, especially the nature of the algebraic integers obtained by evaluating the polynomial defining one such parameter at a different Misiurewicz parameter. In the most challenging such combinations, we describe a connection between such algebraic integers and the multipliers of associated periodic points. As part of our considerations, we also introduce a new class of polynomials we callp-special, which may be of independent number theoretic interest.more » « less
-
Abstract The isodynamic points of a plane triangle are known to be the only pair of its centers invariant under the action of the Möbius group$${\mathcal {M}}$$ on the set of triangles, Kimberling (Encyclopedia of Triangle Centers,http://faculty.evansville.edu/ck6/encyclopedia). Generalizing this classical result, we introduce below theisodynamicmap associating to a univariate polynomial of degree$$d\ge 3$$ with at most double roots a polynomial of degree (at most)$$2d-4$$ such that this map commutes with the action of the Möbius group$${\mathcal {M}}$$ on the zero loci of the initial polynomial and its image. The roots of the image polynomial will be called theisodynamic pointsof the preimage polynomial. Our construction naturally extends from univariate polynomials to binary forms and further to their ratios.more » « less
-
We use the method of characteristic sets with respect to two term orderings to prove the existence and obtain a method of computation of a bivariate dimension polynomial associated with a non-reflexive difference-differential ideal in the algebra of difference-differential polynomials with several basic derivations and one translation. As a consequence, we obtain a new proof and a method of computation of the dimension polynomial of a non-reflexive prime difference ideal in the algebra of difference polynomials over an ordinary difference field. We also discuss applications of our results to systems of algebraic difference-differential equations.more » « less
An official website of the United States government

