Algorithms for the non-monic case of the sparse modular GCD algorithm (original) (raw)

Let G = (4y 2 + 2z)x 2 + (10y 2 + 6z) be the greatest common divisor (gcd) of two polynomials A, B ∈ Z[x, y, z]. Because G is not monic in the main variable x, the sparse modular gcd algorithm of Richard Zippel cannot be applied directly as one is unable to scale univariate images of G in x consistently. We call this the normalization problem. We present two new sparse modular gcd algorithms which solve this problem without requiring any factorizations. The first, a modification of Zippel's algorithm, treats the scaling factors as unknowns to be solved for. This leads to a structured coupled linear system for which an efficient solution is still possible. The second algorithm reconstructs the monic gcd x 2 + (5y 2 + 3z)/(2y 2 + z) from monic univariate images using a sparse, variable at a time, rational function interpolation algorithm.