Sparse Polynomial Interpolation Based on Diversification

TL;DR

This paper introduces a new Monte Carlo algorithm for interpolating sparse multivariate polynomials over finite fields, improving complexity and efficiency by leveraging diversification and independent variable interpolation.

Contribution

The paper develops the first finite field interpolation algorithm with fractional power complexity in D, linear in n,T, using diversification to distinguish coefficients.

Findings

Algorithm has complexity $O^ hicksim(nT ext{log}^2q + nT ext{sqrt}D ext{log}q)$

Uses diversification to distinguish polynomial coefficients effectively

Experimental results demonstrate applicability to large degree sparse polynomials.

Abstract

We consider the problem of interpolating a sparse multivariate polynomial over a finite field, represented with a black box. Building on the algorithm of Ben-Or and Tiwari for interpolating polynomials over rings with characteristic zero, we develop a new Monte Carlo algorithm over the finite field by doing additional probes. To interpolate a polynomial $f\in F_q[x_1,\dots,x_n]$ with a partial degree bound $D$ and a term bound $T$, our new algorithm costs $O^\thicksim(nT\log ^2q+nT\sqrt{D}\log q)$ bit operations and uses $2(n+1)T$ probes to the black box. If $q\geq O(nT^2D)$, it has constant success rate to return the correct polynomial. Compared with previous algorithms over general finite field, our algorithm has better complexity in the parameters $n,T,D$ and is the first one to achieve the complexity of fractional power about $D$, while keeping linear in $n,T$. A key technique is a randomization which makes all coefficients of the unknown polynomial distinguishable, producing a diverse polynomial. This approach, called diversification, was proposed by Giesbrecht and Roche in 2011. Our algorithm interpolates each variable independently using $O(T)$ probes, and then uses the diversification to correlate terms in different images. At last, we get the exponents by solving the discrete logarithms and obtain coefficients by solving a linear system. We have implemented our algorithm in Maple. Experimental results shows that our algorithm can applied to sparse polynomials with large degree. We also analyze the success rate of the algorithm.