Teresa Krick

We study the problem of representing multivariate polynomials with rational coefficients, which are nonnegative and strictly positive on finite semialgebraic sets, using rational sums of squares. We focus on the case of finite semialgebraic sets S S defined by equality constraints, generating a zero-dimensional ideal I I , and by nonnegative sign constraints. First, we obtain existential results. We prove that a strictly positive polynomial f f with coefficients in a subfield K K of R \mathbb {R} has a representation in terms of weighted Sums-of-Squares with coefficients in this field, even if the ideal I I is not radical. We generalize this result to the case where f f is nonnegative on S S and ( f ) + ( I : f ) = 1 (f) + (I:f)=1 . We deduce that nonnegative polynomials with coefficients in K K can be represented in terms of Sum-of-Squares of polynomials with coefficients in K K , when the ideal is radical. Second, we obtain degree bounds for such Sums-of-Squares representations, which depend linearly on the regularity of the ideal and the degree of the defining equations, when they form a graded basis. Finally, we analyze the bit complexity of the Sums-of-Squares representations for polynomials with coefficients in Q \mathbb {Q} , in the case the ideal is radical. The bitsize bounds are linear in the bitsize of the input polynomials, quadratic or cubic in the Bezout bound, and linear in the regularity, generalizing and improving previous results obtained for special zero dimensional ideals. As an application in the context of polynomial optimization, we retrieve and improve results on the finite convergence and exactness of the moment/Sums-of-Squares hierarchy.

Given rational univariate polynomials f and g such that gcd(f, g) and f / gcd(f, g) are relatively prime, we show that g is non-negative on all the real roots of f if and only if g is a sum of squares of rational polynomials modulo f. We complete our study by exhibiting an algorithm that produces a certificate that a polynomial g is non-negative on the real roots of a non-zero polynomial f , when the above assumption is satisfied.

We generalize Sylvester single sums to multisets (sets with repeated elements), and show that these sums compute subresultants of two univariate polyomials as a function of their roots independently of their multiplicity structure. This is the first closed formula for subresultants in terms of roots that works for arbitrary polynomials, previous efforts only handled special cases. Our extension involves in some cases confluent Schur polynomials, and is obtained by using multivariate symmetric interpolation via an Exchange Lemma.

We describe and analyze a numerical algorithm for computing the homology (Betti numbers and torsion coefficients) of real projective varieties. Here numerical means that the algorithm is numerically stable (in a sense to be made precise). Its cost depends on the condition of the input as well as on its size and is singly exponential in the number of variables (the dimension of the ambient space) and polynomial in the condition and the degrees of the defining polynomials. In addition, we show that outside of an exceptional set of measure exponentially small in the size of the data, the algorithm takes exponential time.

We provide explicit formulae for the coefficients of the order-d polynomial subresultant of (x-\alpha)^m and (x-\beta)^n with respect to the set of Bernstein polynomials \{(x-\alpha)^j(x-\beta)^{d-j}, \, 0\le j\le d\}. They are given by hypergeometric expressions arising from determinants of binomial Hankel matrices.

The theory of symmetric multivariate Lagrange interpolation is a beautiful but rather unknown tool that has many applications. Here we derive from it an Exchange Lemma that allows to explain in a simple and natural way the full description of the double sum expressions introduced by Sylvester in 1853 in terms of subresultants and their Bézout coefficients.

The theory of symmetric multivariate Lagrange interpolation is a beautiful but rather unknown tool that has many applications. Here we derive from it an Exchange Lemma that allows to explain in a simple and natural way the full description of the double sum expressions introduced by Sylvester in 1853 in terms of subresultants and their Bezout coefficients.