cs.SCAug 28, 2026

Generic Characteristic-Zero Equivalence Between Derivative Bézout Inversion and Multipoint Evaluation

Authors: Zijian Zeng

Organizations: Institute of Computer Science and Digital Innovation, UCSI University, Kuala Lumpur, 56000, MALAYSIA

Abstract

Let a1,…,ama_1,\ldots,a_m be distinct elements of a field KK, and let Z(X)=∏i=1m(X−ai)Z(X)=\prod_{i=1}^m (X-a_i). We study the arithmetic complexity of computing the unique normalized Bezout pair s,ts,t satisfying sZ+tZ′=1sZ+tZ'=1, with °t<m°t<m and °s<m−1°s<m-1. The classical product-tree approach requires O(MK(m)log⁡m)O(M_K(m)\log m) field operations, where MK(m)M_K(m) denotes the cost of multiplying degree-<m<m polynomials over KK. Thus, even when MK(m)=O(mlog⁡m)M_K(m)=O(m\log m), the resulting bound is O(mlog⁡2m)O(m\log^2 m) rather than O(mlog⁡m)O(m\log m). Over an infinite field of characteristic zero, we prove that, in the generic rational straight-line-program model, computing all coefficients of the canonical Bezout pair is equivalent, up to an additive O(MK(m))O(M_K(m)) cost, to arbitrary-node multipoint polynomial evaluation and to interpolation. The main ingredient is an explicit differential reconstruction that recovers ZZ from (s,t)(s,t) in O(MK(m))O(M_K(m)) arithmetic operations on a nonempty Zariski-open subset. Combining this reconstruction with automatic differentiation and transposition yields the complexity equivalence. We further transfer Strassen's lower bound for the elementary symmetric functions to the Bezout problem, obtaining an Ω(mlog⁡m)Ω(m\log m) nonscalar lower bound. Hence an O(mlog⁡m)O(m\log m) algorithm, if it exists in this model, would be asymptotically optimal. The reconstruction is genuinely characteristic-dependent: in characteristic pp, all squarefree polynomials Xp+cX+dX^p+cX+d, with fixed c≠0c\ne 0, have the same normalized Bezout pair (0,c−1)(0,c^{-1}). Therefore the all-field, all-input O(mlog⁡m)O(m\log m) problem remains open; generically in characteristic zero, however, it is reduced to the corresponding arbitrary-node multipoint evaluation and interpolation problem.

Figures & tables

Explore similar work

Sep 4, 2026

Representation Redundancy and Structural Complexity in Finite-Field Inversion

The representation chosen for a mathematical operation can affect both its algebraic form and its empirical learning difficulty. We study this phenomenon for inversion over F2n\mathbb F_{2^n}, with field elements expressed in varying ordered F2\mathbb F_2-bases. We prove that two ordered bases induce the same coordinate inversion map if and only if they belong to the same Galois orbit. Since every orbit has size nn, the correspondence between ordered bases and distinct inversion maps is exactly nn-to-one. We then analyze three Boolean formulations of inversion. The reference formulation has algebraic degree n−1n-1 and joint ANF leap 11, the mixed representation formulation has degree 2(n−1)2(n-1) and joint ANF leap 22, and the complete raw formulation has degree at most 3(n−1)3(n-1) and joint ANF leap at least nn. Exhaustive computations agree with the theoretical results and bounds in the cases considered. Controlled experiments with multilayer perceptrons show the same ordering in learning difficulty, while Galois orbit redundancy provides only a limited generalization benefit under the tested conditions. These results show that exact redundancy among representations can coexist with changes in Boolean structure and learning behavior when the representation is exposed as part of the input.
Sep 11, 2026

Algebraic Expressivity Certificates for Shallow Polynomial Neural Networks

We study exact representability by bias-free shallow polynomial neural networks using algebraic geometry. Over C\mathbb{C}, a width-rr network with activation z↦zdz\mapsto z^d computes a sum of rr dd-th powers of linear forms, whose Zariski closure is a Veronese secant variety. Ideal elimination therefore yields polynomial certificates of nonrepresentability. We implement this construction as a generic architecture-to-certificate pipeline. For quadratics, we recover the exact symmetric determinantal description and explain its dimension through orthogonal symmetry. In higher degree, the implementation recovers classical catalecticant and secant equations and maps the practical reach of direct elimination across a finite architecture sweep. We also derive the exact population loss floor for a rank-two quadratic network on the sphere, illustrating how an algebraic obstruction induces irreducible approximation error.
Aug 23, 2026

From Symmetry to Invariance: Learning Galois Equivalent Representations in Finite Fields

Neural networks can learn algebraic operations from finite examples, but it remains unclear whether this ability transfers across mathematically equivalent representations of the same operation. We study this question through multiplication in finite fields under changes of basis. The Galois action organizes basis representations into orbits, and bases in the same orbit induce the same coordinate multiplication map. This structure allows us to separate learning multiplication from transferring it to basis representations that are not used for training. We examine several ways of providing or recovering the relevant orbit structure, including invariant labels, basis matrices, orbit recognition, and algebraic decomposition. Our main approach trains a model to predict the Galois action between basis representations. Repeated applications of the learned transformation are then used to construct a canonical representative for each orbit, which supports multiplication on held-out bases through exact canonical matching. This provides a concrete mechanism for converting a learned algebraic symmetry into an invariant representation that can be used for transfer.