  An input- and output-sensitive GCD algorithm for multi-variate polynomials
over finite fields is proposed by combining the modular method with the
Ben-Or/Tiwari sparse interpolation. The bit complexity of the algorithm is
given and is sensitive to the sparse representation, while for previous sparse
GCD algorithms, the complexities were given only in some special cases. It is
shown that the new algorithm is superior both in theory and in practice
comparing with existing GCD algorithms: the complexity in the degree is
decreased from quadratic to linear and the running times are decreased by 1-3
orders of magnitude in various benchmarks.
