An O(M(n) log n) algorithm for the Jacobi symbol (original) (raw)

View PDF

Abstract:The best known algorithm to compute the Jacobi symbol of two n-bit integers runs in time O(M(n) log n), using Schönhage's fast continued fraction algorithm combined with an identity due to Gauss. We give a different O(M(n) log n) algorithm based on the binary recursive gcd algorithm of Stehlé and Zimmermann. Our implementation - which to our knowledge is the first to run in time O(M(n) log n) - is faster than GMP's quadratic implementation for inputs larger than about 10000 decimal digits.

Submission history

From: Richard Brent [view email]
[v1] Tue, 13 Apr 2010 02:15:09 UTC (22 KB)
[v2] Wed, 2 Jun 2010 01:16:52 UTC (20 KB)