Best-of-n (BoN) sampling is a simple yet effective inference-time alignment method, but hard maximization provides only coarse control over the trade-off between reward and distribution shift. Soft Best-of-n (Verdun et al. 2025) provides smoother control and converges to the optimal distribution associated with KL-regularized reward maximization. In this paper, we introduce ExpBoN, an alternative soft BoN method based on the exponential-noise report-noisy-max mechanism. It admits an exact finite-n decomposition, which yields exponentially fast convergence in total variation, expected reward, and both directions of KL divergence. We provide comprehensive theoretical analyses of its convergence and regret behavior. We further integrate ExpBoN into the guided speculative inference (GSI) framework (Geuter, Mroueh, and AlvarezMelis 2025), resulting in ExpGSI, for efficient reward-guided LLM alignment. ExpGSI yields substantial reductions in computational cost while maintaining comparable accuracy. Experiments on MATH500, MMLU-STEM, and Minerva Math with the Qwen2.5-Math and Qwen3 model families show that ExpGSI reduces estimated computation by 14%-39% across candidate budgets for Qwen2.5-Math and by up to 45% at n=16 for Qwen3. Overall, our results provide a theoretical and algorithmic foundation for exponential-noise BoN and efficient test-time LLM alignment.
Figures & tables
Figure 1: Finite- n comparison of ExpBoN and SBoN on the example P=(0.75,0.20,0.05) . Panels (a) and (b) use r=(0.016,0.164,0.820) and λ=0.5 . Panel (a) compares the KL divergences with the ExpBoN bounds in Theorem 3.2 and the O(1/n) SBoN bound. Panel (b) compares the relative centered-reward gaps with the ExpBoN bound ρλn from Corollary 3.4 and the corresponding O(1/n) SBoN bound. Panel (c) uses r as the true reward and rp=r+0.1(1,0,−1) as the proxy reward at λ=0.5 , and compares the exact finite- n true-reward losses from the common proxy-tilted limit with the ExpBoN bound in Theorem 3.9 and the SBoN bound. The results show the exponential rates of ExpBoN in contrast to SBoN’s polynomial bounds.
Figure 2: Compute–accuracy planes at each candidate budget n (Qwen2.5-Math; macro average over MATH500, MMLU-STEM, and Minerva Math; 3 -seed means). Arrows mark the GSI → ExpGSI reduction in estimated compute. ExpGSI maintains accuracy comparable to GSI at every n while its compute saving widens monotonically with n , from 14% at n=2 to 39% at n=16 ; full metrics for both model families are reported in Table 1 .
ri←r(h,yi),ℓB,i←logπB(yi∣h)
Algorithm 1 ExpBoN-based GSI
Model family
n
Method
Acc. (%)
Time/step (s)
Accept. (%)
Est. TFLOPs/prob.
Qwen2.5-Math(1.5b/7b/7b)
4
S-BoN ( πS )
55.9±0.5
0.49±0.01
–
252±6
RSD †
57.3±1.2
0.58±0.01
95.7±0.2
257±7
S-BoN ( πB )
62.1±0.4
1.06±0.01
–
297±6
GSI
59.9±0.6
0.87±0.01
82.1±0.2
370±5
ExpGSI (ours)
59.9±0.8
0.84±0.02
81.9±0.4
285±4
16
S-BoN ( πS )
58.1±0.2
1.05±0.01
–
1182±16
Table 1: Macro-average results over MATH500, MMLU-STEM, and Minerva Math and multiple seeds. ExpGSI maintains accuracy and acceptance comparable to GSI while reducing time and computation. At n=16 it reduces time and TFLOPs by 18% and 39% for Qwen2.5-Math, and by 30% and 45% for Qwen3. Bold values indicate ExpGSI efficiency improvements over GSI. † RSD follows the GSI implementation and hyperparameter configuration.
n
Method
Acc. (%)
Time/step (s)
Accept. (%)
Est. TFLOPs/prob.
2
S-BoN ( πS )
76.5±0.4
0.40±0.00
–
121±5
RSD †
79.3±1.9
0.49±0.02
96.7±0.3
118±3
S-BoN ( πB )
83.2±1.1
0.91±0.00
–
148±1
GSI
80.3±1.2
0.66±0.02
89.6±0.5
174±5
ExpGSI (ours)
80.2±0.7
0.61±0.02
89.5±0.8
147±3
4
S-BoN ( πS )
78.1±0.8
0.52±0.03
–
244±22
Table 2: Complete results on MATH500 across candidate budgets. Reported values are means and standard deviations over three random seeds. Bold efficiency values indicate improvements of ExpGSI over GSI. † RSD uses the implementation and hyperparameter configuration adopted by GSI.
n
Method
Acc. (%)
Time/step (s)
Accept. (%)
Est. TFLOPs/prob.
2
S-BoN ( πS )
59.2±2.1
0.35±0.00
–
107±0
RSD †
67.1±0.8
0.48±0.01
89.4±0.4
92±2
S-BoN ( πB )
68.8±2.3
0.79±0.01
–
100±1
GSI
68.4±1.5
0.81±0.01
62.4±1.0
133±0
ExpGSI (ours)
67.2±1.2
0.77±0.01
62.6±1.0
110±3
4
S-BoN ( πS )
64.0±0.7
0.44±0.01
–
212±9
Table 3: Complete results on MMLU-STEM across candidate budgets. Reported values are means and standard deviations over three random seeds. Bold efficiency values indicate improvements of ExpGSI over GSI. † RSD uses the implementation and hyperparameter configuration adopted by GSI.
n
Method
Acc. (%)
Time/step (s)
Accept. (%)
Est. TFLOPs/prob.
2
S-BoN ( πS )
25.0±1.3
0.39±0.01
–
149±5
RSD †
24.9±1.7
0.50±0.02
94.1±0.2
150±4
S-BoN ( πB )
31.9±1.5
0.90±0.02
–
171±5
GSI
28.6±1.1
0.77±0.01
78.1±0.3
198±6
ExpGSI (ours)
29.7±0.6
0.73±0.02
78.5±0.8
175±5
4
S-BoN ( πS )
25.7±0.4
0.51±0.01
–
301±10
Table 4: Complete results on Minerva Math across candidate budgets. Reported values are means and standard deviations over three random seeds. Bold efficiency values indicate improvements of ExpGSI over GSI. † RSD uses the implementation and hyperparameter configuration adopted by GSI.
Benchmark
n
Method
Acc. (%)
Time/step (s)
Accept. (%)
Est. TFLOPs/prob.
MATH500
4
S-BoN ( πS )
67.6±0.2
0.27±0.00
–
714±13
RSD †
67.2±0.0
0.28±0.00
99.5±0.0
715±5
S-BoN ( πB )
74.0±0.4
0.75±0.00
–
802±1
GSI
68.2±0.7
0.46±0.01
98.4±0.2
1455±1
ExpGSI (ours)
67.9±0.2
0.40±0.01
98.2±0.1
1006±7
16
S-BoN ( πS )
68.4±0.5
0.70±0.00
–
2874±14
Table 5: Complete per-benchmark results for the Qwen3 configuration ( Qwen3-1.7B draft, Qwen3-14B target, thinking mode disabled, same PRM and clipping level as the main experiments) at n∈{4,16} . Reported values are means and standard deviations over two random seeds. Bold efficiency values indicate improvements of ExpGSI over GSI. † RSD uses the implementation and hyperparameter configuration adopted by GSI.
n
Method
Acc. (%)
Time (s)
Accpt. (%)
TFLOPs
4
S-BoN ( πS )
76.5±0.0
.44±.00
–
41±1
RSD †
79.2±1.4
.93±.02
65.5±.8
43±2
S-BoN ( πB )
82.7±0.9
1.09±.02
–
76±0
GSI
83.6±1.0
.98±.01
73.0±1.2
141±3
ExpGSI (ours)
83.8±0.5
1.01±.03
72.2±.7
117±2
16
S-BoN ( πS )
75.1±1.7
.67±.03
–
148±6
Table 6: Results with the Skywork reward model on MATH500. Values are means and standard deviations over three seeds.
C/C
Acc. (%)
Scored cand./step
Est. TFLOPs/prob.
0.5
82.4±1.9
8.9
892
0.75
83.0±0.4
8.8
892
1.0 (default)
82.9±0.9
9.2
941
1.5
82.3±0.7
9.6
951
2.0
80.8±1.4
10.2
993
∞ (no clipping)
81.9±0.5
16.0
1592
Table 7: Sensitivity of ExpGSI to the clipping level on MATH500 with n=16 . Accuracy is reported as mean ± standard deviation over two seeds. When C=∞ , clipping and early exit are disabled.
Figure 3: Exact finite- n convergence of ExpBoN vs. SBoN on real LLM candidate pools (medians with 10–90% bands). Columns: the reward score r/λ of Section 2 at the deployed temperature λ=1/20 and at λ=1/2 , the temperature of the synthetic example in Figure 1; the GSI score of Section 4 without clipping ( βr+d ) and with clipping ( sC=βr+min(d,C) ; β=20 , C=0.45 ). Rows: KL divergence to the tilted target and relative centered-reward gap. Output laws are computed exactly by numerical integration; the dotted line marks numerical precision ( 10−13 ). ExpBoN converges geometrically in every configuration while SBoN decays polynomially; comparing the last two columns, clipping speeds up convergence by capping the right tail of the scores.