Search

In this paper, we compute an algebraic decomposition of black-box rings in the generic ring model. More precisely, we explicitly decompose a black-box ring as a direct product of a nilpotent black-box ring and unital local black-box rings, by computing all its primitive idempotents. The algorithm presented in this paper uses quantum subroutines for the computation of the p -power parts of a black-box ring and then classical algorithms for the computation of the corresponding primitive idempotents. As a by-product, we get that the reduction of a black-box ring is also a black-box ring. The first application of this decomposition is an extension of the work of Maurer and Raub (2007) on representation problem in black-box finite fields to the case of reduced p -power black-box rings. Another important application is an IND-CCA 1 attack for any ring homomorphic encryption scheme in the generic ring model. Moreover, when the plaintext space is a finite reduced black-box ring, we present a plaintext-recovery attack based on representation problem in black-box prime fields. In particular, if the ciphertext space has smooth characteristic, the plaintext-recovery attack is effectively computable in the generic ring model. [ABSTRACT FROM AUTHOR]

In a series of papers we develop symbolic computation algorithms to investigate finiteness of central configurations for the planar n -body problem. Our approach is based on Albouy-Kaloshin's work on finiteness of central configurations for the 5-body problems. In their paper, bicolored graphs called zw -diagrams were introduced for possible scenarios when the finiteness conjecture fails, and proving finiteness amounts to exclusions of central configurations associated to these diagrams. Following their method, the amount of computations becomes enormous when there are more than five bodies. In this paper we introduce matrix algebra for determination of possible diagrams and asymptotic orders, devise several criteria to reduce computational complexity, and determine possible zw -diagrams by automated deductions. For the planar six-body problem, we show that there are at most 86 zw -diagrams. [ABSTRACT FROM AUTHOR]

In this paper we analyze the computational costs of various operations and algorithms in algebraic number fields using exact arithmetic. Let K be an algebraic number field. In the first half of the paper, we calculate the running time and the size of the output of many operations in K in terms of the size of the input and the parameters of K. We include some earlier results about these, but we go further than them, e.g. we also analyze some R -specific operations in K like less-than comparison. In the second half of the paper, we analyze two algorithms: the Bareiss algorithm, which is an integer-preserving version of the Gaussian elimination, and the LLL algorithm, which is for lattice basis reduction. In both cases, we extend the algorithm from Z n to K n , and give a polynomial upper bound on the running time when the computations in K are performed exactly (as opposed to floating-point approximations). [ABSTRACT FROM AUTHOR]

This article develops new techniques to classify critical configurations for 3D scene reconstruction from images taken by unknown cameras. Generally, all information can be uniquely recovered if enough images and image points are provided, but there are certain cases where unique recovery is impossible; these are called critical configurations. In this paper, we use an algebraic approach to study the critical configurations for two projective cameras. We show that all critical configurations lie on quadric surfaces, and classify exactly which quadrics constitute a critical configuration. The paper also describes the relation between the different reconstructions when unique reconstruction is impossible. [ABSTRACT FROM AUTHOR]

In recent years there has been increased use of machine learning (ML) techniques within mathematics, including symbolic computation where it may be applied safely to optimise or select algorithms. This paper explores whether using explainable AI (XAI) techniques on such ML models can offer new insight for symbolic computation, inspiring new implementations within computer algebra systems that do not directly call upon AI tools. We present a case study on the use of ML to select the variable ordering for cylindrical algebraic decomposition. It has already been demonstrated that ML can make the choice well, but here we show how the SHAP tool for explainability can be used to inform new heuristics of a size and complexity similar to those human-designed heuristics currently commonly used in symbolic computation. [ABSTRACT FROM AUTHOR]

Satisfiability modulo theories (SMT) solvers check the satisfiability of quantifier-free first-order logic formulae over different theories. We consider the theory of non-linear real arithmetic where the formulae are logical combinations of polynomial constraints. Here a commonly used tool is the cylindrical algebraic decomposition (CAD) to decompose the real space into cells where the constraints are truth-invariant through the use of projection polynomials. A CAD encodes more information than necessary for checking satisfiability. One approach to address this is to repackage the CAD theory into a search-based algorithm: one that guesses sample points to satisfy the formula, and generalizes guesses that conflict constraints to cylindrical cells around samples which are avoided in the continuing search. This can lead to a satisfying assignment more quickly, or conclude unsatisfiability with far fewer cells. A notable example of this approach is Jovanović and de Moura's NLSAT algorithm. Since these cells are being produced locally to a sample there is scope to use fewer projection polynomials than the traditional CAD projection. The original NLSAT algorithm reduced the set a little; while Brown's single cell construction reduced it much further still. However, it refines a cell polynomial-by-polynomial, meaning the shape and size of the cell produced depends on the order in which the polynomials are considered. The present paper proposes a method to construct such cells levelwise , i.e. built level-by-level according to a variable ordering instead of polynomial-by-polynomial for all levels. We still use a reduced number of projection polynomials, but can now consider a variety of different reductions and use heuristics to select the projection polynomials in order to optimize the shape of the cell under construction. The new method can thus improve the performance of the NLSAT algorithm. We formulate all the necessary theory that underpins the algorithm as a proof system : while not a common presentation for work in this field, it is valuable in allowing an elegant decoupling of heuristic decisions from the main algorithm and its proof of correctness. We expect the symbolic computation community may find uses for it in other areas too. In particular, the proof system could be a step towards formal proofs for non-linear real arithmetic. This work has been implemented in the SMT-RAT solver and the benefits of the levelwise construction are validated experimentally on the SMT-LIB benchmark library. We also compare several heuristics for the construction and observe that each heuristic has strengths offering potential for further exploitation of the new approach. [ABSTRACT FROM AUTHOR]

Let G be a permutation group acting on a set Ω. Best known algorithms for computing the centralizer of G in the symmetric group on Ω are all based on the same general approach that involves solving the following two fundamental problems: given a G -orbit Δ of size n , compute the centralizer of the restriction of G to Δ in the symmetric group on Δ; and given two G -orbits Δ and Δ ′ each of size n , find an equivalence between the action of G restricted to Δ and the action of G restricted to Δ ′ when one exists. If G is given by a generating set X , previous solutions to each of these two problems take O (| X | n 2) time. In this paper, we first solve each fundamental problem in O (δ n + | X | n log ⁡ n) time, where δ is the depth of the breadth-first-search Schreier tree for X rooted at some fixed vertex. For the important class of small-base groups G , we improve the theoretical worst-case time bound of our solutions to O (n log c ⁡ n + | X | n log ⁡ n) for some constant c. Moreover, if ⌈ 20 log 2 ⁡ n ⌉ uniformly distributed random elements of G are given in advance, our solutions have, with probability at least 1 − 1 / n , a running time of O (n log 2 ⁡ n + | X | n log ⁡ n). We then apply our solutions to obtain a more efficient algorithm for computing the centralizer of G in the symmetric group on Ω. In an experimental evaluation we demonstrate that it is substantially faster than previously known algorithms. [ABSTRACT FROM AUTHOR]

We study linear PDE with constant coefficients. The constant rank condition on a system of linear PDEs with constant coefficients is often used in the theory of compensated compactness. While this is a purely linear algebraic condition, the nonlinear algebra concept of primary decomposition is another important tool for studying such system of PDEs. In this paper we investigate the connection between these two concepts. From the nonlinear analysis point of view, we make some progress in the study of weak lower semicontinuity of integral functionals defined on sequences of PDE constrained fields, when the PDEs do not have constant rank. [ABSTRACT FROM AUTHOR]

The Chow ring of the moduli space of stable marked curves is generated by Keel's divisor classes. The top graded part of this Chow ring is isomorphic to the integers, generated by the class of a single point. In this paper, we give an equivalent graphical characterization on the monomials in this Chow ring, as well as the characterization on the algebraic reduction on such monomials. Moreover, we provide an algorithm for computing the intersection degree of tuples of Keel's divisor classes — we call it the forest algorithm; the complexity of which is O (n 3) in the worst case, where n refers to the number of marks on the stable curves in the ambient moduli space. Last but not least, we introduce three identities on multinomial coefficients which naturally came into play, showing that they are all equivalent to the correctness of the base case of the forest algorithm. [ABSTRACT FROM AUTHOR]

For any zero-dimensional polynomial ideal I and any nonzero polynomial F , this paper shows that the union of the multi-set of zeros of the ideal sum I + 〈 F 〉 and that of the ideal quotient I : 〈 F 〉 is equal to the multi-set of zeros of I , where zeros are counted with multiplicities. Based on this zero relation and the computation of Gröbner bases, a complete multiplicity-preserved algorithm is proposed to decompose any zero-dimensional polynomial set into finitely many squarefree normal triangular sets, resulting in a squarefree normal representation for the zeros of the polynomial set. In the representation the multiplicities of the zeros of the triangular sets can be read out directly. Examples and experiments are presented to illustrate the algorithm and its performance. [ABSTRACT FROM AUTHOR]

In algebraic geometry, it is important to provide effective parametrizations for families of curves, both in theory and in practice. In this paper, we present such an effective parametrization for the moduli of genus-5 curves that are neither hyperelliptic nor trigonal. Subsequently, we construct an algorithm for a complete enumeration of non-special genus-5 curves having more rational points than a specified bound, where "non-special curve" means that the curve is non-hyperelliptic and non-trigonal with mild singularities of the associated sextic model that we propose. As a practical application, we implement this algorithm using the computer algebra system MAGMA, specifically for curves over the prime field of characteristic 3. [ABSTRACT FROM AUTHOR]

The MacWilliams' Extension Theorem is a classical result by Florence Jessie MacWilliams. It shows that every linear isometry between linear block-codes endowed with the Hamming distance can be extended to a linear isometry of the ambient space. Such an extension fails to exist in general for rank-metric codes, that is, one can easily find examples of linear isometries between rank-metric codes which cannot be extended to linear isometries of the ambient space. In this paper, we explore to what extent a MacWilliams' Extension Theorem may hold for rank-metric codes. We provide an extensive list of examples of obstructions to the existence of an extension, as well as a positive result. [ABSTRACT FROM AUTHOR]

In this paper we study the Shor relaxation of quadratic programs by fixing a feasible set and considering the space of objective functions for which the Shor relaxation is exact. We first give conditions under which this region is invariant under the choice of generators defining the feasible set. We then describe this region when the feasible set is invariant under the action of a subgroup of the general linear group. We conclude by applying these results to quadratic binary programs. We give an explicit description of objective functions where the Shor relaxation is exact and use this knowledge to design an algorithm that produces candidate solutions for binary quadratic programs. [ABSTRACT FROM AUTHOR]

The isotropic Grassmannian parametrizes isotropic subspaces of a vector space equipped with a quadratic form. In this paper, we show that any maximal isotropic Grassmannian in its Plücker embedding can be defined by pulling back the equations of G r iso (3 , 7) or G r iso (4 , 8). [ABSTRACT FROM AUTHOR]

We develop the first algorithm for computing rational solutions of scalar integro-differential equations with polynomial coefficients. It starts by finding the possible poles of a rational solution. Then, bounding the order of each pole and solving an algebraic linear system, we compute the singular part of rational solutions at each possible pole. Finally, using partial fraction decomposition, the polynomial part of rational solutions is obtained by computing polynomial solutions of a non-homogeneous scalar integro-differential equation with a polynomial right-hand side. The paper is illustrated by examples where the computations are done with our Maple implementation. [ABSTRACT FROM AUTHOR]

Recent work of Ablett (2021) and Kapustka, Kapustka, Ranestad, Schenck, Stillman and Yuan (2021) outlines a number of constructions for singular Gorenstein codimension four varieties. Earlier work of Coughlan, Gołȩbiowski, Kapustka and Kapustka (2016) details a series of nonsingular Gorenstein codimension four constructions with different Betti tables. In this paper we exhibit a number of flat deformations between Gorenstein codimension four varieties in the same Hilbert scheme, realising many of the singular varieties as specialisations of the earlier nonsingular varieties. [ABSTRACT FROM AUTHOR]

Mixed trigonometric-polynomials (MTPs) are functions of the form f (x , sin ⁡ x , cos ⁡ x) where f is a trivariate polynomial with rational coefficients, and the argument x ranges over the reals. In this paper, an algorithm "isolating" all the real roots of an MTP is provided and implemented. It automatically divides the real roots into two parts: one consists of finitely many roots in an interval [ μ − , μ + ] while the other consists of countably many roots in R ﹨ [ μ − , μ + ]. For the roots in [ μ − , μ + ] , the algorithm returns isolating intervals and corresponding multiplicities while for those greater than μ + , it returns finitely many mutually disjoint small intervals I i ⊂ [ − π , π ] , integers c i > 0 and multisets of root multiplicity { m j , i } j = 1 c i such that any root > μ + is in the set (∪ i ∪ k ∈ N (I i + 2 k π)) and any interval I i + 2 k π ⊂ (μ + , ∞) contains exactly c i distinct roots with multiplicities m 1 , i ,... , m c i , i , respectively. The efficiency of the algorithm is shown by experiments. The method used to isolate the roots in [ μ − , μ + ] is applicable to any other bounded interval as well. The algorithm takes advantages of the weak Fourier sequence technique and deals with the intervals period-by-period without scaling the coordinate so to keep the length of the sequence short. The new approaches can easily be modified to decide whether there is any root, or whether there are infinitely many roots in unbounded intervals of the form (− ∞ , a) or (a , ∞) with a ∈ Q. [ABSTRACT FROM AUTHOR]

Answering connectivity queries in real algebraic sets is a fundamental problem in effective real algebraic geometry that finds many applications in e.g. robotics where motion planning issues are topical. This computational problem is tackled through the computation of so-called roadmaps which are real algebraic subsets of the set V under study, of dimension at most one, and which have a connected intersection with all semi-algebraically connected components of V. Algorithms for computing roadmaps rely on statements establishing connectivity properties of some well-chosen subsets of V , assuming that V is bounded. In this paper, we extend such connectivity statements by dropping the boundedness assumption on V. This exploits properties of so-called generalized polar varieties , which are critical loci of V for some well-chosen polynomial maps. [ABSTRACT FROM AUTHOR]

The famous open Vizing conjecture claims that the domination number of the Cartesian product graph of two graphs G and H is at least the product of the domination numbers of G and H. Recently Gaar, Krenn, Margulies and Wiegele used the graph class G of all graphs with n G vertices and domination number k G and reformulated Vizing's conjecture as the problem that for all graph classes G and H the Vizing polynomial is sum-of-squares (SOS) modulo the Vizing ideal. By solving semidefinite programs (SDPs) and clever guessing they derived SOS-certificates for some values of k G , n G , k H , and n H. In this paper, we consider their approach for k G = k H = 1. For this case we are able to derive the unique reduced Gröbner basis of the Vizing ideal. Based on this, we deduce the minimum degree (n G + n H − 1) / 2 of an SOS-certificate for Vizing's conjecture, which is the first result of this kind. Furthermore, we present a method to find certificates for graph classes G and H with n G + n H − 1 = d for general d , which is again based on solving SDPs, but does not depend on guessing and depends on much smaller SDPs. We implement our new method in SageMath and give new SOS-certificates for all graph classes G and H with k G = k H = 1 and n G + n H ≤ 15. [ABSTRACT FROM AUTHOR]

Given n ⩾ 2 and a 1 , ... , a n ∈ N. Let S = 〈 a 1 , ... , a n 〉 be a semigroup. The aim of this paper is to give an effective pseudo-polynomial algorithm on a 1 , which computes the Apéry set and the Frobenius number of S. We also find the Gröbner basis of the toric ideal defined by S , for the weighted degree reverse lexicographical order ≺ w to x 1 , ... , x n , without using Buchberger's algorithm. As an application we introduce and study some special classes of semigroups. Namely, when S is generated by generalized arithmetic progressions and generalized almost arithmetic progressions with the ratio a positive or a negative number. We determine symmetric and almost symmetric semigroups generated by a generalized arithmetic progression. [ABSTRACT FROM AUTHOR]

The reduction of a multidimensional system is closely related to the reduction of a multivariate polynomial matrix, for which the Smith normal form of the matrix plays a key role. In this paper, we investigate the reduction of weakly linear multivariate polynomial matrices to their Smith normal forms. Using hierarchical-recursive method and Quillen-Suslin Theorem, we derive some necessary and sufficient conditions ensuring that such matrices can be reduced to their Smith normal forms, and these conditions are easily checked by computing the reduced Gröbner bases of the relevant polynomial ideals. Based on the new results, we propose an algorithm for reducing weakly linear multivariate polynomial matrices to their Smith normal forms. [ABSTRACT FROM AUTHOR]

We introduce an algorithmic approach based on a generating tree method for enumerating the inversion sequences with various pattern-avoidance restrictions. For a given set of patterns, we propose an algorithm that outputs either an accurate description of the succession rules of the corresponding generating tree or an ansatz. By using this approach, we determine the generating trees for the pattern classes I n (000 , 021) , I n (100 , 021) , I n (110 , 021) , I n (102 , 021) , I n (100 , 012) , I n (011 , 201) , I n (011 , 210) and I n (120 , 210). Then we use the kernel method, obtain generating functions of each class, and find enumerating formulas. Lin and Yan studied the classification of the Wilf-equivalences for inversion sequences avoiding pairs of length-three patterns and showed that there are 48 Wilf classes among 78 pairs. In this paper, we solve six open cases for such pattern classes. Moreover, we extend the algorithm to restricted growth sequences and apply it to several classes. In particular, we present explicit formulas for the generating functions of the restricted growth sequences that avoid either { 12313 , 12323 } , { 12313 , 12323 , 12333 } , or { 123 ⋯ ℓ 1 }. [ABSTRACT FROM AUTHOR]

In this paper, we propose a symbolic-numeric algorithm to count the number of solutions of a zero-dimensional square polynomial system within a local region. We show that the algorithm succeeds under the condition that the region is sufficiently small and well-isolating for a k -fold solution z of the system. In our analysis, we derive a bound on the size of the region that guarantees success. We further argue that this size depends on local parameters such as the norm and multiplicity of z as well as the distances between z and all other solutions. Efficiency of our method stems from the fact that we reduce the problem of counting the roots of the original system to the problem of solving a truncated system of degree k. In particular, if the multiplicity k of z is small compared to the total degrees of the original polynomials, our method considerably improves upon known complete and certified methods. We see a series of applications of our approach. When combined with a numerical solver in the fashion of an a posteriori certification step, we obtain a certified and reliable method for solving polynomial systems while profiting both from the efficiency of the numerical algorithm and the reliability of the symbolic approach. An alternative application results from incorporating our algorithm as inclusion predicate into an elimination method. For the special case of bivariate systems, we experimentally show that this approach leads to a significant improvement over an existing state-of-the-art elimination method. [ABSTRACT FROM AUTHOR]

A singular k -tuple of a tensor T of format (n 1 , ... , n k) is essentially a complex critical point of the distance function from T constrained to the cone of tensors of format (n 1 , ... , n k) of rank at most one. A generic tensor has finitely many complex singular k -tuples, and their number depends only on the tensor format. Furthermore, if we fix the first k − 1 dimensions n i , then the number of singular k -tuples of a generic tensor becomes a monotone non-decreasing function in one integer variable n k , that stabilizes when (n 1 , ... , n k) reaches a boundary format. In this paper, we study the linear span of singular k -tuples of a generic tensor. Its dimension also depends only on the tensor format. In particular, we concentrate on special order three tensors and order- k tensors of format (2 , ... , 2 , n). As a consequence, if again we fix the first k − 1 dimensions n i and let n k increase, we show that in these special formats, the dimension of the linear span stabilizes as well, but at some concise non-sub-boundary format. We conjecture that this phenomenon holds for an arbitrary format with k > 3. Finally, we provide equations for the linear span of singular triples of a generic order three tensor T of some special non-sub-boundary format. From these equations, we conclude that T belongs to the linear span of its singular triples, and we conjecture that this is the case for every tensor format. [ABSTRACT FROM AUTHOR]

In this study, we determine the structure of cyclic codes over the ring Z q [ u ] / 〈 u 2 〉 which is isomorphic to R = Z q + u Z q where q = p s , p is a prime, s is a positive integer, and u 2 = 0. This is equivalent to determining the algebraic structure of ideals of the polynomial quotient ring R [ x ] / 〈 x n − 1 〉 , which is addressed in this paper completely. By establishing the structure of ideals of R [ x ] / 〈 x n − 1 〉 with gcd ⁡ (p , n) = 1 , we present an exact formula that enumerates the number of ideals of this ring that leads to the enumeration of cyclic codes over this ring. Finally, we consider and explore some special families of cyclic codes for some specific q and determine their size. [ABSTRACT FROM AUTHOR]

The new identifiable case appeared in Angelini et al. (2018) , together with the analysis on simultaneous identifiability of pairs of ternary forms recently developed in Beorchia and Galuppi (2022) , suggested the following conjecture towards a complete classification of all simultaneous Waring identifiable cases: for any d ≥ 2 , the general polynomial vectors consisting of d − 1 ternary forms of degree d and a ternary form of degree d + 1 , with rank d 2 + d + 2 2 , are identifiable over C. In this paper, by means of a computer-aided procedure inspired to the one described in Angelini et al. (2018) , we obtain that the case d = 4 contradicts the previous conjecture, admitting at least 36 complex simultaneous Waring decompositions (of length 11) instead of 1. [ABSTRACT FROM AUTHOR]