Power sampling improves large language model (LLM) reasoning at test time by sampling from a sharpened version of the base model's distribution. Its sampler, PowerMH, is slow because its Metropolis--Hastings (MH) updates are sequential: each accept/reject decision sets the state from which the next proposal is generated. Prior work reduces this cost by changing the sampler or the model (e.g., sequential Monte Carlo, approximate lookahead, distillation). We propose Predictive SubTree Prefetching (PreSTO), which accelerates PowerMH by prefetching future proposals and provably leaves its MH transition kernel unchanged. PreSTO views future MH outcomes as a tree of accept/reject decisions: in each model call, it generates the proposals of a prefetched subtree of this tree in one batch, then follows the realized path with the original MH acceptance rule. We also observe that MH reaches only a sparse, nearly unary part of each prefetched subtree, we further design a predictive traversal rule that chooses the subtree to maximize the predicted number of completed transitions. Across seven datasets and five base LLMs, PreSTO achieves median speedups of 1.71×, 1.91×, and 2.33× on PowerMH and its follow-ups EntropyCut and MultiTryMH, respectively; on the first two, it makes a median of 3.9× fewer model calls, each only 2.1× longer. On LCB V6 with Qwen3.5-9B, the predictive traversal rule is a further 1.25×--1.37× faster than the default breadth-first traversal.
Figures & tables
Figure 1: Across 105 sampler–LLM–dataset settings, PreSTO (ours) speeds up PowerMH, EntropyCut, and MultiTryMH by a median of 1.71× , 1.91× , and 2.33× , respectively. Each point is one setting, and faster settings lie toward the upper right.
Figure 2: Our subtree prefetching (PreSTO) accelerates PowerMH for transitions x→y→y . (a) PowerMH performs two sequential rounds of proposal generation with qτ followed by an MH accept/reject decision. (b) Our method collects feasible requests by traversing the prefetching tree, generates proposals (y,z) in one batched call, and then samples a path over the subtree.
Figure 3: The ϵ -rounded subtree (orange) is sparse and mostly has one child per node within the full prefetched subtree. Blue edges trace a sampled MH path.
Figure 4: Most edges in prefetched subtrees are ϵ -certain. For PreSTO-PowerMH with Qwen3.5-9B and the default BFS traversal, 67% – 95% of edges are ϵ -certain across all seven datasets, with ϵ=0.01 and Npf=20 .
(a) PowerMH vs. PreSTO-PowerMH
Base LLM
MATH500
AIME 24&25
MBPP
GPQA
HumanEval
MMLU
LCB V6
Qwen3-4B
2.11×
2.12×
1.10×
1.63×
1.77×
1.59×
2.04×
Qwen3-8B
2.23×
2.29×
1.00×
2.12×
1.81×
1.50×
2.05×
Qwen3.5-4B
1.49×
1.66×
1.67×
1.69×
1.59×
1.76×
2.02×
Qwen3.5-9B
2.06×
2.08×
1.61×
2.05×
1.21×
1.18×
1.44×
Gemma4-12B-it
1.74×
1.76×
1.44×
1.67×
1.71×
1.52×
1.95×
Table 1: Speedup of PreSTO over sequential samplers: the mean cumulative model-call time of the sampler divided by that of its PreSTO variant ( 20 ), on the first 20 prompts with K=100 and Npf=20 .
Figure 5: Under LCB V6 and Qwen3.5-9B, PreSTO ( Npf=20 , default BFS traversal) speeds up PowerMH by 1.44× in mean cumulative model-call time: it makes 4.92× fewer model calls (a ratio of means), each taking 3.42× as long, while the per-prompt transitions per call average 4.96 (Appendix D.2 ).
PowerMH
PreSTO-PowerMH, prefetch budget Npf
8
10
12
16
20
MH transitions per model call
1.00
3.59
3.94
4.27
4.63
4.91
Time per model call (s)
5.39
16.15
16.56
17.58
18.15
18.43
Speedup over PowerMH
1.00 ×
1.20 ×
1.28 ×
1.31 ×
1.38 ×
1.44 ×
Table 2: Under dataset LCB V6 and base LLM Qwen3.5-9B, raising the prefetch budget Npf from 8 to 20 increases the speedup of PreSTO-PowerMH over PowerMH from 1.20× to 1.44× . Transitions per call is K divided by the mean number of calls, and time per call is the total call time divided by the total number of calls.
Figure 6: Under base LLM Qwen3.5-4B across five datasets, PreSTO-PowerMH and PowerMH produce closely overlapping distributions of base-model log-likelihood and confidence, and the paired TOST finds them equivalent in 7 of the 10 comparisons (shaded).
Figure 7: Predictive traversal rule against BFS for PreSTO-PowerMH on LCB V6 with Qwen3.5-9B. (a–c) The predictive rule is 1.25× – 1.37× faster at every prefetch budget because it realizes 19% – 33% more MH transitions per call; points in (a, b) are ratios of totals over the 20 prompts, and error bars in (c) show ±1 SD over prompts. (d) First prefetched subtree of prompt 0, a favorable case: 3 transitions against 2 at every budget (prompts 1–3 in Appendix E.5 ). Label MH 0→k : the sampled path advances the chain from step 0 to k . Solid left edges: acceptance; dashed right edges: rejection; orange: ϵ -rounded subtree ( ϵ=0.01 ); blue: sampled path (endpoint double-bordered); gray: other prefetched branches.
Appendix figures & tables22 assets
Supplementary material from the paper’s appendix.
Appendix
Figure 8: Two batched suffix-resampling requests, their prefetching tree, and their shared KV cache. (a) Both requests reuse the cached prefix x0:c2 , and the request for y also reuses xc2+1:c1 . (b) Each acceptance (left) edge carries the MH acceptance probability computed from the pair of sequences it connects in (a). (c) Radix-tree representation of the shared KV-cache segments.
Figure 9: Visualization of ( 13 ): (a) certain rejection, (b) certain acceptance, and (c) an uncertain rounded decision. A ✗ marks a pruned edge. The proposal at v+ is required in every case.
Figure 10: Budget decomposition visualization for Proposition 2 .
Category
Dataset
Source (Hugging Face)
Instances
Mathematical reasoning
MATH500
HuggingFaceH4/MATH-500
500
AIME 24&25
Maxwell-Jia/AIME_2024 MathArena/aime_2025
45
Code generation
HumanEval
openai/openai_humaneval
164
MBPP
google-research-datasets/mbpp
257
LCB V6
livecodebench/code_generation_lite
175
Scientific and general knowledge
GPQA
Idavidrein/gpqa
198
Appendix
Table 3: Datasets and their full sizes; every run uses only the first 20 prompts of each dataset.
Family
Base LLM
Source (Hugging Face)
Released
Qwen3
Qwen3-4B
Qwen/Qwen3-4B
Apr 2025
Qwen3-8B
Qwen/Qwen3-8B
Apr 2025
Qwen3.5
Qwen3.5-4B
Qwen/Qwen3.5-4B
Feb 2026
Qwen3.5-9B
Qwen/Qwen3.5-9B
Feb 2026
Gemma4
Gemma4-12B-it
google/gemma-4-12B-it
May 2026
Appendix
Table 4: Base LLMs and their sources.
Software
Version / configuration
Python
3.13.0
PyTorch
2.13.0
vLLM
0.27.1 (V1 engine) with custom engine patches
Transformers
5.14.1
(a) Software configuration
Hardware
Specification
Appendix
Table 5: Software and hardware configuration.
Parameter
Setting
Response length
n=1024 tokens (fixed; no early stop at EOS)
Block size
a single block, B=n=1024
Power exponent
α=4.0
Proposal LLM
qτ= base model sampled at temperature τ
Proposal temperature
τ=1/α=0.25
Realized MH transitions
K=100
Appendix
Table 6: Hyperparameter configurations.
Figure 11: PowerMH and PreSTO-PowerMH on MATH500.
Dataset
Base LLM
Total Time (s) ↓
Transitions per call ↑
Time/call (s) ↓
Speedup
MATH500 (Fig. 12 )
Qwen3-4B
277.41±45.48/157.85±20.51
1/4.26±0.54
2.77±0.45/6.62±0.23
1.76×
Qwen3-8B
336.64±58.87/192.00±24.30
1/4.34±0.47
3.37±0.59/8.23±0.27
1.75×
Qwen3.5-4B
365.78±80.61/233.47±73.77
1/3.81±0.99
3.66±0.81/8.22±0.47
1.57×
Qwen3.5-9B
478.02±74.36/272.29±59.03
1/3.84±0.65
4.78±0.74/10.11±0.52
1.76×
Gemma4-12B-it
980.09±240.16/431.85±175.94
1/4.07±2.00
9.80±2.40/15.06±1.37
2.27×
AIME 24&25
Qwen3-4B
255.04±37.63/150.54±16.85
1/4.57±0.49
2.55±0.38/6.81±0.18
1.69×
Appendix
Table 8: EntropyCut execution metrics across base LLMs.
Figure 12: EntropyCut and PreSTO-EntropyCut on MATH500.
Dataset
Base LLM
Total Time (s) ↓
Transitions per call ↑
Time/call (s) ↓
Speedup
MATH500
Qwen3-4B
741.55±35.32/322.90±17.49
1/2.02±0.11
7.42±0.35/6.49±0.17
2.30×
Qwen3-8B
865.84±42.90/417.08±21.21
1/2.04±0.10
8.66±0.43/8.50±0.18
2.08×
Qwen3.5-4B
2174.72±142.35/375.23±14.43
1/2.05±0.08
21.75±1.42/7.68±0.16
5.80×
Qwen3.5-9B
3246.66±262.85/463.26±25.28
1/2.05±0.11
32.47±2.63/9.48±0.15
7.01×
Gemma4-12B-it
1558.70±79.33/686.49±45.49
1/1.86±0.13
15.59±0.79/12.75±0.20
2.27×
AIME 24&25
Qwen3-4B
633.31±31.36/313.16±17.48
1/2.03±0.11
6.33±0.31/6.33±0.18
2.02×
Appendix
Table 9: MultiTryMH execution metrics across base LLMs.
Figure 13: Per-transition traces across prefetch budgets on LCB V6 with Qwen3.5-9B. Solid curves show the PreSTO mean for each budget, from light (small Npf ) to dark (large Npf ). Dashed curves show PowerMH, and shaded bands show ±1 SD over traces.
Figure 14: Likelihood and confidence distributions with and without PreSTO. Boxes show the per-prompt token-normalized log-likelihood (a, c) and confidence (b, d) of terminal traces under p0 for the eight PowerMH pairs (top) and the eleven EntropyCut pairs (bottom), each against its PreSTO variant. Shaded rows are equivalent under the paired TOST ( peq<0.05 ).
Figure 15: Edge certainty under PreSTO-PowerMH across models and datasets.
Figure 16: Edge certainty under PreSTO-EntropyCut across models and datasets.
Model calls
Transitions per call
Time per call (s)
Model-call time (s)
Npf
BFS
Predictive
BFS
Predictive
BFS
Predictive
BFS
Predictive
Speedup
8
27.85
23.40
3.59
4.27
16.15
15.35
449.7±19.1
359.2±28.0
1.25×
10
25.40
20.15
3.94
4.96
16.56
16.00
420.7±25.5
322.3±28.5
1.31×
12
23.40
18.90
4.27
5.29
17.58
17.13
411.4±38.7
323.8±39.5
1.27×
16
21.60
17.10
4.63
5.85
18.15
17.34
392.1±33.8
296.6±45.8
1.32×
20
20.35
15.30
4.91
6.54
18.43
17.95
375.1±35.1
274.7±40.1
1.37×
Appendix
Table 10: Predictive traversal rule against the default BFS traversal. Speedup is the total model-call time of BFS divided by that of the predictive rule.
Figure 17: First prefetched subtree under each traversal rule and prefetch budget. Prompts 1–3 of LCB V6 with Qwen3.5-9B under PreSTO-PowerMH, one row per prompt and rule; labels and colors follow Figure 7 d.
Engine
Sampler
Calls
Time per call (s)
Model-call time (s)
Speedup
vLLM
PowerMH
100
5.39±0.31
539.24±30.92
1.40×
PreSTO-PowerMH
20.60±1.31
18.75±0.43
385.97±21.25
SGLang
PowerMH
100
3.57±0.24
357.40±23.67
1.97×
PreSTO-PowerMH
21.25±1.86
8.56±0.10
181.74±15.25
Appendix
Table 11: PowerMH and PreSTO-PowerMH on vLLM and SGLang.
Figure 18: Prediction of the ϵ -certainty class by feature set. Held-out accuracy against log loss of the logistic regression for PreSTO-PowerMH (a) and PreSTO-EntropyCut (b); upper left is better. Labels give the feature sets used: a single digit is one set, “1+3” combines Features 1 and 3, and “1–3” combines Features 1 through 3; comma-separated labels mark overlapping points. The cross marks the constant baseline, and the bold label is the configuration we use.
Figure 19: Effect of hidden layers on predicting the ϵ -certainty class. Held-out accuracy against log loss with the cut fraction as the only feature; labels give the units per hidden layer, and “none” is the logistic regression. Markers as in Figure 18 .
Figure 20: Speed and memory against test-time baselines. Peak KV-pool use against time per prompt on the first 20 prompts of LCB V6 with Qwen3.5-9B; lower left is better.
Power-sharpened sampling is an inference-time alternative to reinforcement-learning (RL) post-training for enhancing reasoning in large language models (LLMs). High-probability sequences are amplified under the base model without parameter updates or external rewards, avoiding the costly optimization and jagged generalization of RL. However, this approach faces a fundamental exploration--exploitation trade-off, as % strong sharpening restricts exploration, trapping samplers in plausible but incorrect reasoning trajectories, whereas weak sharpening leaves the answer distribution diffuse. To resolve this trade-off, we introduce Parallel Power Tempering (PPT), instantiating power-sharpened LLM sampling via parallel tempering. Running multiple interacting replicas in parallel at different sharpening levels allows lower-power replicas to explore diverse reasoning trajectories and higher-power chains to further exploit higher-likelihood responses favored by the sharpened target. Specifically, we tailor \method{} to inference-time sampling by mitigating a truncation bias, identified in prior power samplers, and investigate effective swap strategies under finite memory and compute budgets. Extensive experimentation shows that \method{} substantially improves single-chain power-sharpened sampling and outperforms RL-post-trained models, producing higher-quality reasoning traces and even achieving performance comparable to frontier models.
Panagiotis Theodoropoulos, Nan Jiang, Xintong Duan +4
Georgia Institute of Technology · Morgan Stanley · University of Texas at El Paso +1
Reinforcement learning (RL) post-training is a dominant approach for improving the reasoning performance of large language models (LLMs), yet growing evidence suggests that its gains arise primarily from distribution sharpening rather than the acquisition of new capabilities. Recent work has shown that sampling from the power distribution of LLMs using Markov chain Monte Carlo (MCMC) can recover performance comparable to RL post-training without relying on external rewards; however, the high computational cost of MCMC makes such approaches impractical for widespread adoption. In this work, we propose a theoretically grounded alternative that eliminates the need for iterative MCMC. We derive a token-level decomposition of the autoregressive power distribution, showing that the global power distribution can be approximated by a token-level scaled low-temperature one, where the scaling factor captures future trajectory quality. Leveraging this insight, we introduce a training-free and verifier-free algorithm that sharpens the base model's generative distribution autoregressively. Empirically, we evaluate our method on math, QA, and code tasks, and show that our method matches or surpasses one-shot GRPO without relying on any external rewards, while reducing inference latency by over 10× compared to MCMC-based sampling.
Xiaotong Ji, Rasul Tutunov, Matthieu Zimmer +1
Huawei Noah’s Ark Lab · UCL Centre for Artificial Intelligence
Sampling from the sequence-level power distribution pα elicits RL-level reasoning from base language models without any parameter updates, but the standard Metropolis--Hastings (MH), a Markov Chain Monte Carlo (MCMC) sampler, is both expensive and slow-mixing. We trace both to a structural mismatch: pα mainly departs from p at a sparse, spatially clustered set of high-entropy decision points, yet MH proposes resampling positions uniformly along the prefix -- wasting compute on near-degenerate conditionals while under-mixing precisely where modes diverge. We propose Entropy-Guided Power Sampling (EGPS), a training-free and verifier-free sampler that re-derives its proposal from token-level entropy already in the forward pass. EGPS skips deterministic blocks, localizes each MCMC move to a high-entropy neighborhood, and applies Multiple-Try Metropolis at decision points -- making sampling cost scale with entropy mass rather than sequence length. On Qwen2.5-Math-7B, EGPS reaches best or tied-best accuracy on all three benchmarks (MATH500 75.8%, HumanEval 62.2%, GPQA 42.4%) at up to a 12.6× wall-clock speedup over the MH baseline.
Hong Guo, Nianhui Guo, Christoph Meinel +1
Hasso Plattner Institut · German University of Digital Science