Computer Algebra in Scientific Computing

Abstract

We give a new probabilistic algorithm for interpolating a "sparse" polynomial f given by a straight-line program. Our algorithm constructs an approximation f * of f , such that f − f * probably has at most half the number of terms of f , then recurses on the difference f −f * . Our approach builds on previous work by Garg andSchost (2009), andRoche (2011), and is asymptotically more efficient in terms of the total cost of the probes required than previous methods, in many cases.