Theory of Computing Report

Herman Chernoff passed away on July 6, 5 days after turning 103. Ravi Boppana wrote a guest post about Chernoff's life for his 100th birthday. 

Let me talk about his most famous work, the Chernoff Bounds themselves.

If you have a coin that will be heads with probability \(p\), and you flip it \(n\) times, the expected number of heads is \(pn\). Informally Chernoff bounds says that for large \(n\) the number of heads will be quite close to \(pn\) with an exponentially small probability of being far away from \(pn\). For example, if you flip a coin with probability 30% chance of being heads 10,000 times, the probability that you will get at most 2500 heads is less than \(10^{-18}\).

More formally, for \(\delta \in (0,1)\), Chernoff's bound shows that

\[\Pr[|X - \mu| \geq \delta\mu] \leq 2e^{-\mu\delta^2/3}\]

for \(\mu = \mathbb{E}[X]\). I find the variation known as Hoeffding's inequality \[\Pr[|X - \mu| \geq t] \leq 2e^{-2t^2/n}\] easier to use for computational complexity.

Chernoff bounds play a major role in computational complexity, for example you can use Chernoff Bounds for probabilistic, quantum and interactive proof algorithms to reduce the error to exponentially small. That means using the union bound, a single random sequence will give the correct answer for all inputs, showing that BPP is in P/poly. 

Chernoff bounds are a key ingredient in the proof that MA (prover then verifier) is contained in AM (verifier then prover). Roughly you make the MA error so small that the same random coins work no matter what the prover might say. MA in AM plays a key role in the proof that NP in BPP implies PH in BPP which itself is necessary for Toda's theorem.

The proof of Chernoff bounds comes from a simple trick: Rather than bounding \(\Pr[X \geq a]\) directly via Markov's inequality on \(X\), you apply Markov to \(e^{tX}\) for a free parameter \(t > 0.\) Since \(e^{tX}\) is nonnegative,

\[\Pr[X \geq a] = \Pr[e^{tX} \geq e^{ta}] \leq \frac{\mathbb{E}[e^{tX}]}{e^{ta}},\]

and then you optimize over \(t\). The quantity \(\mathbb{E}[e^{tX}]\) is the moment generating function, and for a sum \(X = \sum_i X_i\) of independent variables it factorizes into \(\prod_i \mathbb{E}[e^{tX_i}]\). That factorization is what converts a linear tail bound into an exponential one, each independent term contributes multiplicatively, so the deviation probability shrinks like a product rather than a sum.

Chernoff did not think of himself as a theorist. 

People regard me as a theoretical statistician, but I’ve decided in recent years that I’m really an applied statistician. My theoretical insights have relied upon my work in thinking about applied problems.

A good lesson for us all.

By Lance Fortnow

Authors: Vaibhav Krishan, Sundar Vishwanathan

The class $ACC^0$ consists of Boolean functions that can be computed by constant-depth circuits of polynomial size with $AND, NOT$ and $MOD_m$ gates, where $m$ is a natural number. At the frontier of our understanding lies a widely believed conjecture asserting that $MAJORITY$ does not belong to $ACC^0$. A few years ago, Bhrushundi, Hosseini, Lovett and Rao (ITCS 2019) introduced torus polynomial approximations as an approach towards this conjecture. Torus polynomials approximate Boolean functions when the fractional part of their value on Boolean points is close to half the value of the function. They reduced the conjecture that $MAJORITY \notin ACC^0$ to a conjecture concerning the non-existence of low degree torus polynomials that approximate $MAJORITY$. We reduce the non-existence problem further, to a statement about finding feasible solutions for an infinite family of linear programs. The main advantage of this statement is that it allows for incremental progress, which means finding feasible solutions for successively larger collections of these programs. As an immediate first step, we find feasible solutions for a large class of these linear programs, leaving only a finite set for further consideration. Our method is inspired by the method of dual polynomials, which is used to study the approximate degree of Boolean functions. Using our method, we also propose a way to progress further. We prove several additional key results with the same method, including lower bounds for approximating the $AND$ function, lower bounds when the approximating polynomial is symmetric, showcasing the power of our machinery.

Authors: Bruno Cavalar, Susanna F. de Rezende, Matthew Gray, Rahul Santhanam

We show the following hardness results for monotone learning and approximation of monotone circuit size: 1. Under the Randomised Exponential-Time Hypothesis (rETH), it requires time $n^{Ω(\log n)}$ to PAC-learn monotone formulas with $n$ input bits and size $s(n) = n$ by monotone circuits of size $n^{(\log n)^{1-ε}}$, for every $ε> 0$. 2. Under the Randomised Exponential-Time Hypothesis (rETH), for any $δ> 0$, there is a polynomially bounded function $m$ such that $m^{1-δ}$-multiplicatively approximating the minimum monotone circuit size of a monotone function consistent with a sequence of $m(n)$ labelled examples $\{(x_i, b_i)\}$ over $n$-bit inputs requires time $m^{Ω(\log(m))}$. Our results are shown by a novel application of lifting arguments in proof and communication complexity to hardness of monotone learning, by building on the seminal result of Atserias and Müller (J. ACM, 2020) on hardness of automating Resolution proofs.

Authors: Ho-Lin Chen, Xiang Huang

Current analog complexity theory, built on the General-Purpose Analog Computer (GPAC) model and polynomial ODEs, allows unbounded state variables -- an assumption that is physically unrealistic for chemical reaction networks and other laboratory-scale analog computers. We develop a bounded analog complexity theory in which all state variables remain in compact intervals and physical time (wall-clock time) is the only diverging resource. Our main technical contribution is bounded surrogate compilation, a compilation framework that transforms unbounded polynomial ODE systems into bounded ones while preserving computational limits and time-to-precision guarantees. We prove that if a system is compiled into a bounded system through our algorithm, the wall-clock time of the compiled system is polynomial in the arc length and physical time of the original system. We exhibit concrete constructions demonstrating fine-grained bounded time complexity -- a tunable polynomial-degree family, a Lambert-$W$-based system achieving $Θ(r\log r)$ time-to-precision (where $r$ is the desired precision parameter, in nats: $|x(t)-α|

Authors: Mikhail V. Volkov, Yinfeng Zhu

We determine which digraphs admit an edge labeling by letters from a finite alphabet such that the resulting labeled digraph is a completely reachable automaton. Such digraphs are recognizable in polynomial time; however, the problem becomes NP-complete when the size of the label alphabet is fixed. We also classify the digraphs for which every edge labeling results in a completely reachable automaton.

Authors: Avner Kantor

Snooker folklore holds that no single stroke can pocket all twenty-one object balls. We examine the claim in an idealized but fully specified model of billiard dynamics. Within the model we exhibit an admissible configuration of the twenty-two balls and a stroke of the cue ball that pockets all twenty-one object balls, and we show that the set of such strokes has positive Lebesgue measure in the natural shot space: total clearances are not flukes of measure zero but open events. For the regulation opening configuration we conjecture the same and explain both why a simulation cannot settle the conjecture by brute force and what kind of computation could settle it in principle. Monte Carlo experiments in the same model estimate the probability P(k) that a uniformly random stroke pockets exactly k balls; the observed decay of P(k), extrapolated conditionally on the conjecture, places the probability of a total clearance from the break far beyond anything observable. The folk claim is thus right in practice and wrong in principle, and the gap between the two is exactly the distance between measure zero and unobservably small.

Authors: Mark Jones, Jannik Schestag

Phylogenetic Diversity (PD) is a fundamental measure of biodiversity, originally defined on phylogenetic trees and widely used in conservation biology. Phylogenetic trees are often generalised to directed acyclic graphs, called phylogenetic networks. As such, a corresponding generalization of PD is needed. A natural generalization to edge-weighted phylogenetic networks is the all-paths measure, where the diversity of a set S of species (taxa) is defined as the total weight of all edges that lie on a path from the root to at least one species in S. While maximizing PD on trees can be solved in polynomial time, the corresponding problem on networks is NP-hard and difficult to approximate. We undertake a systematic parameterized complexity study of the Max-All-Paths-PD (MapPD) problem. We establish W[2]-hardness when parameterized by the number of species that are included in a solution, and W[1]-hardness for the number of species that are excluded. On the positive side, we show that the problem is fixed-parameter tractable with respect to the threshold of diversity and the acceptable loss of diversity. We further analyze how the network's proximity to a tree influences algorithmic behavior and present single-exponential fixed-parameter algorithms when parameterized by the number of reticulations and by the treewidth of the underlying graph. Finally, we present a polynomial kernelization for MapPD with respect to the number of reticulation edges.

Authors: Ken-ichi Kawarabayashi, Hirotaka Yoneda, Masataka Yoneda

We show that every planar graph with $n \geq 3$ vertices admits a 4-coloring in which each color is used on fewer than $n/2$ vertices. This bound is the best possible. Moreover, such a coloring can be found in $O(n \log n)$ time. We also extend these results to five or more colors and to graphs on general surfaces.

Authors: Stephan A. Fahrenkrog-Petersen, Aleksander Figiel, Darya Melnyk, Tijana Milentijević, Stefan Schmid

The stable marriage problem appears in many privacy-sensitive domains, for example in the National Resident Matching Program in the US. In such applications, preserving the privacy of users' preference lists is essential to prevent strategic manipulation, discourage misreporting, and comply with data protection regulations. In this work, we investigate privacy attacks on stable marriage algorithms. Assuming that the attacker (e.g., the hospitals) can repeatedly interact with the stable marriage algorithm, we demonstrate how such interactions can reveal private preferences of the non-malicious side (e.g., the residents). We show that the widely applied Gale-Shapley Matching Algorithm, where the proposers' side is malicious, is vulnerable to privacy attacks and all honest agents' preferences can be revealed. We further investigate which preference distributions of the honest, non-malicious side are susceptible to privacy attacks and show that the Gale-Shapley Matching Algorithm where the honest side proposes can preserve privacy in non-susceptible preference distributions. We extend our results to the decentralized setting and show that the attacker's side can infer all preference orderings. In an experimental evaluation, we test privacy attacks on synthetic and real-world data and show that real-world data is indeed susceptible to privacy attacks. This work underlines a need for new privacy-preserving stable marriage algorithms.

Authors: Elena Grigorescu, Shreya Nasa, Cameron Seth

The optimal sample complexity of testing if an $n$-vertex graph has an independent set of size $ρn$, or is $\varepsilon$-far from having an independent set of size $ρn$, was established to be $\widetilde{O}(ρ^3/\varepsilon^2)$, in a notable result by Blais and Seth (SICOMP 2025). In contrast, for $q$-uniform hypergraphs, there is a significant gap between the best known upper and lower bounds, and there has been no progress on the problem for the last two decades. In this work, we prove a new upper bound of $\widetilde{O}\!\left(\frac{qρ^{2q-3}}{\varepsilon^2 (q-2)!^2}\right)$ on the sample complexity of testing the $ρ$-independent set property. The previous best known upper bound was $\widetilde{O}\!\left(\frac{2^q q! ρ^{2q}}{\varepsilon^3}\right)$, due to Langberg (RANDOM 2004). This establishes the optimal dependence on $\varepsilon$ and gives an exponential improvement in the dependence on $q$. We prove our result via a new application of the hypergraph container method.

Authors: Sitan Chen, Weiyuan Gong, Qi Ye, Zhihan Zhang

Unlike in finite dimensions, quantum information in continuous-variable systems has the peculiar feature that without imposing physical constraints, the sample complexity of state tomography can be unbounded. Remarkably, this is even the case for state-of-the-art protocols for learning Gaussian states, which have finite-dimensional descriptions: the best known rates scale with $\log \log E$, where $E$ is the energy of the system. We prove this is not an artifact of existing analyses, but a fundamental limitation of the measurements used. We show: (1) Any protocol that uses Gaussian measurements, even entangled or adaptively chosen ones, must incur a $\log \log E$ dependence. This answers an open question posed by a number of previous works. (2) There is a smooth tradeoff between the number of rounds of adaptivity and the energy dependence, and we give a matching protocol achieving this interpolated rate. (3) With highly entangled, non-Gaussian measurements, one can learn $n$-mode pure Gaussian states with $O(n^2 / ε^2)$ samples, independent of $E$. This answers an open question posed by Chen et al. (4) A simple protocol based on the single-copy canonical phase POVM of Holevo and Helstrom learns single-mode pure Gaussian states with $O(1/ε^2)$ samples, again independent of $E$. Our results clarify the role of energy in bosonic state tomography and shed new light on the intriguing interplay between adaptivity, entanglement, and magic in quantum learning.

Authors: Siddharth Mitra, Vishwak Srinivasan, Xiuyuan Wang, Andre Wibisono

We show the Randomized Hamiltonian Monte Carlo (RHMC) algorithm has accelerated mixing time guarantees for sampling from log-concave probability distributions. RHMC proceeds by repeatedly simulating the continuous-time Hamiltonian dynamics for some random integration times, and resetting the velocity to be an independent Gaussian random variable between each simulation. We show that when the target distribution is log-concave and satisfies an $α$-Talagrand inequality (for example, if the target distribution is $α$-strongly log-concave), if we use a random integration time from either the triangular or the exponential distribution with mean $Θ(α^{-1/2})$, then RHMC converges exponentially fast in KL divergence, and the total integration time to reach error $\varepsilon$ in KL divergence scales as $O(α^{-1/2} \log(\varepsilon^{-1}))$. We also show that when the target distribution is log-concave, if we use a sequence of random integration times from the triangular distribution with exponentially increasing means, then the total integration time to reach error $\varepsilon$ in KL divergence scales as $O(\varepsilon^{-1/2})$. Our analysis relies on a bound on the average KL divergence along Hamiltonian dynamics, which is inspired by an analogous result on accelerated optimization methods based on Hamiltonian dynamics.

Authors: Peter Kiss, Mohammad Sharifi

We study the online preemptive matching problem, in which the edges of a graph arrive sequentially and the algorithm must maintain a matching by accepting or rejecting arriving edges and possibly discarding previously accepted ones. We prove a new upper bound of $0.5661$ on the competitive ratio achievable for the problem. This bound applies to arbitrary randomized algorithms, bipartite graphs and if we allow the algorithm to output a fractional solution. Our result improves upon the strongest previously known upper bound of $2-\sqrt{2} \approx 0.585$, due to Huang et al. [SODA'19]. Previous hardness constructions relied on edge sequences described by vertex arrivals where each arriving vertex reveals its edges to yet unvaried vertices. Under such sequences, Huang et al. showed that there exists a non-preemptive online algorithm with competitive ratio $\sim0.567$ (or $2-\sqrt{2}$ for fractional solutions). Consequently, our hardness construction is the first result which shows hardness for instances where the optimal algorithm employs preemption.

Authors: Moses Charikar, Jon Kleinberg, Chirag Pabbaraju

Motivated by the power of large language models, there has been renewed interest in the Gold-Angluin model of language identification in the limit, with an eye toward variants of the model that might overcome the negative results for its original formulation. Recent papers on this question have proposed looking at computational traces and annotations of training strings as a source of additional power for a learner, reflecting empirical regularities such as the way that commented source code is easier to learn from than arbitrary source code, and text annotated with algorithmically generated chain-of-thought tokens can be easier to learn from than the raw text itself. This recent work has shown positive results for language identification in the presence of such computational traces, but the traces in these positive results come from explicit automata-theoretic machine models that generate the language, where the underlying vocabulary of tokens for the traces is very large. In this paper, we address two fundamental issues left open by this line of work: can we achieve positive results with traces that use only a small alphabet, and can we define traces directly from the language itself, without requiring an underlying machine model that generates it? We establish positive results for both of these questions: for an arbitrary collection of languages, we show how to define computational traces that enable identification in the limit, using an alphabet of tokens that is linear in the size of the alphabet that the languages are defined over, and independent of any other properties of the languages.

Authors: David Aleman Espinosa, Sharat Ibrahimpur, Chaitanya Swamy

In recent years, much attention has been devoted to the study of optimization problems under norm-based objectives coming from the rich class of monotone, symmetric norms (and their generalizations). This work has however almost exclusively focused on covering problems, wherein one seeks to minimize the norm of the cost vector induced by a solution. We introduce and study the class of {\em norm-budgeted packing problems}, which are packing problems where the resource constraints underlying the packing problem are modeled via a {\em norm budget constraint} involving a {\em monotone, symmetric norm}. Formally, we have some elements with associated rewards and sizes, a downwards-closed collection of feasible solutions, and a budget $B$. Each solution induces a size vector, and the goal is to maximize the total reward subject to the norm-budget constraint $f(\text{size vector})\leq B$. The versatility of monotone, symmetric norms implies that a variety of classical packing problems can be captured under the umbrella of norm-budgeted packing problems. Moreover, the closure properties of monotone, symmetric norms, also enable one to encode multiple different norm-budget constraints via a single monotone, symmetric norm. We consider the norm-budgeted versions of a variety of canonical packing problems, including knapsack, matching, maximum-weight independent set in a $k$-set system, maximum generalized assignment problem (MaxGAP), and $k$-facility location, and develop a framework that allows us to obtain {\em constant-factor approximation guarantees} for these problems, and {\em PTASes for knapsack, and MaxGAP on identical and related machines}. We also develop constant-factor approximation algorithms for the {\em submodular} versions of some norm-budgeted packing problems, wherein the reward function is now specified by a monotone, submodular function.

Authors: Haripriya Pulyassary, Chaitanya Swamy

In $k$-clustering problems, we are given a metric space $(\mathcal{C}, d)$, and must choose a set $S$ of $k$ centers to open. Each client $j \in \mathcal{C}$ incurs an assignment cost, which is the distance between $j$ and center in $S$ that it has been assigned to. In this work, we study the \emph{minimum-norm $k$-clustering problem}, where we are given an arbitrary monotone symmetric norm $f$, and wish to open $k$ centers so as to minimize $f$(assignment-cost vector). This is a powerful generalization, encompassing many classical $k$-clustering problems including the $k$-median, $k$-means, and $k$-center problems. A simple and efficient algorithmic idea is that of \emph{adaptive sampling}, wherein we randomly choose the location of the next center to open with probability proportional to its ``cost" under the currently chosen set. While this has yielded fast algorithms for some $k$-clustering problem, little is known for settings \emph{without} ``min-sum" objectives. We devise the first adaptive-sampling-based bicriteria constant-factor approximation algorithm for general minimum-norm $k$-clustering, vastly expanding the scope of problems handled by adaptive sampling. For the special case of $\text{Top}_\ell$ norms, which form a building block of monotone symmetric norms, we show that adaptive sampling yields an $O(\log k)$-approximation algorithm.

Authors: Justin Dallant

Recent work of Ameli, Nederlof and Wang and of Dallant and Kozma introduced a framework for improving classical space--time tradeoffs for the Traveling Salesman Problem (TSP) and related permutation problems via extremal set systems with many maximal chains. In this note we observe that, for so called permutation problems whose outer aggregation is a minimum (such as TSP), the same framework admits a simple quantum analogue: instead of iterating over the covering family of set systems, we apply quantum minimum finding over the family. More precisely, let $P_S$ denote the optimal inverse normalized chain density among set systems of normalized size at most $S$. Then TSP admits a bounded-error quantum algorithm using $\widetilde O(S^n)$ QRAM space and \[ \widetilde O\!\left((S\sqrt{P_S})^n\right) \] time. The same argument applies to other minimization problems over permutations with a similar structure to TSP. Combining this observation with improved extremal set-system constructions of Andoni, Dallant, Kozma and Yu gives an explicit quantum space--time tradeoff curve, which beats the known quantum tradeoff by Caroppo et al. for all $1

Authors: Nima Anari, Aniket Das, Alireza Haqi

We study the parallel complexity of sampling from the high-temperature Ising mixed $p$-spin Gibbs measure, a canonical instance of a mean-field spin glass on the hypercube $\{\pm 1\}^n$. We propose two different algorithms for this problem, corresponding to two different regimes of accuracy. Our first algorithm is a parallel implementation of a Markov chain known as block dynamics, combined with an approximate rejection sampling step that uses an Ising model in a novel way as a proposal distribution to approximate the quadratic interaction terms of the $p$-spin Hamiltonian. For any $\varepsilon > 0$, this algorithm runs in $n^{\tfrac{1}{3}}\operatorname{polylog}(\tfrac{n}{\varepsilon})$ parallel time with $\operatorname{poly}(n, \log(\tfrac{1}{\varepsilon}))$ work, and outputs a sample whose law is $\varepsilon$-close to the $p$-spin measure in total variation distance. Our second algorithm uses Picard iterations to parallelize the Algorithmic Stochastic Localization (ASL) process of El Alaoui, Montanari, and Sellke (2025), and for any $\varepsilon > \varepsilon_n$, takes $\operatorname{polylog}(\tfrac{n}{\varepsilon})$ parallel time and $\operatorname{poly}(\tfrac{n}{\varepsilon})$ work to produce a sample that is $\varepsilon$-close to the $p$-spin measure in the normalized 2-Wasserstein metric. Here, $\varepsilon_n > 0$ is a threshold that goes to $0$ as $n \to \infty$. Our result constitutes a doubly exponential improvement in the $\varepsilon$ dependence of the runtime and an exponential improvement in the $\varepsilon$ dependence of the total work when compared to naïve ASL, whose runtime scales as $\exp(\operatorname{poly}(\tfrac{1}{\varepsilon}))$.

Authors: Maria Chudnovsky, Julien Codsi, Ajaykrishnan E S, Daniel Lokshtanov

An independent set in a graph $G$ is a set of pairwise non-adjacent vertices. A tree decomposition of $G$ is a pair $(T, χ)$ where $T$ is a tree and $χ: V(T) \rightarrow 2^{V(G)}$ is a function satisfying two axioms: for every edge $uv \in E(G)$ there is an $x \in V(T)$ such that $\{u,v\} \subseteq χ(x)$, and for every vertex $u \in V(G)$ the set $\{x \in V(T) | u \in χ(x)\}$ induces a non-empty and connected subtree of $T$. The sets $χ(x)$ for $x \in V(T)$ are called the bags of the tree decomposition. The tree-independence number of $G$ is the minimum taken over all tree decompositions of $G$ of the maximum size of an independent set of the graph induced by a bag of the decomposition. A graph $H$ is an induced minor of a graph $G$ if a graph isomorphic to $H$ can be obtained from $G$ by vertex deletions and edge contractions. We prove that for every $t\in\mathbb{N}$ there exists an $ε> 0$ such that every graph $G$ either contains the complete bipartite graph $K_{t,t}$ or the wall $W_{t\times t}$ as an induced minor, or has tree-independence at most $O(2^{O((\log n)^{1-ε})})$. This leads to algorithms with running time $2^{n^{o(1)}}$, for a wide range of problems on $\{K_{t,t}, W_{t\times t}\}$-induced minor free graphs. Our result is a substantial generalization of existing bounds for the tree-independence and tree-width on various graph classes, and a partial resolution of the conjecture of Chudnovsky, E S, and Lokshtanov [Arxiv, 2025] that $\{K_{t,t}, W_{t\times t}\}$-induced minor free graphs have poly-logarithmic tree independence number. The generality comes at the cost of a sub-polynomial, rather than poly-logarithmic upper bound. Our result leads to a complete classification of induced-minor closed classes into ones that have sub-polynomial tree-independence, tree-independence equal to $\tilde{O}(\sqrt{n})$, and linear tree-independence.

Front-running is a pervasive and costly problem on blockchains. Users earn rewards by publishing functional transactions that keeps markets efficient, such as arbitrage. But an attacker can observe such a transaction before it is ordered and publish her own ahead of it, seizing the reward and eroding users’ incentive to issue these transactions at all. The problem is well known and has drawn sustained effort from both industry and academia,...

By Hongyin Chen, Xu Zheng, Jichen Li, Ittay Eyal

Of the many cursed segments of vertical video, few annoy me more than the cottage industry of lunatics ranting over screenshots of PubMed pages to claim authority for whatever therapy, training program, or peptide they are selling. These people are all full of shit. They promise to cure your pain, get you thin, and make you stronger. All through SCIENCE.

These influencers will tell you that the scientific is better than the not scientific, and might even go so far as to say that the remainder is pseudoscience or quackery. But we unfortunately live in a narrow-minded world where scientific too often means “proven efficacious in a systematic review of randomized controlled trials.” You’ll be hard-pressed to find anything in the realm of treating musculoskeletal limitations that fits the bill.

It’s hard to shoehorn these sorts of therapies into the starting requirements for randomized trials. We’d need to start with clean definitions of an intervention and an outcome. What would these be for the management of pain by physical therapy?

Let’s start with the intervention. If we are being dogmatic evidence-based practitioners, the intervention in any physical therapy trial is the invitation to have therapy. According to the intention-to-treat principle, the invitation is the only thing that can be randomized. We can’t look only at the people who comply with every instruction and make it to the end of the rigorous therapy regimen. The model patients who diligently comply might differ from those who don’t, and our statistical signal will be biased if we include only the former. The only way to avoid these selection biases is to count everyone who was randomized, regardless of what happens between randomization and the final assessment. Both committed patients and no-shows contribute to the measured average efficacy.1

But what does it even mean for patients to follow the protocol? Physical therapy is far more complex than taking a drug. There’s no simple unit of treatment applied. Each interaction with a physical therapist involves a conversation about how things have been going, a plan for moving forward, some sort of interactive intervention in the office, and a discussion about what to do once the patient goes home. Each step here introduces a new branch in a deep decision tree. And every PT I’ve interacted with has been different, even when performing similar range-of-motion tests or manual therapies. Moreover, the treatment of any session depends on the entire history of the treatment so far. On top of this, every physical therapist I’ve seen has assigned daily exercises to do between sessions. This is part of the treatment, too! There is no way to perfectly isolate and randomize a single component of these complex treatment protocols.

A multi-stage protocol is an exponentially large collection of interventions. I made this point on the blog a few weeks ago in the context of anticoagulant trials for heart disease: “If you want to compare the effect of three different timings and three different dosages of a single drug, you need nine arms in your trial. If you want to additionally see if a second drug is helpful, you need 18.” Physical therapy is arguably much more complex.2

What about the outcome? In drug trials, we might get grim, unambiguous, objective outcomes like mortality. In vaccine trials, we might get unambiguous outcomes, such as a PCR diagnosis. Unfortunately, in pain management, the outcome is necessarily subjective. You can measure changes in range of pain-free motion, but there is too much heterogeneity to definitively stake out what a good outcome would be. Instead, pain therapies are most commonly evaluated based on improvements on the Numeric Ranking Scale. Studies ask participants at admission how their pain is on a scale of 0 (no pain) to 10 (worst pain imaginable). They ask them again at the follow-up. Statistical protocol then dictates computing the mean of the differences in treatment and control and running a t-test. You can try to remove the heterogeneity in how people respond to these questions, but these adjustments are based on subjective clinician calls. No matter what you do, pain is hard to mathematize. Doing statistics on these “numbers” and coming away with strong conclusions is a fool’s errand.

Beyond the treatment and outcome, all sorts of investigator biases make randomized trials even messier. You can blind the patients and the clinicians who assess outcomes, but you can’t blind the people applying physical therapy. It’s impossible to say what effect this sort of bias has on the scientific record. Even when well-intentioned, a clinician who believes in PT can subtly give away the secret assignment to their patients during a session.

I’ve never read a single study in this space that’s been compelling, and I don’t know why we hope that a narrow view of therapy can help us out of it. This fuels the fire of debate with people using studies to attack each other’s practices. There are countless articles and videos castigating stretching, massage, or cupping as not backed by evidence. These are all denounced as pseudoscience by a medical establishment that prescribed OxyContin like candy for two decades. Boy do I have some bad news for people who think there is great evidence that opioids work for pain management.

If we want to understand best practices for “wellness,” we need a different language around it. Maybe this language will need to lean on biomechanical plausibility or biochemical pathways. That certainly wouldn’t hurt. But more importantly, the language will have to prioritize discussions of craft, practice, and the cultivation of expertise. Whatever the case, the narrow definition of evidence-based needs to be reimagined. Healthcare is far more than a collection of unambiguous interventions with unambiguous outcomes.

Subscribe now

By Ben Recht

Authors: Luyining Gan, Jie Han, Bin Wang

A $4$-uniform $2$-cycle in a $4$-uniform hypergraph of length $t$ is a cyclic ordering of $2t$ vertices $v_1v_2\cdots v_{2t}v_1$ such that $v_{2i+1}v_{2i+2}v_{2i+3}v_{2i+4}$ are edges for $0\le i\le t-1$ while the addition is modulo $2t$. For every $γ>0$ and large $n$, we characterize the $n$-vertex $4$-uniform hypergraphs such that every triple of vertices is contained in at least $(1/3+γ)n$ edges and admits a Hamilton $2$-cycle. Up to the error term $γn$, the assumption on the minimum codegree is best possible and verifies a conjecture of Garbe and Mycroft. As a consequence, this gives a polynomial-time algorithm that decides whether an $n$-vertex $4$-uniform hypergraph with minimum codegree $(1/3+γ)n$ contains a Hamilton $2$-cycle. This stands as a steep contrast to the graph case where such a hardness gap has size $o(n)$.

Authors: Pingshi Yu, Chengsong Tan, Nicolas Wu, Alastair Donaldson

Randomised testing is a widely-used approach to software validation, yet its theoretical foundations remain thin. In particular, the fundamental question of what it means for a set of inputs to be \emph{generable} has gone unanswered in both the literature and folklore. We present the first complexity-theoretic foundations for random generators in software testing. We model generators as Turing transducers that consume random bits and produce string-encoded outputs, and show that the theoretically generable languages coincide exactly with the recursively enumerable languages. This has direct implications for testing at the boundaries of decidability, such as compiler testing. For \emph{efficient} generation, we show that the polynomial-time generable languages lie within \textit{NP}, that certain \textit{NP}-complete languages admit efficient generators, and that -- under standard cryptographic assumptions -- there are languages in \textit{P} for which no efficient generator exists: the complexity of efficienct generation and of efficient decision are not the same. We show space-bounded complexity is the natural framework for generators producing \emph{correlated} samples, capturing methodologies such as coverage-guided fuzzing and symbolic execution. Beyond classification, we characterise efficient generability: a language has a polynomial-time generator iff it admits a \emph{certificate scheme} over a verifier -- so witness planting, the folklore technique behind generators to test SAT solvers, is in a sense the only route to efficient generation. On the design of property-based testing libraries, we prove no library can compositionally derive efficient generators from logical predicates involving conjunction or negation, under standard assumptions. However, restricted classes like \textit{NL} (equivalently, linear Datalog predicates) would admit such a compilation.

Authors: Christoph Hertrich, Moritz Stargalla

We study tropical circuits with scalar multiplication gates, that is, algebraic circuits whose gates implement $\max$, $+$, or multiplication with a positive constant. For such circuits, we prove exponential size lower bounds for computing maximum weight directed spanning trees and maximum weight bipartite perfect matchings. As a corollary, we obtain an exponential size separation between monotone and non-monotone maxout neural networks, which generalize the popularly used ReLU neural networks. One conclusion from this is that neural network models with enforced convexity constraints, such as input-convex neural networks (ICNNs), sometimes need to be exponentially larger than their unrestricted counterparts in order to express the same functions.

Authors: Léo Gayral, Mathieu Hoyrup

Links between physicals systems and computability properties have been an active field of investigation in recent years. Inspired by a previous work in the context of positive temperature Gibbs measures, we prove here that in the context of zero-temperature ergodic optimisation, for a computable potential and provided with several reasonable assumptions, the maximum ergodic average is a computable real number, and the set of maximising measures is a $Π_1$-computable compact set. Then, in the more specific context of symbolic dynamics, with finite-range interactions on subshifts of finite type, we provide an explicit algorithm to compute both the maximum ergodic average and the set of maximising measures in finite time, with a matching code repository.

Authors: Chaithanya Rayudu, Jun Takahashi

The Lee-Yang theorem and its quantum extensions state that, for a broad class of Hamiltonians on any graph, the partition function's zeros in the complex magnetic field plane lie only on the imaginary axis. For these Hamiltonians, we prove that under a uniform Z-field of any strength h, the ground state has a spectral gap of at least h/4, independent of the system size and of the coupling strengths. The proof uses the zero-freeness of the partition function as given by Asano and Suzuki-Fisher to show exponential decay of the imaginary-time correlations for any product of Z-operators. Our result gives a polynomial-time quantum algorithm for computing the ground state energy of any Lee-Yang Hamiltonian.

Authors: Kai Cieliebak, Ángel González-Prieto, Eva Miranda

We prove that, given any contact $3$-manifold and any computable function $f: \mathbb{N} \dashrightarrow \mathbb{N}$, there exists a defining contact form and a Poincaré section of its Reeb flow whose partially defined return map computes $f$.

Authors: Hanlin Ren, Ryan Williams

We prove a near-maximum ($2^n / n$) circuit lower bound for the complexity class $\mathsf{E}^{\mathrm{pr}\mathsf{MA}}/_1$, corresponding to exponential time with access to a promise-$\mathsf{MA}$ oracle and one bit of advice. Our proof incorporates the iterative win-win paradigm (Chen--Lu--Oliveira--Ren--Santhanam, FOCS'23), the reduction from the Range Avoidance problem to circuit lower bounds (Jeřábek, Ann. Pure Appl. Log. '04; Korten, FOCS'21), and the PCP theorem. Crucial to our proof is the analysis of the complexity class $\mathsf{P}^\mathsf{NP}[{\textsf{#rounds}}=r, {\textsf{length}}=s]$, which is $\mathsf{P}^\mathsf{NP}$ with $r(n)$ adaptive rounds of $\mathsf{NP}$ queries, where each $\mathsf{NP}$ query has witness length $s(n)$.

Authors: Christoph Buchheim, Lowig Duer, Eva Ley, Maximilian Merkert, Komal Muluk

We consider bilevel optimization problems in which leader and follower jointly construct a feasible solution for an underlying combinatorial optimization problem. Response problems ask whether the leader can encourage -- or, in the pessimistic setting, enforce -- a reaction of the follower that includes a set of mandatory items while excluding a set of forbidden items. Our investigation focuses on tractability results for various cases which emerge from different combinations of the total number of mandatory, forbidden, and neutral items. After providing some results for response problems that hold for any underlying combinatorial optimization problem, we examine response problems over the maximum-weight matching problem and the minimum-weight perfect matching problem as illustrative and surprisingly varied examples. Among other results, we show that the response problem is hard for even a single given mandatory or forbidden edge. On the other hand, it is fixed-parameter tractable with respect to the total number of non-mandatory edges. If, however, each follower's edge is either mandatory or forbidden, the response problem for the perfect matching problem is solvable in polynomial time while it remains NP-hard for the maximum-weight matching problem.

Authors: Kagwe A. Muchane

Given two dense multivectors of the Clifford algebra $C\ell(V, Q)$ with $N=2^{p+q}$ coefficients, the fastest known classical algorithms compute their geometric product in $O(N^{ω/2})$ arithmetic operations, where $ω$ denotes the matrix multiplication exponent. I show that, under amplitude encoding, a quantum computer executes the geometric product in $O(\operatorname{polylog} N)$ time, using logarithmic space with sublogarithmic circuit depth. This exponential speedup establishes Clifford multiplication as a quantum primitive, providing an efficient computational foundation for quantum geometric algorithms and relativistic simulations.

Authors: Ei Ando

We investigate the problem of computing the distribution function for the shortest and longest path lengths in a directed graph with random edge lengths. Specifically, when these lengths are uniformly distributed, the problem reduces to computing the volume of a polytope defined by the graph structure. We establish that the problem is $\#P$-hard, even under the restricted condition that the random edge lengths are identically and independently distributed (i.i.d.) according to any continuous probability distribution with certain natural conditions, the local uniformity. This hardness result applies broadly: while the uniform distribution provides an essential case for the reduction, other distributions -- such as exponential or normal -- are similarly hard because they contain uniform distributions in every arbitrarily small interval. Furthermore, we show that the problem is contained within $\mathrm{XP}$ with respect to the treewidth $k$ of the underlying undirected graph. For the specific case of i.i.d. uniform edge lengths, we present a novel dynamic programming algorithm that processes a tree decomposition by iteratively performing convolutions to propagate distribution functions. Our approach achieves a time complexity of $n^{O(k^2)}$ for any fixed treewidth $k$.

Authors: Amirmahdi Mirfakhar, Maria-Florina Balcan, Hedyeh Beyhaghi

We study the problem of efficient online proportional sampling from a high-dimensional domain under a $σ$-smoothed adversary, where the sampling distribution is induced by a dynamically evolving weight function defined over a sequence of piecewise-structured partitions. This setting captures a broad range of applications, including principal-agent games (e.g., pricing and contract design), and algorithm configuration and parameter tuning. The central challenge is maintaining an efficient data structure as the induced partition grows increasingly complex over time -- naively, the number of subregions can grow as $O(t^d)$ by round $t$ in $d$ dimensions. We design a data structure that supports efficient updates and proportional sampling while avoiding the cost of explicitly maintaining this exponential growth, where the discontinuities are structured from axis-parallel hyperplanes. Under a $σ$-smoothed adaptive adversary, we prove a tight $O(\sqrt{σT})$ bound on the depth of our data structure, and an $O(\log T)$ bound under a random-order adversary -- to our knowledge, the first such results for this class of problems. We apply this framework to online learning with piecewise-structured rewards, obtaining efficient no-regret algorithms under both full-information and bandit feedback, with provable sublinear regret guarantees.

Authors: Sariel Har-Peled, Elfarouk Harb, Qizheng He

Given a positive integer $k$, we study the problem of finding a convex polygon of minimum perimeter that encloses exactly $k$ points of $\mathbf{Z}^2$. We show that an optimal polygon is contained in a circular annulus of width $O(k^{1/6})$, has $Θ(k^{1/3})$ boundary grid points, and its longest edge has length $Θ(k^{1/4})$. Using these structural bounds, we present a deterministic algorithm that computes an optimal polygon in $O(k^{29/18+o(1)})$ time, improving over the previous $O(k^3)$-time algorithm.

Authors: Chaya Keller, Shakhar Smorodinsky

The VC-dimension is a fundamental measure of the complexity of a set system. In this paper, we introduce and study a colorful variant of VC-dimension that captures the behavior of set systems on colored ground sets. By studying this new notion, we obtain a variety of geometric results. First, we prove that separable abstract convexity spaces with Radon number $D$ admit a Tverberg theorem with Tverberg number $O(D^2 r \log r)$. This bound significantly improves the $O(Dr^2\log r)$ bound of Alon and Smorodinsky from SODA'26 and is the first quasi-linear bound in $r$, in which the dependence on $D$ is not super-exponential. Second, we prove the first colorful $k$-wise Tverberg theorem for separable abstract convexity spaces. Using this theorem, we obtain a colorful selection lemma with $O(D^3)$ colors, an uncolored selection lemma for subsets of size $O(D^3)$, a weak $\varepsilon$-net theorem with nets of size $O_D(\varepsilon^{-O(D^3)})$, and a $(p,q)$-theorem with exponent of $\mathrm{poly}(D)$. All these quantitative bounds are significantly better than the best previously known general bounds for abstract convexity spaces. Finally, we extend our method to obtain a colorful Tverberg theorem for unions of convex sets, generalizing the uncolored theorem of Alon and Smorodinsky (SODA'26).

Authors: Oscar Chacón-Rivera

We study maximum-sum red-blue matchings and matching equilibrium for finite planar point sets. For a red-blue perfect matching $M = \{(a_i,b_i) : 1 \le i \le n\}$, we define the gain of a directed red cycle as the change in total weight produced by cyclically shifting the corresponding blue partners. We prove that $M$ is maximum-sum if and only if every directed red cycle has nonpositive gain, and we derive a geometric sufficient condition for optimality from cyclic intersections of distance-difference regions. We then characterize balanced matchings, in which all red-blue perfect matchings have the same total weight. Equilibrium is shown to be equivalent to vanishing cycle gains, to an additive form of the distance matrix, and to a common level-set condition for distance-difference functions. In the squared Euclidean case this yields an orthogonality classification, while in the Euclidean case it yields a hyperbolic level-set description and a collinear-separation classification in the nondegenerate setting.

Authors: Theodore Fung, Csaba D. Tóth

For an edge-weighted graph $G=(V,E)$ and a stretch parameter $t\geq 1$, a $t$-spanner is a subgraph $H\subseteq G$ such that the shortest path distances in $G$ and $H$ satisfy $δ_H(u,v)\leq t\, δ_G(u,v)$ for all $u,v\in V$. In metric spanners, $V$ is a finite metric space, and $G$ is the complete graph with edge weights corresponding to the distances between the endpoints. When $G$ is the complete graph on $n$ points in the plane, $O(n)$-size $t$-spanners are possible for any $t>1$: For every $\varepsilon>0$, there is an $(1+\varepsilon)$-spanner with $O(n/\varepsilon)$ edges (i.e., the stretch can be arbitrarily close to 1). When $G=K(R,B)$ is the complete bipartite graph on $n$ bichromatic points in the plane, in general, no spanner construction can guarantee stretch $t<3$ with $o(n^2)$ edges. Bose et al.~(SICOMP 2009) constructed a $(3+\varepsilon)$-spanner with $O(n\log n)$ edges for any constant $\varepsilon>0$. Our main result is a new construction for a $(3+\varepsilon)$-spanner with $O(\sqrt{1/\varepsilon}\cdot n)$ edges. Eliminating the $O(\log n)$ factor resolves a problem left open for more than 17 years, and raises a new research problem about optimizing the dependence on $\varepsilon$. We also study spanners for $G=K(R,B)$ on $n$ bichromatic points on the real line: In this case, we show that the MST of $K(R,B)$ is a 7-spanner, and we construct a 3-spanner with at most $2n-3$ edges.

Authors: Reilly Browne, Prahlad Narasimhan Kasthurirangan

We consider covering and partitioning a simple polygon into pieces which either have unit geodesic radius or unit geodesic diameter, using the $\ell_2$-metric for distances. There is no known method for finding an exact solution to these problems, even when the input size is constant, and the problem is known to be NP-hard in the case of polygons with holes. With this in mind, we instead devote our attention to developing simple approximation algorithms that run in polynomial time. For the radius problem, we present the first known approximation algorithms for both covering and partitioning, achieving a factor of 9. For the diameter problem, we are only able to give a positive result for the partition version of the problem, where we improve upon a complicated 72-approximation from Abrahamsen and Rasmussen [SODA '25], achieving a simple 15-approximation.

Authors: Lotte Blank, Anne Drieme, Sariel Har-Peled, Marena Richter

We give linear-time, and thus optimal, $(1+\varepsilon)$-approximation algorithms for numerous variants of the Frechet distance between $c$-packed curves (where $c \in O(1)$), removing an additional log factor that was present in previous algorithms. The key to our new algorithms is a linear-size approximation of the elevation function, which uses a decomposition of the domain into rectangles, and a careful implicit dynamic programming on this decomposition. The algorithm extends to the strong, weak, discrete, and continuous Frechet distances with a running time of roughly $O(cn/\varepsilon)$. The $c$-packedness assumption is used only in the analysis, and the algorithm is simple and should work efficiently for other inputs.

Authors: Marcin Kozik, Stanislav Živný

For a relational structure A, the Minimum Cost Constraint Satisfaction Problem is the following problem denoted by MinCostCSP(A): Given an instance of CSP(A) with rational costs on variable-value pairs, find a solution to the instance minimizing the sum of the chosen costs. For the exact minimization, a classification of MinCostCSP(A) in terms of A was established by Takhanov [STACS'10]. We focus on constant-factor approximations of MinCostCSP(A). DeHaan, Huang, and Lee recently showed that if A fails to admit a conservative near-unanimity polymorphism then MinCostCSP(A) is not constant-factor approximable [APPROX'25]. We provide a first step towards a classification, by proving a dichotomy for structures A admitting a conservative majority (also known as 3-near-unanimity) polymorphism. Our dichotomy criterion is not in terms of an algebraic condition on A but we show that this is unavoidable. We include a simple argument proving that no such condition exists.

Authors: Jason Fong, Jeffrey Kam, Steven Wong

The subgraph reconfiguration problem asks whether one subgraph can be transformed into another via a sequence of local changes while maintaining a specified graph property. In this work, we focus on the setting where the subgraph is specified by its set of edges. Our contributions in this paper are twofold. First, motivated by the contrast that path reconfiguration is $\textsf{NP}$-hard while tree reconfiguration is solvable in linear time, we prove two generalizations: (1) for any fixed $k$ at least one, reconfiguring connected graphs with pathwidth at most $k$ is $\textsf{NP}$-hard, and (2) for any fixed $k$ at least two, reconfiguring graphs with pathwidth at most $k$ is also $\textsf{NP}$-hard. En route to proving (2), we show a general hardness result that applies to a range of minor-closed graph classes, which we use to show planar graph reconfiguration is also $\textsf{NP}$-hard. Second, given our negative results, we extend the problem to a resource-focused setting, asking how much additional buffer space is needed to turn a non-reconfigurable instance into a reconfigurable one. We show that $Ω(n)$ extra buffer space is needed for planar graphs and graphs with bounded pathwidth and treewidth, while $O(1)$ extra buffer space is sufficient for cactus graphs in a restricted setting.

Authors: Yu Gao

Rectilinear matching to the integer grid asks to assign each of $n$ points in $\mathbb R^2$ to a distinct point of $\mathbb Z^2$, minimizing total $\ell_1$ movement. The main difficulty is that the target set is infinite: one must first identify a finite set of relevant grid points without losing optimality. We prove a geometric compression theorem for this infinite-target problem. In $O(n\log^2 n)$ time, we construct a set $\mathcal{C}$ of asymptotically optimal size $O(n)$ such that, simultaneously for every $p\in[1,\infty]$, some optimal $\ell_p$ assignment uses only points of $\mathcal{C}$. The construction is independent of the subsequent optimization algorithm and of the coordinate spread. For the rectilinear case, we combine this candidate set with a linear-size sparse network representation of $\ell_1$ distances. In the word-RAM model with $O(1)$-word dyadic coordinates and $O(\log n)$ fractional bits, a nearly-linear time minimum-cost flow algorithm then gives a randomized exact algorithm with expected running time $\widetilde O(n)$. This improves the standard $\widetilde O(n^2)$ approach. Combined with existing finite geometric matching algorithms, the same candidate set also gives an $\widetilde O(n\sqrt n\log(1/\varepsilon))$-time $(1+\varepsilon)$ approximation for every fixed integer $p\ge1$.

Authors: Fran Sučić, Leo Vitasović, Nikola Petrušić

We present Need for Speed Sort (NFS Sort), a recursive distribution-based sorting algorithm designed for numeric arrays. The algorithm partitions elements into equal-width value intervals, recursively refines dense buckets, and propagates analytical interval bounds between recursive calls, avoiding repeated scans for local minima and maxima. NFS Sort combines a fragment-based, cache-conscious scatter procedure for large subarrays with a lower-overhead auxiliary-array approach for smaller inputs. Small buckets are deferred to a final insertion-sort cleanup, while a comparison-based fallback is activated when recursive partitioning repeatedly fails to reduce the problem size. This mechanism guarantees a worst-case running time of O(n log n) and auxiliary space usage of O(log n). Experimental evaluation on synthetic inputs and real-world datasets from the SOSD benchmark suite compares NFS Sort with Balanced Learned Sort, IPS4o, Boost Spreadsort, PDQSort, and std::sort. The results show that NFS Sort is competitive or better than established state-of-the-art sorting methods across dataset sizes and distributions, outperforming the learned baseline particularly on smaller inputs while retaining strong performance at larger scales. Overall, NFS Sort combines efficient recursive distribution, practical memory management, and robust worst-case guarantees for high-performance numeric sorting.

Authors: Thiago Oliveira, Mohit Singh, Sahil Singla

We study online welfare maximization with divisible resources. A sequence of $n$ players arrive one by one; upon arrival, each player draws a valuation function over $m$ divisible items from a known distribution, reveals this valuation, and must be allocated an irrevocable fractional bundle subject to unit supply constraints. While online welfare maximization has been extensively studied for indivisible items and combinatorial valuations, much less is known when the resources are divisible and players have multi-dimensional concave valuations. We give approximation algorithms for monotone concave valuations satisfying diminishing returns. Our main result is a $2/3$-approximation to the optimal online policy, also known as the philosopher benchmark. The algorithm is guided by a low-dimensional concave relaxation of the online benchmark and rounds it via a new single-item capped online contention resolution scheme. This Capped-OCRS problem allocates to each realized type no more than its prescribed fractional bundle while preserving a $2/3$-fraction of that bundle in expectation. Its analysis uses a submartingale potential for the remaining side, we show that computing the optimal online policy is #P-hard even for a single divisible item. We also obtain a tight prophet inequality against the offline hindsight optimum. We show that a fixed-price auction with one linear per-unit price for each original divisible item achieves a $1/2$-approximation to the offline/prophet benchmark. The prices are obtained by aggregating Aumann--Shapley supporting prices, a continuous analogue of supporting prices for submodular/XOS set functions, and yield simple item prices rather than copy-dependent prices arising from discretization. The factor $1/2$ for the prophet benchmark is information-theoretically tight even for one item with linear valuations.

Authors: Aljaž Krpan, Janez Povh

We study the interaction between reduction rules and upper-bound functions for the Maximum Clique Problem (MCP). We show how MCP upper-bound functions can strengthen classical core and truss reductions by replacing local size conditions with upper-bound tests. This leads to the \((k,ω^u)\)-core, the \((k,ω^u)\)-truss, and the more general \((k,d,ω^u)\)-truss, where the parameter \(d\) controls the trade-off between stronger reductions and additional computational cost. For each of these notions, we prove clique-preservation properties, correctness of the corresponding peeling algorithm, and running-time bounds. Based on these reductions, we introduce a general framework for improving upper-bound values for MCP. We give two concrete instantiations of the framework: one that uses only the combined truss and core reductions, and one that combines the truss and core reductions with repeated applications of structions. Computational experiments on 73 benchmark graphs show that the proposed reductions can substantially improve several standard upper-bound functions and that combining multiple reduction methods can be beneficial in practice. In particular, the combination of structions, truss and core reductions with a DSatur-based bound often reached SDP-level upper-bound values faster than direct SDP computation; on the tested graphs with edge density below \(0.7\), it did so in every case. Using the truss and core reduction with the Lovász theta upper-bound function, we also improve the previously best certified integer upper-bound values for three difficult DIMACS instances whose exact clique numbers are not known. In particular, we improve upper-bound values for graph \texttt{C500.9} from 83 to 73, for graph \texttt{C1000.9} from 122 to 115, and for graph \texttt{C2000.9} from 177 to 168.

Authors: Alexander E. Black, Lei Xue

In 1992, Billera and Sturmfels introduced coherent monotone paths on polytopes as part of their description of the fiber polytope construction, and later in 1994 showed with Kapranov that these coherent monotone paths capture the topology of the space of all monotone paths, paths from a minimum to a maximum, in the directed graph of a polytope with orientation induced by a linear function. Those results motivate the following analog of the polynomial Hirsch conjecture: Does there always exist a coherent monotone path of polynomial length on a polytope for any choice of orientation induced by a linear function? We show this is not the case by exhibiting a family of polytopes and corresponding linear functions for which every coherent monotone path is exponentially long. As applications, we strengthen longstanding results pertaining to lower bounds for the shadow simplex method, geometric transversals in discrete geometry, and parametric linear optimization.

Authors: Moses Charikar, Jon Kleinberg, Chirag Pabbaraju

We study how little extra information is needed to make adversarial language learning possible. In Gold's model of language identification in the limit, a learner is given an enumeration of the strings from an unknown language chosen from a countable language collection. The learner guesses the identity of the language over the course of the enumeration, and it succeeds if, eventually, all of its guesses are the correct language. Classical results of Gold and Angluin show that many natural collections cannot be learned in this way. Recent work on trace colorings, motivated by the success of thinking-trace strategies in language learning, overcomes this obstruction by annotating every symbol of every string with a color. We ask whether the learner really needs this whole sequence of colors, or whether one color at the end of each string (a terminal coloring) is enough for language identification. We show that just one terminal bit per string is enough for every countable collection of infinite languages. In fact, the colorings can be chosen collection-independently: there is a single assignment of a two-color terminal coloring to every infinite language such that the same preassigned colorings identify every countable subcollection. Thus, in this model, an entire color trace can be compressed to one bit attached to the end of each example. Our global construction uses transfinite recursion, and we prove that this kind of nonconstructivity is unavoidable for any bounded number of colors. As a notion of constructivity, we use the formalism of Borel maps (a regularity condition satisfied by natural explicit constructions); we show that no global terminal coloring with a finite number of colors defined by a Borel map can identify all countable subcollections. By contrast, known trace-coloring constructions are Borel when encoded as terminal colorings, but require infinitely many colors.

Authors: Sayan Bhattacharya, Ermiya Farokhnejad, Thatchaphol Saranurak, Haoze Wang

The minimum degree spanning tree problem is a classic NP-hard problem whose optimal approximation guarantee was established since the early 1990s: Fürer and Raghavachari [FR92] gave an $\tilde O(mn)$-time algorithm that computes a spanning tree with maximum degree $Δ^\star+1$, where $Δ^\star$ denotes the optimum value. Whether similarly strong guarantees can be achieved in near-linear time has remained open for over three decades. We give the first near-linear-time algorithm that computes a spanning tree with maximum degree $\lceil (1+ε)Δ^\star\rceil+1$ in $\tilde O(m/ε^2)$ time. Prior near-linear-time algorithms either achieved the weaker bound $\lceil (1+ε)Δ^\star\rceil + O(\log n/ε^2)$ [DHZ20] or required dense graphs with $m\ge n^{7/4}$ [CQT21,BFW26]. Using the same framework, our algorithm can also compute a spanning tree with maximum degree $Δ^\star+1$ in $\tilde O(mn^{2/3})$ time, improving upon the recent $\tilde O(mn^{3/4})$-time algorithm of [BFW26]. These two results strictly improve all previous construction algorithms for the minimum degree spanning tree problem.

Authors: Tamas Visegrady

Montgomery's trick accelerates simultaneous modular inversion of $N$ inputs by amortizing a single shared inversion, but auxiliary multiplications for complement products are typically scheduled in a linear, serial form. We construct a maximally parallelizable data-flow graph (DFG) that computes all $\overline{x}$ complement~products by scheduling auxiliary multiplications into idle multiplier slots during accumulation of the product of all inputs, and that of the shared inversion. This scheduling ensures the post-inversion phase adds exactly one multiplication layer of latency regardless of $N$, yielding a critical path latency of $\lceil \log_2 N \rceil$ multiply layers, one inversion, and one final parallel multiply layer.

Authors: Salil Vadhan, Jiyu Zhang

Inspired by recent advances in notions of spectral approximation of digraphs [Ahm+20], we study spectral algorithms for finding periodic structures in digraphs via the spectrum of a class of rotated Laplacian matrices. This class of Laplacian matrices was previously studied by Lange, Liu, Peyerimhoff, and Post [Lan+15]. We consider a notion of periodicity ratio that generalizes the bipartiteness ratio of Trevisan [Tre09], and show that it is closely related to the spectrum of rotated Laplacian matrices. In particular, if the digraph is strongly connected and represents a Markov chain, this periodicity ratio for a given $p \in \mathbb{N}$ is a quantitative measure of how close this Markov chain is to having periodicity $p$. We propose and analyze a periodicity-ratio variant of the spectral algorithm by Louis, Raghavendra, Tetali and Vempala [Lou+12]. We show that the algorithm runs in randomized polynomial time and can find many nearly periodic components (i.e, components with small periodicity ratio). This also implies a new higher-order Cheeger-type inequality for periodicity in the spirit of that in [Lou+12; LOT14]. As part of our analysis, we prove a new theorem that upper bounds the probability that the largest magnitudes of two sequences of coordinate-wise correlated complex Gaussian random variables occur at different indices, which may be of independent interest. Previously, an analogous result was known only for real Gaussian random variables.

Authors: Alexandr Andoni, Justin Dallant, László Kozma, Hantao Yu

We nearly settle a natural extremal question about set systems over $[n]$: the tradeoff between the {size} (number of sets) and the number of {full chains}. This question was initially raised by Johnson, Leader, and Russell [Combin.~Probab.~Comp., 2015] as a counterpart to Sperner-type results in combinatorics. Recently, a framework introduced by Ameli, Nederlof, and Wang, and independently by Dallant and Kozma [FOCS 2026] linked this question to the space- and time-complexity of Bellman-Held-Karp-style dynamic programming algorithms for permutation problems such as the traveling salesman (TSP). Precisely, they showed that a space-time product $γ^{n+o(n)}$ is feasible for the TSP, whenever a set system of (normalized) size $S$ and chain density $D$ exists, with $ γ= S^2/D$. In this paper we show an essentially {optimal} bound of $γ\approx 3.1819$ for this quantity, closing the gap between the previous best lower and upper bounds of $γ\geq 3.015$ and $ γ\leq 3.572$ respectively. This implies a TSP algorithm with space-time product $O(3.1819^n)$ for input size $n$, as well as a limit to further improvements in this broad framework. More generally, we can obtain close to optimal values $D$ for any feasible value $S$, effectively settling the question of the number of full chains at every size. The crucial step towards our results is casting the extremal combinatorics question as an {information~vs.~entropy} tradeoff involving two random variables. This reformulation {exactly} captures the optimal tradeoff for the combinatorial problem, leading to a framework in which primal-dual certificates can be derived, proving rigorous upper and lower bounds on $γ$. We also give a further application of our techniques, improving a bound of Duffus, Sands, and Winkler on the minimum size of fibres in the Boolean lattice.