cs.NEJun 11, 2026

Runtime Analysis of the (μ+1)(μ+ 1)-ES in a Homogenous Progress Model

Authors: Johannes LenglerRaghu Raman Ravi

Abstract

We introduce a new simple model to study the fitness progress of Evolution Strategies (ES) in generic problems. In this model, we bypass the underlying fitness landscape and assume that the mutation of any individual produces an offspring whose fitness relative to the parent is given by an invariant distribution ZZ, such as a mean-shifted Gaussian. This serves as a prototypical model for the optimisation landscape when an evolution algorithm operates far from the global optimum. This simple model can be used to approximate the optimisation process for problems where it is intractable to model the exact fitness function, including tasks such as hyperparameter tuning in machine learning models. We rigorously analyse the expected growth rate Rμ\mathcal{R}_μ of the continuous steady-state (μ+1)(μ+1)-ES in this model. Unlike comma-selection strategies, the steady-state (μ+1)(μ+1)-ES maintains overlapping generations, introducing complex mathematical dependencies among surviving parents that make it harder to analyse. We give a general technique to analyse the the (μ+1)(μ+ 1)-ES by constructing modified processes whose growth rates provably sandwich that of the original process. These modified processes are then easier to analyse but still close enough to the true process to give a tight bound on the expected growth rate. When Z=N(δ,1)Z = \mathcal{N}(-δ, 1) and μeδμ\le e^δ, we show that Rμ=log1+o(1)μμR1\mathcal{R}_μ = \frac{\log^{1 + o(1)} μ}μ \mathcal{R}_1.

Explore similar work

Sep 14, 2026cs.NE

Signed Sensitivity of Expected Hitting Time to Mutation Rate in the (1+1) EA: Per-State Sign Theorems and Verifiable Certificates for Non-Lumpable Families

For the (1+1) evolutionary algorithm with standard bit mutation, we study the sensitivity of the expected hitting time Hp=ExTH_p=\mathbb{E}_x T to the mutation rate. We first point out an easily overlooked formalization pitfall: the improvement event is not monotone in the mutation mask, so the unsigned (total-influence) form of the Margulis-Russo formula does not apply; the correct object is the signed endpoint difference. Second, we give an exact three-dimensional separation: two fitness functions share the entire one-step success-rate curve, yet their expected hitting times are two different exact rational numbers; hence one-step success-rate quantities do not determine the expected hitting time. Building on the runtime derivative Hp=(IQp)1QpHpH'_p=(I-Q_p)^{-1}Q'_p H_p, we construct computable double-residual sign certificates, prove a per-initial-state sign theorem on OneMax (for every non-optimal initial state, cH<0\partial_c H<0 on 0<c<10<c<1, where p=c/np=c/n; at c=1c=1 only the distance-one state is stationary), and extend the framework to non-lumpable positive linear families: an explicit non-lumpability witness, a block-interval double-residual certificate that covers all states without enumerating them, a uniform sign bound cET9n/16\partial_c\mathbb{E}T\le -9n/16 over the whole interval c[1/4,1/2]c\in[1/4,1/2] for an explicit family at all even scales n8n\ge 8, and a heterogeneous instance certificate Hx1/6H'_x\le -1/6 on 57 of 63 states across c=1c=1. All finite verifications use exact rational arithmetic. A bounded systematic literature search did not uncover this exact combination, although the underlying tools are well established; we therefore make no novelty claim beyond the stated combination.
RenKai Wang
Jun 11, 2026cs.NE

The (1+1)(1 + 1)-EA in Dynamic Environments

We study the (1+1)(1 + 1)-EA in dynamic linear environments, where in every generation selection is performed with respect to a freshly sampled linear function with positive weights. We consider the Dynamic Binary Value problem, where each generation uses a uniformly random permutation of 1,2,4,,2n11,2,4,\dots,2^{n-1}, and a Uniform weight variant, where the weights are drawn independently from Unif(0,1)\mathrm{Unif}(0,1). Both of them have recently been integrated into the IOHprofiler platform and empirically studied. For both models we prove a sharp threshold in the mutation parameter χχ for mutation rate χ/nχ/n. Below the threshold, the expected optimisation time is O(nlogn)\mathcal{O}(n\log n), whereas above it the runtime becomes 2Ω(n)2^{Ω(n)}. For the Dynamic Binary Value problem in the exponential regime, we also quantify at what distance from the optimum the optimisation process stagnates. We show that there is a second threshold: a distance that is efficiently reached, but reaching any smaller distance takes exponential time. This quantifies and proves previous empirical findings.
Georg Hasebe, Johannes Lengler, Raghu Raman Ravi
May 20, 2026cs.NE

Convergence Analysis of Evolution Strategies for Mixed-Integer Optimization

Mixed-integer extensions of evolution strategies (ES) that discretize selected coordinates of sampled continuous vectors often impose a lower bound on the standard deviation of integer variables to prevent premature convergence. While these methods show promising empirical results, this handling can slow the convergence of continuous variables, and its impact has lacked a clear theoretical account. In this paper, we provide a convergence analysis of evolution strategies for mixed-integer optimization, inspired by the drift analysis of the (1+1)-ES in the continuous domain. Specifically, we consider two (1+1)-ES variants for mixed-integer domains: (1+1)-LB-ES, which introduces a lower bound on the standard deviation for integer variables, and (1+1)-LUB-ES, which combines both lower and upper bounds to enhance the convergence of the continuous variables. Focusing on the optimization phase after the integer variables have been optimized, we rigorously analyze their convergence behavior on a benchmark function designed for mixed-integer domains. Our results show that (1+1)-LB-ES can suffer from premature convergence when the number of integer variables is large, while (1+1)-LUB-ES achieves linear convergence under suitable parameter settings. These findings provide theoretical insights into the impact of integer handling on convergence performance and guidance for the design of mixed-integer ES.
Ryoki Hamano, Kento Uchida, Shinichi Shirakawa