Maximal quotient rational reconstruction: an almost optimal algorithm for rational reconstruction

Let n/d ∈ Q, m be a positive integer and let u = n/d mod m. Thus $u$ is the image of a rational number modulo m. The rational reconstruction problem is; given u and m find n/d. A solution was first given by Wang in 1981. Wang's algorithm outputs n/d when m > 2 M2 where M = max(|n|,d). Because of the wide application of this algorithm in computer algebra, several authors have investigated its practical efficiency and asymptotic time complexity.In this paper we present a new solution which is almost optimal in the following sense; with controllable high probability, our algorithm will output n/d when m is a modest number of bits longer than 2 |n| d. This means that in a modular algorithm where m is a product of primes, the modular algorithm will need one or two primes more than the minimum necessary to reconstruct n/d; thus if |n| ⇐ d or d ⇐ |n| the new algorithm saves up to half the number of primes. Further, our algorithm will fail with high probability when m < 2 |n| d.

[1]  Joachim von zur Gathen,et al.  Modern Computer Algebra , 1998 .

[2]  Erich Kaltofen,et al.  Computing greatest common divisors and factorizations in quadratic number fields , 1989 .

[3]  Victor Y. Pan,et al.  Acceleration of Euclidean algorithm and extensions , 2002, ISSAC '02.

[4]  Paul S. Wang,et al.  Early detection of true factors in univariate polynominal factorization , 1983, EUROCAL.

[5]  Elizabeth A. Arnold,et al.  Modular algorithms for computing Gröbner bases , 2003, J. Symb. Comput..

[6]  Scott McCallum,et al.  The Computation of Polynomial Greatest Common Divisors Over an Algebraic Number Field , 1989, J. Symb. Comput..

[7]  James H. Davenport,et al.  P-adic reconstruction of rational numbers , 1982, SIGS.

[8]  P. L. Montgomery,et al.  An FFT extension of the elliptic curve method of factorization , 1992 .

[9]  Mark van Hoeij,et al.  A modular GCD algorithm over number fields presented with multiple extensions , 2002, ISSAC '02.

[10]  N. S. Barnett,et al.  Private communication , 1969 .

[11]  Arnold Schönhage,et al.  Schnelle Berechnung von Kettenbruchentwicklungen , 1971, Acta Informatica.

[12]  Tudor Jebelean Improving the Multiprecision Euclidian Algorithm , 1993, DISCO.

[13]  Donald E. Knuth,et al.  The art of computer programming. Vol.2: Seminumerical algorithms , 1981 .

[14]  Donald Ervin Knuth,et al.  The Art of Computer Programming , 1968 .

[15]  Paul S. Wang,et al.  A p-adic algorithm for univariate partial fractions , 1981, SYMSAC '81.

[16]  Lars Langemyr An Asymptotically Fast Probabilitstic Algorithm for Computing Polynominal GCD's over an Algebraic Number Field , 1990, AAECC.

[17]  George E. Collins,et al.  Efficient Rational Number Reconstruction , 1995, J. Symb. Comput..

[18]  Mark J. Encarnación Computing GCDs of Polynomials over Algebraic Number Fields , 1995, J. Symb. Comput..