I am a Ph.D. student at USTC, currently under the guidance of Shuai Shao.
Prior to this, I earned my Bachelor's degree in Computer Science from UESTC while conducting theoretical computer science research under the guidance of Chao Xu.
During the first three years of my Ph.D., I focused on theoretical computer science, including algorithms, combinatorial optimization, and approximate counting. I have since shifted my focus and am now exploring more practical, application-oriented directions. Here is my CV.
You can contact me through email self.ke.shi@gmail.com , ke.shi@ustc.edu or chihaya.anon@ustc.edu.
We investigate zero-free regions of the partition function of the hard-core model on finite graphs and their implications for the analyticity of free energy on infinite lattices. Classical results establish zero-freeness up to the tree uniqueness threshold based on maximum degree, but the connective constant offers a tighter structural measure for graphs like regular lattices. We bridge this gap by defining the connective constant for finite graphs via a lower bound on the number of self-avoiding walks. We prove zero-freeness in a complex neighborhood for graph families with lower connective constant \mu when \lambda < \lambda_c(\mu), establishing uniqueness and analyticity of the free energy density for infinite lattices beyond maximum-degree bounds. The proof uses a block contraction technique that lifts the correlation decay property from a real interval to a strip-like complex neighborhood.
We study spatial correlation decay, focusing on strong spatial mixing (SSM), for continuous spin systems with multi-interactions on hypergraphs. For several classes of discrete spin systems, zero-freeness of the partition function and SSM have been related in both directions. For continuous spin systems, this relationship is less understood. In this paper, we prove that a uniform pinned zero-free disk about \beta = 0 implies SSM throughout that disk for continuous spin systems with bounded interactions. Combined with the zero-freeness result of Barvinok, this gives SSM for systems whose interactions are 1-Lipschitz in the Hamming metric, on hypergraphs with maximum degree \Delta and maximum edge size r, for all complex \beta with |\beta| < 1/(3\Delta\sqrt{r-1}). The SSM constants are uniform over the system size, the product reference measures, and all valid product boundary events.
A key technical contribution of this paper is a divisibility relation for partition functions of continuous spin systems. We follow the local dependence of coefficients (LDC) framework developed for discrete systems and extend its divisibility step to integral partition functions. The Taylor coefficients in our divisibility identity are sums of integrals of products of interactions, with repeated hyperedges allowed. A measure-preserving swap of two configurations, together with a reindexing of the interaction factors, cancels all contributions below the distance between the boundary supports. This shows how the LDC and zero-freeness approaches fit together; algorithmic running times additionally depend on how the continuous interactions and the required integrals are represented and computed.
We present a Weitz-type FPTAS for the ferromagnetic Ising model across the entire Lee–Yang zero-free region, without relying on the strong spatial mixing (SSM) property. Our algorithm is Weitz-type for two reasons. First, it expresses the partition function as a telescoping product of ratios, with the key being to approximate each ratio. Second, it uses Weitz’s self-avoiding walk tree, and truncates it at logarithmic depth to give a good and efficient approximation. The key difference from the standard Weitz algorithm is that we approximate a carefully designed edge-deletion ratio instead of the marginal probability of a vertex’s spin, ensuring our algorithm does not require SSM.
Furthermore, by establishing local dependence of coefficients (LDC), we indeed prove a novel form of SSM for these edge-deletion ratios, which, in turn, implies the standard SSM for the random cluster model. This is the first SSM result for the random cluster model on general graphs, beyond lattices. We prove LDC using a new division relation, and remarkably, such relations hold quite universally. As a result, we establish LDC for a variety of models. Combined with existing zero-freeness results for these models, we derive new SSM results for them. Our work suggests that both Weitz-type FPTASes and SSM can be derived from zero-freeness, while zero-freeness alone suffices for Weitz-type FPTASes, SSM additionally requires LDC, a combinatorial property independent of zero-freeness.
The Stacker Crane Problem (SCP) is a variant of the Traveling Salesman Problem. In SCP, pairs of pickup and delivery points are designated on a graph, and a crane must visit these points to move objects from each pickup location to its respective delivery point. The goal is to minimize the total distance traveled. SCP is known to be NP-hard, even on tree structures. The only positive results, in terms of polynomial-time solvability, apply to graphs that are topologically equivalent to a path or a cycle. We propose an algorithm that is optimal for each fixed topology, running in near-linear time. This is achieved by demonstrating that the problem is fixed-parameter tractable (FPT) when parameterized by both the cycle rank and the number of branch vertices.
COCOON .
A subset B of ring \Z_n is called a \ell-covering set if \{ ab \pmod n \mid 0\leq a \leq \ell, b\in B\} = \Z_n. We show there exists a \ell-covering set of \Z_n of size O(\frac{n}{\ell}\log n) for all n and \ell, and how to construct such set. We also show examples where any \ell-covering set must have size \Omega(\frac{n}{\ell}\frac{\log n}{\log \log n}). The proof uses a refined bound for relative totient function obtained through sieve theory, and existence of a large divisor with linear divisor sum.
In this paper, we study the minimum k-partition problem of submodular functions, i.e., given a finite set V and a submodular function f:2^V\to \R, computing a k-partition \{ V_1, \ldots, V_k \} of V with minimum \sum_{i=1}^k f(V_i). The problem is a natural generalization of the minimum k-cut problem in graphs and hypergraphs. It is known that the problem is NP-hard for general k, and solvable in polynomial time for k \leq 3. In this paper, we construct the first polynomial-time algorithm for the minimum 4-partition problem.
Mathematical Programming .
In this paper, we study the minimum k-partition problem of submodular functions, i.e., given a finite set V and a submodular function f:2^V\to \R, computing a k-partition \{ V_1, \ldots, V_k \} of V with minimum \sum_{i=1}^k f(V_i). The problem is a natural generalization of the minimum k-cut problem in graphs and hypergraphs. It is known that the problem is NP-hard for general k, and solvable in polynomial time for k \leq 3. In this paper, we construct the first polynomial-time algorithm for the minimum 4-partition problem.
2026, Manuscript.
We establish complex zero-freeness and deterministic fully polynomial-time approximation schemes (FPTASes) for spin systems on the Erdős–Rényi random graph G(n,d/n). For the antiferromagnetic q-state Potts model (q\ge3) with d<d_2(q)\sim2q\log q, the partition function is with high probability zero-free on a fixed-width neighborhood of [0,1], including proper q-colorings at z=0; this yields a deterministic FPTAS for q\ge(1/2+\eta)d/\log d (fixed \eta>0, large d), and, combined with the Achlioptas–Naor noncolorability bound, a deterministic FPTAS for counting proper q-colorings for all but at most one integer q at each fixed d>0. An appendix extends the zero-freeness and the FPTAS up to the Erdős–Rényi condensation threshold of the Potts model, including proper colorings. For the zero-field Ising antiferromagnet and the zero-field Viana–Bray diluted spin glass on G(n,d/n), we likewise obtain fixed-width zero-free neighborhoods and deterministic FPTASes up to their Kesten–Stigum thresholds. Dividing out each graph's forest weight removes the order-\sqrt n fluctuations of the log-partition function caused by the random number of edges; complex moment bounds and Poisson short-cycle Taylor truncations then transfer nonvanishing from an explicit random cycle model, compensated where necessary, while connected-set Möbius inversion and Cayley's tree bound handle unbounded degrees in Barvinok interpolation.
2026, Submitted.
For \Delta \geq 2 and q \geq 11\Delta/6, we prove that the antiferromagnetic q-state Potts partition function on finite simple graphs of maximum degree at most \Delta has no Fisher zeros in a graph-uniform complex neighbourhood of [0, 1]. For q > 11\Delta/6, the proof establishes coupling independence throughout [0, 1] using a soft version of Vigoda's flip dynamics. Our main tool is a separator-shell transfer theorem for the Potts model on induced-subgraph closed classes of graphs of maximum degree at most \Delta, with q \geq \Delta + 1. It yields a graph-uniform zero-free neighbourhood of [0, 1] from Hamming coupling independence at 0 and a uniform coupling-independence bound on each interval [\delta, 1], \delta \in (0, 1]. We also obtain zero-free Lee--Yang polydiscs around the uniform field for vertex- and edge-colour fields. For Boolean Holant problems on bounded-degree graphs whose signatures come from a fixed finite family of log-concave symmetric signatures f with f(0) > 0, such as b-matchings, we obtain graph-uniform zero-free polytubes around every bounded box of nonnegative activities; their union is an open zero-free neighbourhood of the nonnegative orthant. An appendix summarizes further coupling-independence inputs and the zero-free regions they yield.
2026, Submitted.
We establish a complete complexity dichotomy for deciding whether a graph admits an Eulerian orientation satisfying local Boolean constraints at vertices. We denote this problem by \mathrm{EO}. For every fixed finite set of constraints, the problem is either polynomial-time solvable or NP-complete. The classification identifies two new tractable classes. In particular, one such class is obtained by exploiting the promise polymorphism condition used for proving tractability in counting \mathrm{EO} [GSS26]. For the hardness proof, we develop a nonconstructive reduction based on unsatisfiable instances and show how failures of the promise conditions yield obstructions to affine and bijunctive tractability, leading to NP-hardness through Schaefer's dichotomy.
Combining our \mathrm{EO} dichotomy with the classification of bipartite edge constraint satisfaction problems, we also classify the complexity of general graph orientation problems assuming both unary pinning constraints are available. Every such problem is polynomial-time solvable or NP-complete, except for a residual family of unresolved \Delta-matroid parity cases. These results advance the classification of constraint satisfaction problems with restricted variable occurrences and connect their tractability to promise polymorphisms.
2026, Manuscript.
The weighted Eulerian orientation counting problem (\#\mathrm{EO}) plays a key role in the complexity classification program for Holant problems. A recent result established an \mathrm{FP}^{\mathrm{NP}} versus \#\mathrm P-hard dichotomy for \#\mathrm{EO} problems. On the tractable side of this dichotomy, the supports of all functions admit the ternary XOR operation as a promise polymorphism relaxed in the same direction. It was left open whether all cases on the \mathrm{FP}^{\mathrm{NP}} side are in fact in FP. In this paper, we settle this question by giving a polynomial-time algorithm for all cases on the \mathrm{FP}^{\mathrm{NP}} side. Consequently, we obtain a complete FP versus \#\mathrm P-hard dichotomy for counting weighted Eulerian orientations, and further for complex-valued Holant problems with a nonzero signature of odd arity.
Our algorithm is based on a linear programming relaxation, but we use it in a nonstandard way. Instead of proving that the relaxation is integral and solving the problem directly from an optimal LP solution, we use the relaxation as a structural tool to lift the promise polymorphism condition to an ordinary polymorphism condition. This reveals an affine local structure of the constraint functions, which leads to tractability. A key technical innovation of our paper is to exploit the one-sided promise polymorphism through cubic forms of LP constraints, which interact precisely with the desired ternary XOR polymorphism.
2025, Submitted.
We study the connections between strong spatial mixing and zero-freeness, the two main notions used to devise deterministic approximation algorithms for counting problems. Following a recent framework introduced by [Regts, G. (2023). Absence of zeros implies strong spatial mixing. Probability Theory and Related Fields, 186(1), 621-641.] for the vertex model, we show that the implication from zero-freeness of the partition function to strong spatial mixing also works for the matching polynomial, a typical counting problem in the edge model. Based on Heilmann and Lieb's results, a zero-free region C_\Delta=\C \backslash (-\infty,-\frac{1}{4(\Delta-1)}] of the partition function and a Christoffel--Darboux type identity, we give a very simple proof for strong spatial mixing of the matching polynomial on graphs of degree bounded by \Delta with any complex edge parameter z \in C_\Delta. This proof does not rely on the tree recursion which is commonly used for proving strong spatial mixing in previous results.