I found an answer for how to determine the coefficients of a polynomial $P(x)$ with all nonnegative integer coefficients, knowing its value on one or two points.
How could that method be extended to a multivariate polynomial on $n$ variables with all nonnegative integer coefficients? Or what is the most efficient method to determine the coefficients if it is possible to get the evaluation of the polynomial in any needed point?
In a special case I am currently interested in polynomials of the form:
$$P(x_1,\ldots,x_n)=\prod_{1 \le i \lt j \le n}{\left(1+x_{i}x_{j}\right)}$$