-
Maximizing the Maximum Degree in Ordered Yao Graphs
Authors:
Péter Ágoston,
Adrian Dumitrescu,
Arsenii Sagdeev,
Karamjeet Singh,
Ji Zeng
Abstract:
For an ordered point set in a Euclidean space or, more generally, in an abstract metric space, the ordered Yao graph is obtained by connecting each of the points to its closest predecessor by a directed edge. We show that for every set of $n$ points in $\mathbb{R}^d$, there exists an order such that the corresponding ordered Yao graph has maximum degree at least $\log{n}/(4d)$. Apart from the…
▽ More
For an ordered point set in a Euclidean space or, more generally, in an abstract metric space, the ordered Yao graph is obtained by connecting each of the points to its closest predecessor by a directed edge. We show that for every set of $n$ points in $\mathbb{R}^d$, there exists an order such that the corresponding ordered Yao graph has maximum degree at least $\log{n}/(4d)$. Apart from the $1/(4d)$ factor, this bound is the best possible. As for the abstract setting, we show that for every $n$-element metric space, there exists an order such that the corresponding ordered Yao graph has maximum degree $Ω(\sqrt{\log{n}/\log\log{n}})$.
△ Less
Submitted 13 June, 2024;
originally announced June 2024.
-
On the structure of envy-free orientations on graphs
Authors:
Jinghan A Zeng,
Ruta Mehta
Abstract:
Fair division is the problem of allocating a set of items among agents in a fair manner. One of the most sought-after fairness notions is envy-freeness (EF), requiring that no agent envies another's allocation. When items are indivisible, it ceases to exist, and envy-freeness up to any good (EFX) emerged as one of its strongest relaxations. The existence of EFX allocations is arguably the biggest…
▽ More
Fair division is the problem of allocating a set of items among agents in a fair manner. One of the most sought-after fairness notions is envy-freeness (EF), requiring that no agent envies another's allocation. When items are indivisible, it ceases to exist, and envy-freeness up to any good (EFX) emerged as one of its strongest relaxations. The existence of EFX allocations is arguably the biggest open question within fair division. Recently, Christodoulou, Fiat, Koutsoupias, and Sgouritsa (EC 2023) introduced showed that EFX allocations exist for the case of graphical valuations where an instance is represented by a graph: nodes are agents, edges are goods, and each agent values only her incident edges. On the other hand, they showed NP-hardness for checking the existence of EFX orientation where every edge is allocated to one of its incident vertices, and asked for a characterization of graphs that exhibit EFX orientation regardless of the assigned valuations. In this paper, we make significant progress toward answering their question. We introduce the notion of strongly EFX orientable graphs -- graphs that have EFX orientations regardless of how much agents value the edges. We show a surprising connection between this property and the chromatic number of the graph, namely $χ(G)$ for graph $G$. In particular, we show that graphs with $χ(G)\le 2$ are strongly EFX orientable, and those with $χ(G)>3$ are not strongly EFX orientable. We provide examples of strongly EFX orientable and non-strongly EFX orientable graphs of $χ(G)=3$ to prove tightness. Finally, we give a complete characterization of strong EFX orientability when restricted to binary valuations.
△ Less
Submitted 21 April, 2024;
originally announced April 2024.
-
Gamma positivity of variations of $(α,t)$-Eulerian polynomials
Authors:
Chao Xu,
Jiang Zeng
Abstract:
In 1977 Carlitz and Scoville introduced the cycle $(α,t)$-Eulerian polynomials $A^{\mathrm{cyc}}_n(x,y, t\,|\,α)$ by enumerating permutations with respect to the number of excedances, drops, fixed points and cycles. In this paper, we introduce a nine-variable generalization of the Eulerian polynomials $A_n(u_1,u_2,u_3,u_4, f, g, t\,|\,α, β)$ in terms of descent based statistics of permutations and…
▽ More
In 1977 Carlitz and Scoville introduced the cycle $(α,t)$-Eulerian polynomials $A^{\mathrm{cyc}}_n(x,y, t\,|\,α)$ by enumerating permutations with respect to the number of excedances, drops, fixed points and cycles. In this paper, we introduce a nine-variable generalization of the Eulerian polynomials $A_n(u_1,u_2,u_3,u_4, f, g, t\,|\,α, β)$ in terms of descent based statistics of permutations and prove a connection formula between these two kinds of generalized Eulerian polynomials. By exploring the connection formula, we derive plainly the exponential generating function of the latter polynomials and various $γ$-positive formulas for variants of Eulerian polynomials. In particular, our results unify and strengthen the recent results by Ji and Ji-Lin. In related work to the transition matrix between the Specht and web bases, Hwang, Jang and Oh recently introduced the web permutations, which can be characterised by cycle André permutations. We show that enumerating the latter permutations with respect to the number of drops, fixed points and cycles gives rise to the normalised $γ$-vectors of the $(α,t)$-Eulerian polynomials. Our result generalizes and unifies several known results in the literature.
△ Less
Submitted 25 June, 2024; v1 submitted 12 April, 2024;
originally announced April 2024.
-
Mahonian-Stirling statistics for partial permutations
Authors:
Ming-Jian Ding,
Jiang Zeng
Abstract:
Recently Cheng et al. (Adv. in Appl. Math. 143 (2023) 102451) generalized the inversion number to partial permutations, which are also known as Laguerre digraphs, and asked for a suitable analogue of MacMahon's major index. We provide such a major index, namely, the corresponding maj and inv statistics are equidistributed, and exhibit a Haglund-Remmel-Wilson type identity. We then interpret some J…
▽ More
Recently Cheng et al. (Adv. in Appl. Math. 143 (2023) 102451) generalized the inversion number to partial permutations, which are also known as Laguerre digraphs, and asked for a suitable analogue of MacMahon's major index. We provide such a major index, namely, the corresponding maj and inv statistics are equidistributed, and exhibit a Haglund-Remmel-Wilson type identity. We then interpret some Jacobi-Rogers polynomials in terms of Laguerre digraphs generalizing Deb and Sokal's alternating Laguerre digraph interpretation of some special Jacobi-Rogers polynomials.
△ Less
Submitted 1 April, 2024;
originally announced April 2024.
-
DG singular equivalence and singular locus
Authors:
Leilei Liu,
Jieheng Zeng
Abstract:
For a commutative Gorenstein Noetherian ring $R$, we construct an affine scheme $X$ solely from DG singularity category $S_{dg}(R)$ of $R$ such that there is a finite surjective morphism $X \rightarrow \mathrm{Spec}(R /I)$, where $\mathrm{Spec}(R /I)$ is the singular locus in $\mathrm{Spec}(R)$. As an application, for two such rings with equivalent DG singularity categories, we prove that the sing…
▽ More
For a commutative Gorenstein Noetherian ring $R$, we construct an affine scheme $X$ solely from DG singularity category $S_{dg}(R)$ of $R$ such that there is a finite surjective morphism $X \rightarrow \mathrm{Spec}(R /I)$, where $\mathrm{Spec}(R /I)$ is the singular locus in $\mathrm{Spec}(R)$. As an application, for two such rings with equivalent DG singularity categories, we prove that the singular loci in their affine schemes have the same dimension.
△ Less
Submitted 31 March, 2024; v1 submitted 20 March, 2024;
originally announced March 2024.
-
Paintability of $r$-chromatic graphs
Authors:
Peter Bradshaw,
Jinghan A Zeng
Abstract:
The paintability of a graph is a coloring parameter defined in terms of an online list coloring game. In this paper we ask, what is the paintability of a graph $G$ of maximum degree $Δ$ and chromatic number $r$? By considering the Alon-Tarsi number of $G$, we prove that the paintability of $G$ is at most $\left(1 - \frac{1}{4r+1} \right ) Δ+ 2$. We also consider the DP-paintability of $G$, which i…
▽ More
The paintability of a graph is a coloring parameter defined in terms of an online list coloring game. In this paper we ask, what is the paintability of a graph $G$ of maximum degree $Δ$ and chromatic number $r$? By considering the Alon-Tarsi number of $G$, we prove that the paintability of $G$ is at most $\left(1 - \frac{1}{4r+1} \right ) Δ+ 2$. We also consider the DP-paintability of $G$, which is defined in terms of an online DP-coloring game. By considering the strict type $3$ degeneracy parameter recently introduced by Zhou, Zhu, and Zhu, we show that when $r$ is fixed and $Δ$ is sufficiently large, the DP-paintability of $G$ is at most $Δ- Ω( \sqrt{Δ\log Δ})$.
△ Less
Submitted 4 July, 2024; v1 submitted 18 March, 2024;
originally announced March 2024.
-
Saturation results around the Erdős--Szekeres problem
Authors:
Gábor Damásdi,
Zichao Dong,
Manfred Scheucher,
Ji Zeng
Abstract:
In this paper, we consider saturation problems related to the celebrated Erdős--Szekeres convex polygon problem. For each $n \ge 7$, we construct a planar point set of size $(7/8) \cdot 2^{n-2}$ which is saturated for convex $n$-gons. That is, the set contains no $n$ points in convex position while the addition of any new point creates such a configuration. This demonstrates that the saturation nu…
▽ More
In this paper, we consider saturation problems related to the celebrated Erdős--Szekeres convex polygon problem. For each $n \ge 7$, we construct a planar point set of size $(7/8) \cdot 2^{n-2}$ which is saturated for convex $n$-gons. That is, the set contains no $n$ points in convex position while the addition of any new point creates such a configuration. This demonstrates that the saturation number is smaller than the Ramsey number for the Erdős--Szekeres problem. The proof also shows that the original Erdős--Szekeres construction is indeed saturated.
Our construction is based on a similar improvement for the saturation version of the cups-versus-caps theorem. Moreover, we consider the generalization of the cups-versus-caps theorem to monotone paths in ordered hypergraphs. In contrast to the geometric setting, we show that this abstract saturation number is always equal to the corresponding Ramsey number.
△ Less
Submitted 2 December, 2023;
originally announced December 2023.
-
On cliques in three-dimensional dense point-line arrangements
Authors:
Andrew Suk,
Ji Zeng
Abstract:
As a variant of the celebrated Szemerédi--Trotter theorem, Guth and Katz proved that $m$ points and $n$ lines in $\mathbb{R}^3$ with at most $\sqrt{n}$ lines in a common plane must determine at most $O(m^{1/2}n^{3/4})$ incidences for $n^{1/2}\leq m\leq n^{3/2}$. This upper bound is asymptotically tight and has an important application in Erd{\H o}s distinct distance problem. We characterize the ex…
▽ More
As a variant of the celebrated Szemerédi--Trotter theorem, Guth and Katz proved that $m$ points and $n$ lines in $\mathbb{R}^3$ with at most $\sqrt{n}$ lines in a common plane must determine at most $O(m^{1/2}n^{3/4})$ incidences for $n^{1/2}\leq m\leq n^{3/2}$. This upper bound is asymptotically tight and has an important application in Erd{\H o}s distinct distance problem. We characterize the extremal constructions towards the Guth--Katz bound by proving that such a large dense point-line arrangement must contain a $k$-clique in general position provided $m \ll n$. This is an analog of a result by Solymosi for extremal Szemerédi--Trotter constructions in the plane.
△ Less
Submitted 8 November, 2023;
originally announced November 2023.
-
Random Turán and counting results for general position sets over finite fields
Authors:
Yaobin Chen,
Xizhi Liu,
Jiaxi Nie,
Ji Zeng
Abstract:
Let $α(\mathbb{F}_q^d,p)$ denote the maximum size of a general position set in a $p$-random subset of $\mathbb{F}_q^d$. We determine the order of magnitude of $α(\mathbb{F}_q^2,p)$ up to polylogarithmic factors for all possible values of $p$, improving the previous results obtained by Roche-Newton--Warren and Bhowmick--Roche-Newton. For $d \ge 3$ we prove upper bounds for $α(\mathbb{F}_q^d,p)$ tha…
▽ More
Let $α(\mathbb{F}_q^d,p)$ denote the maximum size of a general position set in a $p$-random subset of $\mathbb{F}_q^d$. We determine the order of magnitude of $α(\mathbb{F}_q^2,p)$ up to polylogarithmic factors for all possible values of $p$, improving the previous results obtained by Roche-Newton--Warren and Bhowmick--Roche-Newton. For $d \ge 3$ we prove upper bounds for $α(\mathbb{F}_q^d,p)$ that are essentially tight within certain ranges for $p$.
We establish the upper bound $2^{(1+o(1))q}$ for the number of general position sets in $\mathbb{F}_q^d$, which matches the trivial lower bound $2^{q}$ asymptotically in the exponent. We also refine this counting result by proving an asymptotically tight (in the exponent) upper bound for the number of general position sets with a fixed size. The latter result for $d=2$ improves a result of Roche-Newton--Warren.
Our proofs are grounded in the hypergraph container method, and additionally, for $d=2$ we also leverage the pseudorandomness of the point-line incidence graph of $\mathbb{F}_{q}^2$.
△ Less
Submitted 25 February, 2024; v1 submitted 14 September, 2023;
originally announced September 2023.
-
Real-rootedness of the type A minuscule polynomials
Authors:
Ming-Jian Ding,
Jiang Zeng
Abstract:
We prove two recent conjectures of Bourn and Erickson (2023) regarding the real-rootedness of a certain family of polynomials $N_n(t)$ as well as the sum of their coefficients. These polynomials arise as the numerators of generating functions in the context of the discrete one-dimensional earth mover's distance (EMD) and have also connection to the Wiener index of minuscule lattices. We also prove…
▽ More
We prove two recent conjectures of Bourn and Erickson (2023) regarding the real-rootedness of a certain family of polynomials $N_n(t)$ as well as the sum of their coefficients. These polynomials arise as the numerators of generating functions in the context of the discrete one-dimensional earth mover's distance (EMD) and have also connection to the Wiener index of minuscule lattices. We also prove that the coefficients of $N_n(x)$ are asymptotically normal, the coefficient matrix of $N_n(x)$ is totally positive and the polynomial sequence $N_n(x)$'s is $x$-log-concave.
△ Less
Submitted 6 July, 2024; v1 submitted 31 August, 2023;
originally announced August 2023.
-
Note on disjoint faces in simple topological graphs
Authors:
Ji Zeng
Abstract:
We prove that every $n$-vertex complete simple topological graph generates at least $Ω(n)$ pairwise disjoint $4$-faces. This improves upon a recent result by Hubard and Suk. As an immediate corollary, every $n$-vertex complete simple topological graph drawn in the unit square generates a $4$-face with area at most $O(1/n)$. This can be seen as a topological variant of Heilbronn's problem for $4$-f…
▽ More
We prove that every $n$-vertex complete simple topological graph generates at least $Ω(n)$ pairwise disjoint $4$-faces. This improves upon a recent result by Hubard and Suk. As an immediate corollary, every $n$-vertex complete simple topological graph drawn in the unit square generates a $4$-face with area at most $O(1/n)$. This can be seen as a topological variant of Heilbronn's problem for $4$-faces. We construct examples showing that our result is asymptotically tight. We also discuss the similar problem for $k$-faces with arbitrary $k\geq 3$.
△ Less
Submitted 17 September, 2023; v1 submitted 9 August, 2023;
originally announced August 2023.
-
Some identities involving $q$-Stirling numbers of the second kind in type B
Authors:
Ming-Jian Ding,
Jiang Zeng
Abstract:
The recent interest in $q$-Stirling numbers of the second kind in type B prompted us to give a type B analogue of a classical identity connecting the $q$-Stirling numbers of the second kind and Carlitz's major $q$-Eulerian numbers, which turns out to be a $q$-analogue of an identity due to Bagno, Biagioli and Garber. We provide a combinatorial proof of this identity and an analytical proof of a mo…
▽ More
The recent interest in $q$-Stirling numbers of the second kind in type B prompted us to give a type B analogue of a classical identity connecting the $q$-Stirling numbers of the second kind and Carlitz's major $q$-Eulerian numbers, which turns out to be a $q$-analogue of an identity due to Bagno, Biagioli and Garber. We provide a combinatorial proof of this identity and an analytical proof of a more general identity for colored permutations. In addition, we prove some $q$-identities about the $q$-Stirling numbers of the second kind in types A, B and D.
△ Less
Submitted 12 January, 2024; v1 submitted 2 July, 2023;
originally announced July 2023.
-
Proof of an explicit formula for a series from Ramanujan's Notebooks via tree functions
Authors:
Ming-Jian Ding,
Jiang Zeng
Abstract:
We prove a recent conjecture, due to Vigren and Dieckmann, about an explicit triple sum formula for a series from Ramanujan's Notebooks. We shall give two proofs: the first one is by evaluation and based on the identity \begin{equation*}
\sum_{k=0}^\infty \frac{(x+k)^{m+k}}{k!}e^{-u(x+k)} u^k
= \sum_{j=0}^\infty \sum_{i=0}^{m}\binom{m+j}{i}
\stirl{m+j-i}{j}x^iu^j, \end{equation*} where…
▽ More
We prove a recent conjecture, due to Vigren and Dieckmann, about an explicit triple sum formula for a series from Ramanujan's Notebooks. We shall give two proofs: the first one is by evaluation and based on the identity \begin{equation*}
\sum_{k=0}^\infty \frac{(x+k)^{m+k}}{k!}e^{-u(x+k)} u^k
= \sum_{j=0}^\infty \sum_{i=0}^{m}\binom{m+j}{i}
\stirl{m+j-i}{j}x^iu^j, \end{equation*} where $\genfrac\{\}{0pt}{}{n}{k}$ is a Stirling number of the second kind, and the second one is combinatorial in nature and by induction.
△ Less
Submitted 2 July, 2023;
originally announced July 2023.
-
Nonsmooth Control Barrier Functions for Obstacle Avoidance between Convex Regions
Authors:
Akshay Thirugnanam,
Jun Zeng,
Koushil Sreenath
Abstract:
In this paper, we focus on non-conservative obstacle avoidance between robots with control affine dynamics with strictly convex and polytopic shapes. The core challenge for this obstacle avoidance problem is that the minimum distance between strictly convex regions or polytopes is generally implicit and non-smooth, such that distance constraints cannot be enforced directly in the optimization prob…
▽ More
In this paper, we focus on non-conservative obstacle avoidance between robots with control affine dynamics with strictly convex and polytopic shapes. The core challenge for this obstacle avoidance problem is that the minimum distance between strictly convex regions or polytopes is generally implicit and non-smooth, such that distance constraints cannot be enforced directly in the optimization problem. To handle this challenge, we employ non-smooth control barrier functions to reformulate the avoidance problem in the dual space, with the positivity of the minimum distance between robots equivalently expressed using a quadratic program. Our approach is proven to guarantee system safety. We theoretically analyze the smoothness properties of the minimum distance quadratic program and its KKT conditions. We validate our approach by demonstrating computationally-efficient obstacle avoidance for multi-agent robotic systems with strictly convex and polytopic shapes. To our best knowledge, this is the first time a real-time QP problem can be formulated for general non-conservative avoidance between strictly convex shapes and polytopes.
△ Less
Submitted 22 June, 2023;
originally announced June 2023.
-
On Induced Subgraph of Cartesian Product of Paths
Authors:
Jiasheng Zeng,
Xinmin Hou
Abstract:
Chung, Füredi, Graham, and Seymour (JCTA, 1988) constructed an induced subgraph of the hypercube $Q^n$ with $α(Q^n)+1$ vertices and with maximum degree smaller than $\lceil \sqrt{n} \rceil$. Subsequently, Huang (Annals of Mathematics, 2019) proved the Sensitivity Conjecture by demonstrating that the maximum degree of such an induced subgraph of hypercube $Q^n$ is at least $\lceil \sqrt{n} \rceil$,…
▽ More
Chung, Füredi, Graham, and Seymour (JCTA, 1988) constructed an induced subgraph of the hypercube $Q^n$ with $α(Q^n)+1$ vertices and with maximum degree smaller than $\lceil \sqrt{n} \rceil$. Subsequently, Huang (Annals of Mathematics, 2019) proved the Sensitivity Conjecture by demonstrating that the maximum degree of such an induced subgraph of hypercube $Q^n$ is at least $\lceil \sqrt{n} \rceil$, and posed the question: Given a graph $G$, let $f(G)$ be the minimum of the maximum degree of an induced subgraph of $G$ on $α(G)+1$ vertices, what can we say about $f(G)$? In this paper, we investigate this question for Cartesian product of paths $P_m$, denoted by $P_m^k$. We determine the exact values of $f(P_{m}^k)$ when $m=2n+1$ by showing that $f(P_{2n+1}^k)=1$ for $n\geq 2$ and $f(P_3^k)=2$, and give a nontrivial lower bound of $f(P_{m}^k)$ when $m=2n$ by showing that $f(P_{2n}^k)\geq \lceil \sqrt{β_nk}\rceil$. In particular, when $n=1$, we have $f(Q^k)=f(P_{2}^k)\ge \sqrt{k}$, which is Huang's result. The lower bounds of $f(P_{3}^k)$ and $f(P_{2n}^k)$ are given by using the spectral method provided by Huang.
△ Less
Submitted 6 June, 2023;
originally announced June 2023.
-
Velocity Obstacle for Polytopic Collision Avoidance for Distributed Multi-robot Systems
Authors:
Jihao Huang,
Jun Zeng,
Xuemin Chi,
Koushil Sreenath,
Zhitao Liu,
Hongye Su
Abstract:
Obstacle avoidance for multi-robot navigation with polytopic shapes is challenging. Existing works simplify the system dynamics or consider it as a convex or non-convex optimization problem with positive distance constraints between robots, which limits real-time performance and scalability. Additionally, generating collision-free behavior for polytopic-shaped robots is harder due to implicit and…
▽ More
Obstacle avoidance for multi-robot navigation with polytopic shapes is challenging. Existing works simplify the system dynamics or consider it as a convex or non-convex optimization problem with positive distance constraints between robots, which limits real-time performance and scalability. Additionally, generating collision-free behavior for polytopic-shaped robots is harder due to implicit and non-differentiable distance functions between polytopes. In this paper, we extend the concept of velocity obstacle (VO) principle for polytopic-shaped robots and propose a novel approach to construct the VO in the function of vertex coordinates and other robot's states. Compared with existing work about obstacle avoidance between polytopic-shaped robots, our approach is much more computationally efficient as the proposed approach for construction of VO between polytopes is optimization-free. Based on VO representation for polytopic shapes, we later propose a navigation approach for distributed multi-robot systems. We validate our proposed VO representation and navigation approach in multiple challenging scenarios including large-scale randomized tests, and our approach outperforms the state of art in many evaluation metrics, including completion rate, deadlock rate, and the average travel distance.
△ Less
Submitted 10 June, 2024; v1 submitted 16 April, 2023;
originally announced April 2023.
-
i2LQR: Iterative LQR for Iterative Tasks in Dynamic Environments
Authors:
Yifan Zeng,
Suiyi He,
Han Hoang Nguyen,
Yihan Li,
Zhongyu Li,
Koushil Sreenath,
Jun Zeng
Abstract:
This work introduces a novel control strategy called Iterative Linear Quadratic Regulator for Iterative Tasks (i2LQR), which aims to improve closed-loop performance with local trajectory optimization for iterative tasks in a dynamic environment. The proposed algorithm is reference-free and utilizes historical data from previous iterations to enhance the performance of the autonomous system. Unlike…
▽ More
This work introduces a novel control strategy called Iterative Linear Quadratic Regulator for Iterative Tasks (i2LQR), which aims to improve closed-loop performance with local trajectory optimization for iterative tasks in a dynamic environment. The proposed algorithm is reference-free and utilizes historical data from previous iterations to enhance the performance of the autonomous system. Unlike existing algorithms, the i2LQR computes the optimal solution in an iterative manner at each timestamp, rendering it well-suited for iterative tasks with changing constraints at different iterations. To evaluate the performance of the proposed algorithm, we conduct numerical simulations for an iterative task aimed at minimizing completion time. The results show that i2LQR achieves an optimized performance with respect to learning-based MPC (LMPC) as the benchmark in static environments, and outperforms LMPC in dynamic environments with both static and dynamics obstacles.
△ Less
Submitted 6 September, 2023; v1 submitted 27 February, 2023;
originally announced February 2023.
-
Learning Cut Selection for Mixed-Integer Linear Programming via Hierarchical Sequence Model
Authors:
Zhihai Wang,
Xijun Li,
Jie Wang,
Yufei Kuang,
Mingxuan Yuan,
Jia Zeng,
Yongdong Zhang,
Feng Wu
Abstract:
Cutting planes (cuts) are important for solving mixed-integer linear programs (MILPs), which formulate a wide range of important real-world applications. Cut selection -- which aims to select a proper subset of the candidate cuts to improve the efficiency of solving MILPs -- heavily depends on (P1) which cuts should be preferred, and (P2) how many cuts should be selected. Although many modern MILP…
▽ More
Cutting planes (cuts) are important for solving mixed-integer linear programs (MILPs), which formulate a wide range of important real-world applications. Cut selection -- which aims to select a proper subset of the candidate cuts to improve the efficiency of solving MILPs -- heavily depends on (P1) which cuts should be preferred, and (P2) how many cuts should be selected. Although many modern MILP solvers tackle (P1)-(P2) by manually designed heuristics, machine learning offers a promising approach to learn more effective heuristics from MILPs collected from specific applications. However, many existing learning-based methods focus on learning which cuts should be preferred, neglecting the importance of learning the number of cuts that should be selected. Moreover, we observe from extensive empirical results that (P3) what order of selected cuts should be preferred has a significant impact on the efficiency of solving MILPs as well. To address this challenge, we propose a novel hierarchical sequence model (HEM) to learn cut selection policies via reinforcement learning. Specifically, HEM consists of a two-level model: (1) a higher-level model to learn the number of cuts that should be selected, (2) and a lower-level model -- that formulates the cut selection task as a sequence to sequence learning problem -- to learn policies selecting an ordered subset with the size determined by the higher-level model. To the best of our knowledge, HEM is the first method that can tackle (P1)-(P3) in cut selection simultaneously from a data-driven perspective. Experiments show that HEM significantly improves the efficiency of solving MILPs compared to human-designed and learning-based baselines on both synthetic and large-scale real-world MILPs, including MIPLIB 2017. Moreover, experiments demonstrate that HEM well generalizes to MILPs that are significantly larger than those seen during training.
△ Less
Submitted 31 January, 2023;
originally announced February 2023.
-
On higher dimensional point sets in general position
Authors:
Andrew Suk,
Ji Zeng
Abstract:
A finite point set in $\mathbb{R}^d$ is in general position if no $d + 1$ points lie on a common hyperplane. Let $α_d(N)$ be the largest integer such that any set of $N$ points in $\mathbb{R}^d$ with no $d + 2$ members on a common hyperplane, contains a subset of size $α_d(N)$ in general position. Using the method of hypergraph containers, Balogh and Solymosi showed that $α_2(N) < N^{5/6 + o(1)}$.…
▽ More
A finite point set in $\mathbb{R}^d$ is in general position if no $d + 1$ points lie on a common hyperplane. Let $α_d(N)$ be the largest integer such that any set of $N$ points in $\mathbb{R}^d$ with no $d + 2$ members on a common hyperplane, contains a subset of size $α_d(N)$ in general position. Using the method of hypergraph containers, Balogh and Solymosi showed that $α_2(N) < N^{5/6 + o(1)}$. In this paper, we also use the container method to obtain new upper bounds for $α_d(N)$ when $d \geq 3$. More precisely, we show that if $d$ is odd, then $α_d(N) < N^{\frac{1}{2} + \frac{1}{2d} + o(1)}$, and if $d$ is even, we have $α_d(N) < N^{\frac{1}{2} + \frac{1}{d-1} + o(1)}$.
We also study the classical problem of determining the maximum number $a(d,k,n)$ of points selected from the grid $[n]^d$ such that no $k + 2$ members lie on a $k$-flat. For fixed $d$ and $k$, we show that \begin{equation*}
a(d,k,n)\leq O\left(n^{\frac{d}{2\lfloor (k+2)/4\rfloor}(1-\frac{1}{2\lfloor(k+2)/4\rfloor d+1})}\right), \end{equation*} which improves the previously best known bound of $O\left(n^{\frac{d}{\lfloor (k + 2)/2\rfloor}}\right)$ due to Lefmann when $k+2$ is congruent to 0 or 1 mod 4.
△ Less
Submitted 14 March, 2023; v1 submitted 29 November, 2022;
originally announced November 2022.
-
Iterative Convex Optimization for Model Predictive Control with Discrete-Time High-Order Control Barrier Functions
Authors:
Shuo Liu,
Jun Zeng,
Koushil Sreenath,
Calin A. Belta
Abstract:
Safety is one of the fundamental challenges in control theory. Recently, multi-step optimal control problems for discrete-time dynamical systems were formulated to enforce stability, while subject to input constraints as well as safety-critical requirements using discrete-time control barrier functions within a model predictive control (MPC) framework. Existing work usually focus on the feasibilit…
▽ More
Safety is one of the fundamental challenges in control theory. Recently, multi-step optimal control problems for discrete-time dynamical systems were formulated to enforce stability, while subject to input constraints as well as safety-critical requirements using discrete-time control barrier functions within a model predictive control (MPC) framework. Existing work usually focus on the feasibility or the safety for the optimization problem, and the majority of the existing work restrict the discussions to relative-degree one control barrier functions. Additionally, the real-time computation is challenging when a large horizon is considered in the MPC problem for relative-degree one or high-order control barrier functions. In this paper, we propose a framework that solves the safety-critical MPC problem in an iterative optimization, which is applicable for any relative-degree control barrier functions. In the proposed formulation, the nonlinear system dynamics as well as the safety constraints modeled as discrete-time high-order control barrier functions (DHOCBF) are linearized at each time step. Our formulation is generally valid for any control barrier function with an arbitrary relative-degree. The advantages of fast computational performance with safety guarantee are analyzed and validated with numerical results.
△ Less
Submitted 13 July, 2023; v1 submitted 9 October, 2022;
originally announced October 2022.
-
Enumeration of permutations by the parity of descent positions
Authors:
Qiongqiong Pan,
Jiang Zeng
Abstract:
Noticing that some recent variations of descent polynomials are special cases of Carlitz and Scoville's four-variable polynomials, which enumerate permutations by the parity of descent and ascent positions, we prove a $q$-analogue of Carlitz-Scoville's generating function by counting the inversion number and a type B analogue by enumerating the signed permutations with respect to the parity of des…
▽ More
Noticing that some recent variations of descent polynomials are special cases of Carlitz and Scoville's four-variable polynomials, which enumerate permutations by the parity of descent and ascent positions, we prove a $q$-analogue of Carlitz-Scoville's generating function by counting the inversion number and a type B analogue by enumerating the signed permutations with respect to the parity of desecnt and ascent positions. As a by-product of our formulas, we obtain a $q$-analogue of Chebikin's formula for alternating descent polynomials, an alternative proof of Sun's gamma-positivity of her bivariate Eulerian polynomials and a type B analogue, the latter refines Petersen's gamma-positivity of the type B Eulerian polynomials.
△ Less
Submitted 13 June, 2023; v1 submitted 30 September, 2022;
originally announced September 2022.
-
On asymptotic packing of convex geometric and ordered graphs
Authors:
Jiaxi Nie,
Erlang Surya,
Ji Zeng
Abstract:
A convex geometric graph $G$ is said to be packable if there exist edge-disjoint copies of $G$ in the complete convex geometric graph $K_n$ covering all but $o(n^2)$ edges. We prove that every convex geometric graph with cyclic chromatic number at most $4$ is packable. With a similar definition of packability for ordered graphs, we prove that every ordered graph with interval chromatic number at m…
▽ More
A convex geometric graph $G$ is said to be packable if there exist edge-disjoint copies of $G$ in the complete convex geometric graph $K_n$ covering all but $o(n^2)$ edges. We prove that every convex geometric graph with cyclic chromatic number at most $4$ is packable. With a similar definition of packability for ordered graphs, we prove that every ordered graph with interval chromatic number at most $3$ is packable. Arguments based on the average length of edges imply these results are best possible. We also identify a class of convex geometric graphs that are packable due to having many "long" edges.
△ Less
Submitted 23 July, 2022;
originally announced July 2022.
-
Semi-restricted Rock, Paper, Scissors
Authors:
Sam Spiro,
Erlang Surya,
Ji Zeng
Abstract:
Consider the following variant of Rock, Paper, Scissors (RPS) played by two players Rei and Norman. The game consists of $3n$ rounds of RPS, with the twist being that Rei (the restricted player) must use each of Rock, Paper, and Scissors exactly $n$ times during the $3n$ rounds, while Norman is allowed to play normally without any restrictions. Answering a question of Spiro, we show that a certain…
▽ More
Consider the following variant of Rock, Paper, Scissors (RPS) played by two players Rei and Norman. The game consists of $3n$ rounds of RPS, with the twist being that Rei (the restricted player) must use each of Rock, Paper, and Scissors exactly $n$ times during the $3n$ rounds, while Norman is allowed to play normally without any restrictions. Answering a question of Spiro, we show that a certain greedy strategy is the unique optimal strategy for Rei in this game, and that Norman's expected score is $Θ(\sqrt{n})$. Moreover, we study semi-restricted versions of general zero sum games and prove a number of results concerning their optimal strategies and expected scores, which in particular implies our results for semi-restricted RPS.
△ Less
Submitted 22 July, 2022;
originally announced July 2022.
-
Counting derangements with signed right-to-left minima and excedances
Authors:
Yanni Pei,
Jiang Zeng
Abstract:
Recently Alexandersson and Getachew proved some multivariate generalizations of a formula for enumerating signed excedances in derangements. In this paper we first relate their work to a recent continued fraction for permutations and confirm some of their observations. Our second main result is two refinements of their multivariate identities, which clearly explain the meaning of each term in thei…
▽ More
Recently Alexandersson and Getachew proved some multivariate generalizations of a formula for enumerating signed excedances in derangements. In this paper we first relate their work to a recent continued fraction for permutations and confirm some of their observations. Our second main result is two refinements of their multivariate identities, which clearly explain the meaning of each term in their main formulas.
We also explore some similar formulas for permutations of type B.
△ Less
Submitted 8 October, 2023; v1 submitted 22 June, 2022;
originally announced June 2022.
-
Unavoidable patterns in complete simple topological graphs
Authors:
Andrew Suk,
Ji Zeng
Abstract:
We show that every complete $n$-vertex simple topological graph contains a topological subgraph on at least $(\log n)^{1/4 - o(1)}$ vertices that is weakly isomorphic to the complete convex geometric graph or the complete twisted graph. This is the first improvement on the bound $Ω(\log^{1/8}n)$ obtained in 2003 by Pach, Solymosi, and Tóth. We also show that every complete $n$-vertex simple topolo…
▽ More
We show that every complete $n$-vertex simple topological graph contains a topological subgraph on at least $(\log n)^{1/4 - o(1)}$ vertices that is weakly isomorphic to the complete convex geometric graph or the complete twisted graph. This is the first improvement on the bound $Ω(\log^{1/8}n)$ obtained in 2003 by Pach, Solymosi, and Tóth. We also show that every complete $n$-vertex simple topological graph contains a plane path of length at least $(\log n)^{1 -o(1)}$.
△ Less
Submitted 6 September, 2022; v1 submitted 8 April, 2022;
originally announced April 2022.
-
Learning to Reformulate for Linear Programming
Authors:
Xijun Li,
Qingyu Qu,
Fangzhou Zhu,
Jia Zeng,
Mingxuan Yuan,
Kun Mao,
Jie Wang
Abstract:
It has been verified that the linear programming (LP) is able to formulate many real-life optimization problems, which can obtain the optimum by resorting to corresponding solvers such as OptVerse, Gurobi and CPLEX. In the past decades, a serial of traditional operation research algorithms have been proposed to obtain the optimum of a given LP in a fewer solving time. Recently, there is a trend of…
▽ More
It has been verified that the linear programming (LP) is able to formulate many real-life optimization problems, which can obtain the optimum by resorting to corresponding solvers such as OptVerse, Gurobi and CPLEX. In the past decades, a serial of traditional operation research algorithms have been proposed to obtain the optimum of a given LP in a fewer solving time. Recently, there is a trend of using machine learning (ML) techniques to improve the performance of above solvers. However, almost no previous work takes advantage of ML techniques to improve the performance of solver from the front end, i.e., the modeling (or formulation). In this paper, we are the first to propose a reinforcement learning-based reformulation method for LP to improve the performance of solving process. Using an open-source solver COIN-OR LP (CLP) as an environment, we implement the proposed method over two public research LP datasets and one large-scale LP dataset collected from practical production planning scenario. The evaluation results suggest that the proposed method can effectively reduce both the solving iteration number ($25\%\downarrow$) and the solving time ($15\%\downarrow$) over above datasets in average, compared to directly solving the original LP instances.
△ Less
Submitted 16 January, 2022;
originally announced January 2022.
-
An Improved Reinforcement Learning Algorithm for Learning to Branch
Authors:
Qingyu Qu,
Xijun Li,
Yunfan Zhou,
Jia Zeng,
Mingxuan Yuan,
Jie Wang,
Jinhu Lv,
Kexin Liu,
Kun Mao
Abstract:
Most combinatorial optimization problems can be formulated as mixed integer linear programming (MILP), in which branch-and-bound (B\&B) is a general and widely used method. Recently, learning to branch has become a hot research topic in the intersection of machine learning and combinatorial optimization. In this paper, we propose a novel reinforcement learning-based B\&B algorithm. Similar to offl…
▽ More
Most combinatorial optimization problems can be formulated as mixed integer linear programming (MILP), in which branch-and-bound (B\&B) is a general and widely used method. Recently, learning to branch has become a hot research topic in the intersection of machine learning and combinatorial optimization. In this paper, we propose a novel reinforcement learning-based B\&B algorithm. Similar to offline reinforcement learning, we initially train on the demonstration data to accelerate learning massively. With the improvement of the training effect, the agent starts to interact with the environment with its learned policy gradually. It is critical to improve the performance of the algorithm by determining the mixing ratio between demonstration and self-generated data. Thus, we propose a prioritized storage mechanism to control this ratio automatically. In order to improve the robustness of the training process, a superior network is additionally introduced based on Double DQN, which always serves as a Q-network with competitive performance. We evaluate the performance of the proposed algorithm over three public research benchmarks and compare it against strong baselines, including three classical heuristics and one state-of-the-art imitation learning-based branching algorithm. The results show that the proposed algorithm achieves the best performance among compared algorithms and possesses the potential to improve B\&B algorithm performance continuously.
△ Less
Submitted 16 January, 2022;
originally announced January 2022.
-
A positive fraction Erdos-Szekeres theorem and its applications
Authors:
Andrew Suk,
Ji Zeng
Abstract:
A famous theorem of Erdos and Szekeres states that any sequence of $n$ distinct real numbers contains a monotone subsequence of length at least $\sqrt{n}$. Here, we prove a positive fraction version of this theorem. For $n > (k-1)^2$, any sequence $A$ of $n$ distinct real numbers contains a collection of subsets $A_1,\ldots, A_k \subset A$, appearing sequentially, all of size $s=Ω(n/k^2)$, such th…
▽ More
A famous theorem of Erdos and Szekeres states that any sequence of $n$ distinct real numbers contains a monotone subsequence of length at least $\sqrt{n}$. Here, we prove a positive fraction version of this theorem. For $n > (k-1)^2$, any sequence $A$ of $n$ distinct real numbers contains a collection of subsets $A_1,\ldots, A_k \subset A$, appearing sequentially, all of size $s=Ω(n/k^2)$, such that every subsequence $(a_1,\ldots, a_k)$, with $a_i \in A_i$, is increasing, or every such subsequence is decreasing. The subsequence $S = (A_1,\ldots, A_k)$ described above is called block-monotone of depth $k$ and block-size $s$. Our theorem is asymptotically best possible and follows from a more general Ramsey-type result for monotone paths, which we find of independent interest. We also show that for any positive integer $k$, any finite sequence of distinct real numbers can be partitioned into $O(k^2\log k)$ block-monotone subsequences of depth at least $k$, upon deleting at most $(k-1)^2$ entries. We apply our results to mutually avoiding planar point sets and biarc diagrams in graph drawing.
△ Less
Submitted 1 April, 2023; v1 submitted 3 December, 2021;
originally announced December 2021.
-
Mixed displacement-pressure-phase field framework for finite strain fracture of nearly incompressible hyperelastic materials
Authors:
Fucheng Tian,
Jun Zeng,
Mengnan Zhang,
Liangbin Li
Abstract:
The favored phase field method (PFM) has encountered challenges in the finite strain fracture modeling of nearly or truly incompressible hyperelastic materials. We identified that the underlying cause lies in the innate contradiction between incompressibility and smeared crack opening. Drawing on the stiffness-degradation idea in PFM, we resolved this contradiction through loosening incompressible…
▽ More
The favored phase field method (PFM) has encountered challenges in the finite strain fracture modeling of nearly or truly incompressible hyperelastic materials. We identified that the underlying cause lies in the innate contradiction between incompressibility and smeared crack opening. Drawing on the stiffness-degradation idea in PFM, we resolved this contradiction through loosening incompressible constraint of the damaged phase without affecting the incompressibility of intact material. By modifying the perturbed Lagrangian approach, we derived a novel mixed formulation. In numerical aspects, the finite element discretization uses the classical Q1/P0 and high-order P2/P1 schemes, respectively. To ease the mesh distortion at large strains, an adaptive mesh deletion technology is also developed. The validity and robustness of the proposed mixed framework are corroborated by four representative numerical examples. By comparing the performance of Q1/P0 and P2/P1, we conclude that the Q1/P0 formulation is a better choice for finite strain fracture in nearly incompressible cases. Moreover, the numerical examples also show that the combination of the proposed framework and methodology has vast potential in simulating complex peeling and tearing problems
△ Less
Submitted 1 December, 2021;
originally announced December 2021.
-
A propagation of chaos result for weakly interacting nonlinear Snell envelopes
Authors:
Boualem Djehiche,
Roxana Dumitrescu,
Jia Zeng
Abstract:
In this article, we establish a propagation of chaos result for weakly interacting nonlinear Snell envelopes which converge to a class of mean-field reflected backward stochastic differential equations (BSDEs) with jumps and right-continuous and left-limited obstacle, where the mean-field interaction in terms of the distribution of the $Y$-component of the solution enters both the driver and the l…
▽ More
In this article, we establish a propagation of chaos result for weakly interacting nonlinear Snell envelopes which converge to a class of mean-field reflected backward stochastic differential equations (BSDEs) with jumps and right-continuous and left-limited obstacle, where the mean-field interaction in terms of the distribution of the $Y$-component of the solution enters both the driver and the lower obstacle. Under mild Lipschitz and integrability conditions on the coefficients, we prove existence and uniqueness of the solution to both the mean-field reflected BSDEs with jumps and the corresponding system of weakly interacting particles and provide a propagation of chaos result for the whole solution $(Y,Z,U,K)$, which requires new technical results due to the dependence of the obstacle on the solution and the presence of jumps.
△ Less
Submitted 7 May, 2022; v1 submitted 28 November, 2021;
originally announced November 2021.
-
Safety-Critical Control and Planning for Obstacle Avoidance between Polytopes with Control Barrier Functions
Authors:
Akshay Thirugnanam,
Jun Zeng,
Koushil Sreenath
Abstract:
Obstacle avoidance between polytopes is a challenging topic for optimal control and optimization-based trajectory planning problems. Existing work either solves this problem through mixed-integer optimization, relying on simplification of system dynamics, or through model predictive control with dual variables using distance constraints, requiring long horizons for obstacle avoidance. In either ca…
▽ More
Obstacle avoidance between polytopes is a challenging topic for optimal control and optimization-based trajectory planning problems. Existing work either solves this problem through mixed-integer optimization, relying on simplification of system dynamics, or through model predictive control with dual variables using distance constraints, requiring long horizons for obstacle avoidance. In either case, the solution can only be applied as an offline planning algorithm. In this paper, we exploit the property that a smaller horizon is sufficient for obstacle avoidance by using discrete-time control barrier function (DCBF) constraints and we propose a novel optimization formulation with dual variables based on DCBFs to generate a collision-free dynamically-feasible trajectory. The proposed optimization formulation has lower computational complexity compared to existing work and can be used as a fast online algorithm for control and planning for general nonlinear dynamical systems. We validate our algorithm on different robot shapes using numerical simulations with a kinematic bicycle model, resulting in successful navigation through maze environments with polytopic obstacles.
△ Less
Submitted 30 May, 2022; v1 submitted 25 September, 2021;
originally announced September 2021.
-
On average hitting time and Kemeny's constant for weighted trees
Authors:
Ji Zeng
Abstract:
For a connected graph $G$, the average hitting time $α(G)$ and the Kemeny's constant $κ(G)$ are two similar quantities, both measuring the time for the random walk on $G$ to travel between two randomly chosen vertices. We prove that, among all weighted trees whose edge-weights form a fixed multiset, $α$ is maximized by a special type of "polarized" paths and is minimized by a unique weighted star…
▽ More
For a connected graph $G$, the average hitting time $α(G)$ and the Kemeny's constant $κ(G)$ are two similar quantities, both measuring the time for the random walk on $G$ to travel between two randomly chosen vertices. We prove that, among all weighted trees whose edge-weights form a fixed multiset, $α$ is maximized by a special type of "polarized" paths and is minimized by a unique weighted star graph. We also give a short proof of the fact that, among all simple trees of a fixed size, $κ$ is maximized by the path and is minimized by the star graph. Our proofs are based on the forest formulas for the average hitting time and the Kemeny's constant.
△ Less
Submitted 19 February, 2022; v1 submitted 19 September, 2021;
originally announced September 2021.
-
Forest formulas of discrete Green's functions
Authors:
Fan Chung,
Ji Zeng
Abstract:
The discrete Green's functions are the pseudoinverse (or the inverse) of the Laplacian (or its variations) of a graph. In this paper, we will give combinatorial interpretations of Green's functions in terms of enumerating trees and forests in a graph that will be used to derive further formulas for several graph invariants. For example, we show that the trace of the Green's function $\mathbf{G}$ a…
▽ More
The discrete Green's functions are the pseudoinverse (or the inverse) of the Laplacian (or its variations) of a graph. In this paper, we will give combinatorial interpretations of Green's functions in terms of enumerating trees and forests in a graph that will be used to derive further formulas for several graph invariants. For example, we show that the trace of the Green's function $\mathbf{G}$ associated with the combinatorial Laplacian of a connected simple graph $Γ$ on $n$ vertices satisfies $\text{Tr}(\mathbf{G})=\sum_{λ_i \neq 0} \frac 1 {λ_i}= \frac{1}{nτ}|\mathbb{F}^*_2|$, where $λ_i$ denotes the eigenvalues of the combinatorial Laplacian, $τ$ denotes the number of spanning trees and $\mathbb{F}^*_2$ denotes the set of rooted spanning $2$-forests in $Γ$. We will prove forest formulas for discrete Green's functions for directed and weighted graphs and apply them to study random walks on graphs and digraphs. We derive a forest expression of the hitting time for digraphs, which gives combinatorial proofs to old and new results about hitting times, traces of discrete Green's functions, and other related quantities.
△ Less
Submitted 17 August, 2022; v1 submitted 3 September, 2021;
originally announced September 2021.
-
Cycles of even-odd drop permutations and continued fractions of Genocchi numbers
Authors:
Qiongqiong Pan,
Jiang Zeng
Abstract:
Recently, Lazar and Wachs (arXiv:1910.07651) showed that the (median) Genocchi numbers play a fundamental role in the study of the homogenized Linial arrangement and obtained two new permutation models (called D-permutations and E-permutations) for (median) Genocchi numbers. They further conjecture that the distributions of cycle numbers over the two models are equal. In a follow-up, Eu et al. (ar…
▽ More
Recently, Lazar and Wachs (arXiv:1910.07651) showed that the (median) Genocchi numbers play a fundamental role in the study of the homogenized Linial arrangement and obtained two new permutation models (called D-permutations and E-permutations) for (median) Genocchi numbers. They further conjecture that the distributions of cycle numbers over the two models are equal. In a follow-up, Eu et al. (arXiv:2103.09130) further proved the gamma-positivity of the descent polynomials of even-odd descent permutations, which are in bijection with E-permutations by Foata's fundamental transformation. This paper merges the above two papers by considering a general moment sequence which encompasses the number of cycles and number of drops of E-permutations. Using the combinatorial theory of continued fraction, the moment connection enables us to confirm Lazar-Wachs' conjecture and obtain a natural $(p,q)$-analogue of Eu et al's descent polynomials. Furthermore, we show that the $γ$-coefficients of our $(p,q)$-analogue of descent polynomials have the same factorization flavor as the $γ$-coeffcients of Brändén's $(p,q)$-Eulerian polynomials.
△ Less
Submitted 10 August, 2021; v1 submitted 6 August, 2021;
originally announced August 2021.
-
Duality-based Convex Optimization for Real-time Obstacle Avoidance between Polytopes with Control Barrier Functions
Authors:
Akshay Thirugnanam,
Jun Zeng,
Koushil Sreenath
Abstract:
Developing controllers for obstacle avoidance between polytopes is a challenging and necessary problem for navigation in tight spaces. Traditional approaches can only formulate the obstacle avoidance problem as an offline optimization problem. To address these challenges, we propose a duality-based safety-critical optimal control using nonsmooth control barrier functions for obstacle avoidance bet…
▽ More
Developing controllers for obstacle avoidance between polytopes is a challenging and necessary problem for navigation in tight spaces. Traditional approaches can only formulate the obstacle avoidance problem as an offline optimization problem. To address these challenges, we propose a duality-based safety-critical optimal control using nonsmooth control barrier functions for obstacle avoidance between polytopes, which can be solved in real-time with a QP-based optimization problem. A dual optimization problem is introduced to represent the minimum distance between polytopes and the Lagrangian function for the dual form is applied to construct a control barrier function. We validate the obstacle avoidance with the proposed dual formulation for L-shaped (sofa-shaped) controlled robot in a corridor environment. We demonstrate real-time tight obstacle avoidance with non-conservative maneuvers on a moving sofa (piano) problem with nonlinear dynamics.
△ Less
Submitted 18 April, 2022; v1 submitted 18 July, 2021;
originally announced July 2021.
-
Moments of Orthogonal Polynomials and Exponential Generating Functions
Authors:
Ira M. Gessel,
Jiang Zeng
Abstract:
Starting from the moment sequences of classical orthogonal polynomials we derive the orthogonality purely algebraically. We consider also the moments of ($q=1$) classical orthogonal polynomials, and study those cases in which the exponential generating function has a nice form. In the opposite direction, we show that the generalized Dumont-Foata polynomials with six parameters are the moments of r…
▽ More
Starting from the moment sequences of classical orthogonal polynomials we derive the orthogonality purely algebraically. We consider also the moments of ($q=1$) classical orthogonal polynomials, and study those cases in which the exponential generating function has a nice form. In the opposite direction, we show that the generalized Dumont-Foata polynomials with six parameters are the moments of rescaled continuous dual Hahn polynomials. Finally we show that one of our methods can be applied to deal with the moments of Askey-Wilson polynomials.
△ Less
Submitted 9 January, 2022; v1 submitted 1 July, 2021;
originally announced July 2021.
-
Tilting objects in singularity categories of toric Gorenstein varieties
Authors:
Xiaojun Chen,
Leilei Liu,
Jieheng Zeng
Abstract:
We study certain toric Gorenstein varieties with isolated singularities which are the quotient spaces of generic unimodular representations by the one-dimensional torus, or by the product of the one-dimensional torus with a finite abelian group. Based on the works of Špenko and Van den Bergh [Invent. Math. 210 (2017), no. 1, 3-67] and Mori and Ueyama [Adv. Math. 297 (2016), 54-92], we show that th…
▽ More
We study certain toric Gorenstein varieties with isolated singularities which are the quotient spaces of generic unimodular representations by the one-dimensional torus, or by the product of the one-dimensional torus with a finite abelian group. Based on the works of Špenko and Van den Bergh [Invent. Math. 210 (2017), no. 1, 3-67] and Mori and Ueyama [Adv. Math. 297 (2016), 54-92], we show that the singularity categories of these varieties admit tilting objects, and hence are triangle equivalent to the perfect categories of some finite dimensional algebras.
△ Less
Submitted 23 February, 2022; v1 submitted 27 June, 2021;
originally announced June 2021.
-
Enhancing Feasibility and Safety of Nonlinear Model Predictive Control with Discrete-Time Control Barrier Functions
Authors:
Jun Zeng,
Zhongyu Li,
Koushil Sreenath
Abstract:
Safety is one of the fundamental problems in robotics. Recently, one-step or multi-step optimal control problems for discrete-time nonlinear dynamical system were formulated to offer tracking stability using control Lyapunov functions (CLFs) while subject to input constraints as well as safety-critical constraints using control barrier functions (CBFs). The limitations of these existing approaches…
▽ More
Safety is one of the fundamental problems in robotics. Recently, one-step or multi-step optimal control problems for discrete-time nonlinear dynamical system were formulated to offer tracking stability using control Lyapunov functions (CLFs) while subject to input constraints as well as safety-critical constraints using control barrier functions (CBFs). The limitations of these existing approaches are mainly about feasibility and safety. In the existing approaches, the feasibility of the optimization and the system safety cannot be enhanced at the same time theoretically. In this paper, we propose two formulations that unifies CLFs and CBFs under the framework of nonlinear model predictive control (NMPC). In the proposed formulations, safety criteria is commonly formulated as CBF constraints and stability performance is ensured with either a terminal cost function or CLF constraints. Slack variables with relaxing technique are introduced on the CBF constraints to resolve the tradeoff between feasibility and safety so that they can be enhanced at the same. The advantages about feasibility and safety of proposed formulations compared with existing methods are analyzed theoretically and validated with numerical results.
△ Less
Submitted 1 October, 2021; v1 submitted 21 May, 2021;
originally announced May 2021.
-
Batalin-Vilkovisky algebra structure on Poisson manifolds with diagonalizable modular symmetry
Authors:
Xiaojun Chen,
Leilei Liu,
Sirui Yu,
Jieheng Zeng
Abstract:
We study the ``twisted" Poincaré duality of smooth Poisson manifolds, and show that, if the modular vector field is diagonalizable, then there is a mixed complex associated to the Poisson complex, which, combining with the twisted Poincaré duality, gives a Batalin-Vilkovisky algebra structure on the Poisson cohomology. This generalizes the previous results obtained by Xu for unimodular Poisson man…
▽ More
We study the ``twisted" Poincaré duality of smooth Poisson manifolds, and show that, if the modular vector field is diagonalizable, then there is a mixed complex associated to the Poisson complex, which, combining with the twisted Poincaré duality, gives a Batalin-Vilkovisky algebra structure on the Poisson cohomology. This generalizes the previous results obtained by Xu for unimodular Poisson manifolds. We also show that the Batalin-Vilkovisky algebra structure is preserved under Kontsevich's deformation quantization, and in the case of polynomial algebras it is also preserved by Koszul duality.
△ Less
Submitted 2 April, 2023; v1 submitted 29 April, 2021;
originally announced April 2021.
-
Variational approach to the existence of solutions for non-instantaneous impulsive differential equations with perturbation
Authors:
Wangjin Yao,
Liping Dong,
Jing Zeng
Abstract:
In this paper, we study the existence of solutions for second-order non-instantaneous impulsive differential equations with a perturbation term. By variational approach, we obtain the problem has at least one solution under assumptions that the nonlinearities are super-quadratic at infinity, and sub-quadratic at the origin.
In this paper, we study the existence of solutions for second-order non-instantaneous impulsive differential equations with a perturbation term. By variational approach, we obtain the problem has at least one solution under assumptions that the nonlinearities are super-quadratic at infinity, and sub-quadratic at the origin.
△ Less
Submitted 30 March, 2021;
originally announced March 2021.
-
Equidistributions around special kinds of descents and excedances
Authors:
Bin Han,
Jianxi Mao,
Jiang Zeng
Abstract:
We consider a sequence of four variable polynomials by refining Stieltjes' continued fraction for Eulerian polynomials. Using combinatorial theory of Jacobi-type continued fractions and bijections we derive various combinatorial interpretations in terms of permutation statistics for these polynomials, which include special kinds of descents and excedances in a recent paper of Baril and Kirgizov. A…
▽ More
We consider a sequence of four variable polynomials by refining Stieltjes' continued fraction for Eulerian polynomials. Using combinatorial theory of Jacobi-type continued fractions and bijections we derive various combinatorial interpretations in terms of permutation statistics for these polynomials, which include special kinds of descents and excedances in a recent paper of Baril and Kirgizov. As a by-product, we derive several equidistribution results for permutation statistics, which enables us to confirm and strengthen a recent conjecture of Vajnovszki and also to obtain several compagnion permutation statistics for two bistatistics in a conjecture of Baril and Kirgizov.
△ Less
Submitted 8 September, 2021; v1 submitted 24 March, 2021;
originally announced March 2021.
-
Generalized Opinion Dynamics Model for Social Trust Networks in Opinion Maximization
Authors:
Changxiang He,
Jiayuan Zeng,
Shuting Liu,
Guang Zhang,
Xiaofei Qin,
Xuedian Zhang,
Lele Liu
Abstract:
In this paper, we propose a generalized opinion dynamics model (GODM), which can dynamically compute each person's expressed opinion, to solve the internal opinion maximization problem for social trust networks. In the model, we propose a new, reasonable and interpretable confidence index, which is determined by both person's social status and the evaluation around him. By using the theory of diag…
▽ More
In this paper, we propose a generalized opinion dynamics model (GODM), which can dynamically compute each person's expressed opinion, to solve the internal opinion maximization problem for social trust networks. In the model, we propose a new, reasonable and interpretable confidence index, which is determined by both person's social status and the evaluation around him. By using the theory of diagonally dominant, we obtain the optimal analytic solution of the Nash equilibrium with maximum overall opinion. We design a novel algorithm to maximize the overall with given budget by modifying the internal opinions of people in the social trust network, and prove its optimality both from the algorithm itself and the traditional optimization algorithm-ADMM algorithms with $l_1$-regulations. A series of experiments are conducted, and the experimental results show that our method is superior to the state-of-the-art in four datasets. The average benefit has promoted $67.5\%$, $83.2\%$, $31.5\%$, and $33.7\%$ on four datasets, respectively.
△ Less
Submitted 9 March, 2021;
originally announced March 2021.
-
Moreau Envelope Augmented Lagrangian Method for Nonconvex Optimization with Linear Constraints
Authors:
Jinshan Zeng,
Wotao Yin,
Ding-Xuan Zhou
Abstract:
The augmented Lagrangian method (ALM) is one of the most useful methods for constrained optimization. Its convergence has been well established under convexity assumptions or smoothness assumptions, or under both assumptions. ALM may experience oscillations and divergence when the underlying problem is simultaneously nonconvex and nonsmooth. In this paper, we consider the linearly constrained prob…
▽ More
The augmented Lagrangian method (ALM) is one of the most useful methods for constrained optimization. Its convergence has been well established under convexity assumptions or smoothness assumptions, or under both assumptions. ALM may experience oscillations and divergence when the underlying problem is simultaneously nonconvex and nonsmooth. In this paper, we consider the linearly constrained problem with a nonconvex (in particular, weakly convex) and nonsmooth objective. We modify ALM to use a Moreau envelope of the augmented Lagrangian and establish its convergence under conditions that are weaker than those in the literature. We call it the Moreau envelope augmented Lagrangian (MEAL) method. We also show that the iteration complexity of MEAL is $o(\varepsilon^{-2})$ to yield an $\varepsilon$-accurate first-order stationary point. We establish its whole sequence convergence (regardless of the initial guess) and a rate when a Kurdyka-Lojasiewicz property is assumed. Moreover, when the subproblem of MEAL has no closed-form solution and is difficult to solve, we propose two practical variants of MEAL, an inexact version called iMEAL with an approximate proximal update, and a linearized version called LiMEAL for the constrained problem with a composite objective. Their convergence is also established.
△ Less
Submitted 8 December, 2021; v1 submitted 21 January, 2021;
originally announced January 2021.
-
MFG model with a long-lived penalty at random jump times: application to demand side management for electricity contracts
Authors:
Clémence Alasseur,
Luciano Campi,
Roxana Dumitrescu,
Jia Zeng
Abstract:
We consider an energy system with $n$ consumers who are linked by a Demand Side Management (DSM) contract, i.e. they agreed to diminish, at random times, their aggregated power consumption by a predefined volume during a predefined duration. Their failure to deliver the service is penalised via the difference between the sum of the $n$ power consumptions and the contracted target. We are led to an…
▽ More
We consider an energy system with $n$ consumers who are linked by a Demand Side Management (DSM) contract, i.e. they agreed to diminish, at random times, their aggregated power consumption by a predefined volume during a predefined duration. Their failure to deliver the service is penalised via the difference between the sum of the $n$ power consumptions and the contracted target. We are led to analyse a non-zero sum stochastic game with $n$ players, where the interaction takes place through a cost which involves a delay induced by the duration included in the DSM contract. When $n \to \infty$, we obtain a Mean-Field Game (MFG) with random jump time penalty and interaction on the control. We prove a stochastic maximum principle in this context, which allows to compare the MFG solution to the optimal strategy of a central planner. In a linear quadratic setting we obtain an semi-explicit solution through a system of decoupled forward-backward stochastic differential equations with jumps, involving a Riccati Backward SDE with jumps. We show that it provides an approximate Nash equilibrium for the original $n$-player game for $n$ large. Finally, we propose a numerical algorithm to compute the MFG equilibrium and present several numerical experiments.
△ Less
Submitted 15 January, 2021;
originally announced January 2021.
-
On Stochastic Variance Reduced Gradient Method for Semidefinite Optimization
Authors:
Jinshan Zeng,
Yixuan Zha,
Ke Ma,
Yuan Yao
Abstract:
The low-rank stochastic semidefinite optimization has attracted rising attention due to its wide range of applications. The nonconvex reformulation based on the low-rank factorization, significantly improves the computational efficiency but brings some new challenge to the analysis. The stochastic variance reduced gradient (SVRG) method has been regarded as one of the most effective methods. SVRG…
▽ More
The low-rank stochastic semidefinite optimization has attracted rising attention due to its wide range of applications. The nonconvex reformulation based on the low-rank factorization, significantly improves the computational efficiency but brings some new challenge to the analysis. The stochastic variance reduced gradient (SVRG) method has been regarded as one of the most effective methods. SVRG in general consists of two loops, where a reference full gradient is first evaluated in the outer loop and then used to yield a variance reduced estimate of the current gradient in the inner loop. Two options have been suggested to yield the output of the inner loop, where Option I sets the output as its last iterate, and Option II yields the output via random sampling from all the iterates in the inner loop. However, there is a significant gap between the theory and practice of SVRG when adapted to the stochastic semidefinite programming (SDP). SVRG practically works better with Option I, while most of existing theoretical results focus on Option II. In this paper, we fill this gap via exploiting a new semi-stochastic variant of the original SVRG with Option I adapted to the semidefinite optimization. Equipped with this, we establish the global linear submanifold convergence (i.e., converging exponentially fast to a submanifold of a global minimum under the orthogonal group action) of the proposed SVRG method, given a provable initialization scheme and under certain smoothness and restricted strongly convex assumptions. Our analysis includes the effects of the mini-batch size and update frequency in the inner loop as well as two practical step size strategies, the fixed and stabilized Barzilai-Borwein step sizes. Some numerical results in matrix sensing demonstrate the efficiency of proposed SVRG method outperforming Option II counterpart as well as others.
△ Less
Submitted 1 January, 2021;
originally announced January 2021.
-
Weighted and maximally hypoelliptic estimates for the Fokker-Planck Operator with electromagnetic fields
Authors:
Wei-Xi Li,
Juan Zeng
Abstract:
We consider a Fokker-Planck operator with electric potential and electromagnetic fields. We establish the sharp weighted and subelliptic estimates, involving the control of the derivatives of electric potential and electromagnetic fields. Our proof relies on a localization argument as well as a careful calculation on commutators.
We consider a Fokker-Planck operator with electric potential and electromagnetic fields. We establish the sharp weighted and subelliptic estimates, involving the control of the derivatives of electric potential and electromagnetic fields. Our proof relies on a localization argument as well as a careful calculation on commutators.
△ Less
Submitted 18 December, 2020;
originally announced December 2020.
-
Moments of q-Jacobi Polynomials and q-Zeta Values
Authors:
Frédéric Chapoton,
Christian Krattenthaler,
Jiang Zeng
Abstract:
We explore some connections between moments of rescaled little q-Jacobi polynomials, q-analogues of values at negative integers for some Dirichlet series, and the q-Eulerian polynomials of wreath products of symmetric groups.
We explore some connections between moments of rescaled little q-Jacobi polynomials, q-analogues of values at negative integers for some Dirichlet series, and the q-Eulerian polynomials of wreath products of symmetric groups.
△ Less
Submitted 3 December, 2020;
originally announced December 2020.
-
Block Form of Frobenius Groups
Authors:
Jiwen Zeng,
Jiping Zhang
Abstract:
The aim of this paper is to apply character properties of Frobenius group to a local block form of an group algebra. We start by establishing a block form of Brauer permutation Lemma by using block participation of conjugate classes of a group $G$. Then we can define a pair of Frobenius corresponding blocks between a group $G$ and its normal subgroup $N$. A near group condition is given to determi…
▽ More
The aim of this paper is to apply character properties of Frobenius group to a local block form of an group algebra. We start by establishing a block form of Brauer permutation Lemma by using block participation of conjugate classes of a group $G$. Then we can define a pair of Frobenius corresponding blocks between a group $G$ and its normal subgroup $N$. A near group condition is given to determine a pair of Frobenius corresponding blocks. With a pair of Frobenius corresponding blocks, we study its group structure. At last we prove connections between nilpotent properties and Frobenius corresponding blocks.
△ Less
Submitted 28 October, 2020; v1 submitted 28 October, 2020;
originally announced October 2020.
-
Minor summation formula of hyperpfaffians and Selberg integrals
Authors:
Masao Ishikawa,
Jiang Zeng
Abstract:
In the previous paper (J. Combin. Theory Ser. A, 120, 2013, 1263--1284) H. Tagawa and the two authors proposed an algebraic method to compute certain Pfaffians whose form resemble to Hankel determinants associated with moment sequences of the classical orthogonal polynomials. At the end of the paper they offered several conjectures. In this work we employ a completely different method to evaluate…
▽ More
In the previous paper (J. Combin. Theory Ser. A, 120, 2013, 1263--1284) H. Tagawa and the two authors proposed an algebraic method to compute certain Pfaffians whose form resemble to Hankel determinants associated with moment sequences of the classical orthogonal polynomials. At the end of the paper they offered several conjectures. In this work we employ a completely different method to evaluate this type of Pfaffians. The idea is to apply certain de Bruijn type formula and convert the evaluation of the Pfaffians to the certain Selberg type integrals. This approach works not only for Pfaffians but also for hyperpfaffians. Hence it enables us to establish much more generalized identities than those conjectured in the previous paper. We also investigate some Pfaffians related to classical $q$-orthogonal polynomials.
△ Less
Submitted 20 October, 2022; v1 submitted 22 August, 2020;
originally announced August 2020.
-
Equidistributions of mesh patterns of length two and Kitaev and Zhang's conjectures
Authors:
Bin Han,
Jiang Zeng
Abstract:
A systematic study of avoidance of mesh patterns of length 2 was conducted by Hilmarsson et al. in 2015. In a recent paper Kitaev and Zhang examined the distribution of the aforementioned patterns. The aim of this paper is to prove more equidistributions of mesh pattern and confirm Kitaev and Zhang's four conjectures by constructing two involutions on permutations.
A systematic study of avoidance of mesh patterns of length 2 was conducted by Hilmarsson et al. in 2015. In a recent paper Kitaev and Zhang examined the distribution of the aforementioned patterns. The aim of this paper is to prove more equidistributions of mesh pattern and confirm Kitaev and Zhang's four conjectures by constructing two involutions on permutations.
△ Less
Submitted 25 July, 2020;
originally announced July 2020.