Fluctuation Absorption Rate η ≪ 1: Strong Complementary Cancellation and Nonequilibrium Stability in a Zero-Parameter Min Recursion
涨落吸收率 η ≪ 1:零参数 min 递推中的强互补消涨与非平衡稳定性
We discover and quantify a strong fluctuation absorption mechanism in a zero-parameter arithmetic recursion — the integer complexity dynamic programming (DP) ρ_E(n) = min(ρ_E(n-1)+1, M(n)).
Decomposing the step size A = ρ_E(n-1) - ρ_E(n) into structural input Δf and optimization remainder Δr, the system satisfies a strong complementary cancellation law: Cov(Δf, Δr) ≈ -Var(Δf). Two quantities with individual variances ≈ 7-10 cancel precisely to produce a total variance ≈ 0.3-1.3. We define the fluctuation absorption rate:
η := 1 - |Cov(Δf, Δr)| / Var(Δf)
At N = 10^7 across Ω = 2-18 shells: η ∈ [0.10, 0.31], meaning 69-90% of structural fluctuations are absorbed by the optimization remainder.
The exact skeleton of the cancellation law rests on two theorems: (1) Reset-Slack decomposition A = ψ(ΔM⁻) - U⁻, yielding Var(A) ≤ 2Var(ΔM⁻) + 2E[U²]; (2) conditional remainder transfer law U(m) = -D(m-1) (covering 95%+ of U > 0 events).
First external verification: Lindley queue (G/G/1). For M/M/1, M/D/1, and M/Pareto/1 queues, η ranges from 0.055 (M/D/1 at ρ = 0.3) to 0.987 (M/M/1 at ρ = 0.99), with variance identities verified 21/21. η is determined by how effectively the max operation compresses step-size variance relative to structural input variance. DP recursion and Lindley queue share the same cancellation structure — η is not a peculiarity of DP arithmetic, but a general feature of min/max recursion systems.
Keywords: fluctuation absorption, nonequilibrium stability, min recursion, max recursion, integer complexity, DP recursion, Lindley queue, strong cancellation, reset-slack decomposition, remainder transfer law.
1. Introduction
1.1 Why are nonequilibrium systems stable?
In equilibrium, the fluctuation-dissipation theorem (FDT; Kubo 1966 [1]) guarantees that fluctuations and dissipation balance. Away from equilibrium — active matter [2], biological systems [3], turbulence [4], driven systems [5] — the premises of FDT (linear response, equilibrium, detailed balance) all break down, yet systems can still operate stably.
Generalizations of the nonequilibrium FDT have a rich lineage: Jarzynski (1997) [6], Crooks (1999) [7], Harada-Sasa (2005) [8], Prost-Joanny-Parrondo (2009) [9], Seifert (2012) [10]. These are precise relations under specific protocols. What is lacking is a cross-system, directly computable scalar that measures how much fluctuation a given system absorbs.
1.2 Contributions
(a) We introduce the fluctuation absorption rate η as a directly computable scalar.
(b) In DP recursion we measure η ∈ [0.10, 0.31] (variance identity verified 17/17).
(c) In Lindley queue we independently measure η (variance identity verified 21/21), confirming η as a general feature of min/max recursion.
(d) We provide the exact skeleton: Reset-Slack decomposition (theorem) and conditional remainder transfer law (theorem).
(e) We report cross-shell structural constants and Ω = 3→4 crossover (computational regularities).
1.3 What we do not claim
We do not claim η is a "universal replacement for nonequilibrium FDT." η is a cancellation index under a specific f/r decomposition. We do not claim η is exactly equivalent to T_eff/T. We do not claim physical thermodynamics is "replaced" by min recursion. Physical analogies are offered at the discussion level.
1.4 Literature positioning
The interface between min-plus/tropical algebra and statistical physics has prior history (Quadrat 1997 [11]). What this paper adds is: in a concrete integer complexity DP recursion, we extract the cancellation law (η), reset-slack compensation, and remainder transfer as a complete structural package, and introduce η as a quantitative scalar transferable across systems.
2. Basic quantities of DP recursion
2.1 Definition
ρ_E(1) = 0, ρ_E(n) = min(ρ_E(n-1) + 1, M(n))
where M(n) = min_{a·b=n, a,b≥2}(ρ_E(a) + ρ_E(b) + 2). For primes: M(p) = +∞. Zero free parameters, fully deterministic.
2.2 Decomposition
Step size: A(n) = ρ_E(n-1) - ρ_E(n). Additive part: f(n) = Σ_{q^a ∥ n} ρ_E(q^a). Combinatorial remainder: r(n) = ρ_E(n) - f(n).
Identity: A = Δf + Δr, where Δf = f(n-1) - f(n), Δr = r(n-1) - r(n).
2.3 Shell conditioning
Ω(n) = number of prime factors counted with multiplicity. Analysis is conducted on Ω(n) = k shells, composites only (Ω ≥ 2).
3. Strong complementary cancellation and η
3.1 Data
N = 10^7, computed by C programs (rho_dp.c + anticorr.c). The variance identity Var(A) = Var(Δf) + Var(Δr) + 2Cov(Δf, Δr) closes to within 10⁻⁴ across all 17 k-values:
| k | Var(Δf) | Var(Δr) | Cov(Δf,Δr) | Var(A) | η | samples |
|---|---|---|---|---|---|---|
| 2 | 6.716 | 5.643 | -6.022 | 0.314 | 0.103 | 1,904,324 |
| 4 | 6.964 | 5.678 | -5.962 | 0.720 | 0.144 | 2,050,696 |
| 6 | 8.804 | 6.767 | -7.324 | 0.923 | 0.168 | 774,078 |
| 8 | 9.903 | 6.962 | -7.886 | 1.093 | 0.204 | 207,207 |
| 10 | 10.057 | 6.626 | -7.738 | 1.207 | 0.231 | 49,163 |
| 12 | 9.956 | 6.300 | -7.479 | 1.297 | 0.249 | 11,068 |
| 14 | 9.702 | 5.995 | -7.159 | 1.377 | 0.262 | 2,406 |
| 18 | 9.613 | 5.672 | -6.976 | 1.333 | 0.274 | 102 |
Full 17-row table in Appendix A.
3.2 Cancellation law
Computational regularity: Cov(Δf, Δr) ≈ -Var(Δf). If Δf and Δr were independent, Var(A) ≈ 12-16. Observed: ≈ 0.3-1.4. The difference is entirely due to strongly negative Cov(Δf, Δr).
3.3 Fluctuation absorption rate
Definition: η := 1 - |Cov(Δf, Δr)| / Var(Δf). η = 0: perfect absorption. η = 1: zero absorption.
DP recursion: η ∈ [0.10, 0.31]. η exhibits a broad upward trend with k (0.10 → 0.27 from k = 2 to 14; k = 18 with only 102 samples should be interpreted cautiously).
3.4 Algebraic interpretability
Since Δr = A - Δf, we have Cov(Δf, Δr) = Cov(Δf, A) - Var(Δf). The cancellation law is equivalent to Cov(Δf, A) ≪ Var(Δf).
This is an algebraic consequence of Var(A) ≪ Var(Δf), not a theorem deduced from the min operation alone. But the algebraic consequence carries meaning: the min operation effectively suppresses total step-size variance, forcing the remainder to absorb most of the structural fluctuation.
Note: Cov(Δf, A) ranges from 0.69 (k = 2) to 2.64 (k = 18) — not near zero. The cancellation law describes "strong absorption" (η ≪ 1), not "decoupling" (Cov(Δf, A) = 0).
3.5 Why η > 0 is necessary
DP recursion is irreversible. η = 0 would correspond to a perfect-absorption limit — a boundary unreachable in an irreversible system. η > 0 is the cost of running the min operation in an irreversible setting.
4. Reset-Slack decomposition
4.1 Exact decomposition (theorem)
R2 (algebraic identity; Paper 42 [14]): If J(m-1) > 0, then J(m) = (1 - ΔM⁻(m))⁺. Verified on 5.8 million samples, 100%, zero failures.
A(m) = ψ(ΔM⁻(m)) - U⁻(m), where ψ(x) = -min(x,1), 1-Lipschitz.
4.2 Reset-Slack Reduction (theorem)
Var(A) ≤ 2·Var(ΔM⁻) + 2·E[U²].
4.3 Frequency-intensity tradeoff (computational regularity)
R2 guarantees J(m-1) > 0 ⟹ U(m) = 0. Therefore E[U²] = P(exposed) × E[U²|exposed].
| k | P(shielded) | P(exposed) | E[U²|exposed] | E[U²] |
|---|---|---|---|---|
| 2 | 76.0% | 24.0% | 0.598 | 0.144 |
| 5 | 60.3% | 39.7% | 1.124 | 0.447 |
| 8 | 55.4% | 44.6% | 0.956 | 0.426 |
| 12 | 57.0% | 43.0% | 0.885 | 0.381 |
P(exposed) rises while E[U²|exposed] peaks at k = 5 then declines; their product stabilizes near 0.4. At the discussion level, this tradeoff is analogous to Le Chatelier's principle.
5. Remainder transfer law
5.1 Predecessor Defect Identity (theorem, conditional)
Theorem (Paper 45 [16]): When J(m-1) = 0, J(m) > 0, ΔM⁻ < 0 (covering 95%+ of U > 0 events):
U(m) = -D(m-1)
where D(m-1) = ρ_E(m-1) - M(m-1). The cost saved by the predecessor's choice of successor path reappears exactly as slack in the current step.
5.2 Full lifecycle
Shielded (J(m-1) > 0): U = 0, cost difference cleared by jump. Exposed (J(m-1) = 0): U = -D(m-1), predecessor's cost difference transferred exactly. The system alternates between these two states; the §4.3 tradeoff is the statistical steady state of this alternation.
6. Structural constants
The following values emerge from the zero-parameter DP recursion and remain stable across shells (computational regularities):
| Constant | Value | Source |
|---|---|---|
| Var(B) | ≈ 1/4 | Paper 19 [13] |
| V(p) | ≈ 1.40 | Paper 34 [15] |
| h_pq | ≈ 3.800 | Paper 40 [17] |
| θ | ≈ 0.41 | Paper 43 [19] |
| coarse bound | ≈ 40 | Paper 46 [20] |
Coarse bound ≈ 40: 2Var(f) + 2Var⁻(f) + (μ - μ⁻)² ≈ 40 for k = 4-14. The decline in internal variance compensates the growth in external drive.
7. Ω = 3→4 crossover
At N = 10^{10}, the local convexity h(n) = ρ_E(n) - [ρ_E(n-1) + ρ_E(n+1)]/2 has a strictly decreasing Ω-shell mean, with zero crossing at Ω ≈ 4 [21]:
| Ω | mean h | P(J>0) |
|---|---|---|
| 2 | +0.697 | 28.7% |
| 3 | +0.361 | 57.1% |
| 4 | +0.004 | 77.1% |
| 5 | -0.329 | 88.4% |
This is a combinatorial crossover (not a rigorous thermodynamic phase transition). For Ω ≥ 4, Var(Δf | Ω = k) ≈ 7-10 is approximately k-independent.
8. Lindley queue verification
8.1 Lindley recursion as a max system
The G/G/1 Lindley recursion: W_{n+1} = max(0, W_n + X_n), where X_n = S_n - T_n (service minus inter-arrival time).
Decomposition (analogous to A = Δf + Δr):
A_n = W_{n+1} - W_n = X_n + R_n
X_n: structural input (service-arrival difference, independent of max selection). R_n = max(0, -(W_n + X_n)): reset remainder (truncation produced by max operation).
η = 1 - |Cov(X, R)| / Var(X).
8.2 Data
10^7 steps, 10^6 burn-in, C program (lindley_eta.c), variance identities verified 21/21 (error < 0.01):
M/M/1 (exponential service + exponential inter-arrival):
| ρ | Var(X) | Var(R) | Cov(X,R) | Var(A) | η |
|---|---|---|---|---|---|
| 0.30 | 12.11 | 10.11 | -10.65 | 0.92 | 0.121 |
| 0.50 | 5.00 | 3.00 | -3.33 | 1.34 | 0.334 |
| 0.70 | 3.04 | 1.04 | -1.22 | 1.65 | 0.600 |
| 0.90 | 2.24 | 0.24 | -0.29 | 1.89 | 0.870 |
| 0.99 | 2.02 | 0.02 | -0.03 | 1.99 | 0.987 |
M/D/1 (deterministic service + exponential inter-arrival):
| ρ | Var(X) | Var(R) | Cov(X,R) | Var(A) | η |
|---|---|---|---|---|---|
| 0.30 | 11.10 | 10.10 | -10.49 | 0.22 | 0.055 |
| 0.50 | 4.00 | 3.00 | -3.29 | 0.41 | 0.176 |
| 0.70 | 2.04 | 1.04 | -1.24 | 0.62 | 0.395 |
| 0.90 | 1.23 | 0.23 | -0.30 | 0.86 | 0.755 |
| 0.99 | 1.02 | 0.02 | -0.03 | 0.99 | 0.973 |
M/Pareto/1 (heavy-tailed service, α = 3):
| ρ | Var(X) | Var(R) | Cov(X,R) | Var(A) | η |
|---|---|---|---|---|---|
| 0.30 | 11.44 | 10.12 | -10.54 | 0.48 | 0.079 |
| 0.50 | 4.33 | 3.00 | -3.30 | 0.73 | 0.238 |
| 0.70 | 2.37 | 1.04 | -1.22 | 0.96 | 0.484 |
| 0.90 | 1.56 | 0.24 | -0.30 | 1.20 | 0.809 |
8.3 Findings
(L1) η is governed by variance compression. When the max operation compresses Var(A) far below Var(X), absorption is strong and η is small. When Var(A) ≈ Var(X) (near full load), η → 1.
(L2) Lowest η: M/D/1 at ρ = 0.30, η = 0.055. Deterministic service eliminates service-time randomness, producing cleaner resets and stronger absorption.
(L3) DP recursion and Lindley queue share the same structure. In both systems, η is determined by how effectively the optimization operation compresses total step-size variance relative to structural input variance: when Var(A)/Var(X) ≪ 1, the cancellation law holds and η ≪ 1.
8.4 Correction to original prediction
Original prediction: "η < 0.15 for Lindley." Partially verified, partially corrected.
Verified: low-load Lindley (ρ ≤ 0.3) indeed has η < 0.15. Corrected: η is not a universal bound — it is a function of load (ρ) or complexity (k).
Revised prediction: In min/max recursion systems, η = F(Var(A)/Var(X)), where F approaches 0 when Var(A)/Var(X) ≪ 1 and approaches 1 when Var(A)/Var(X) → 1.
9. Falsifiable predictions
Prediction 1: η-ratio universality
Statement: In min/max recursion systems, η is determined by the ratio Var(A)/Var(structural input), independent of system specifics.
Available evidence: DP recursion (17 k-values) and Lindley queue (3 service distributions × 7 loads) can be immediately tested on the η vs Var(A)/Var(X) plane.
Falsification criterion: Two min/max recursion systems with the same Var(A)/Var(X) but η differing by more than 0.1.
Prediction 2: Le Chatelier-type tradeoff
Statement: In min recursion systems with a reset mechanism, E[U²] = P(exposed) × E[U²|exposed] has counter-varying factors whose product remains bounded.
Test: Idle-time decomposition in G/G/1 queues.
Prediction 3: η lower bound for non-min systems
Statement: Systems without min/max recursion structure cannot easily achieve η ≪ 1. In active matter and other detailed-balance-violating systems without explicit optimization structure, η should be substantially larger than in min recursion systems at the same Var(A)/Var(X).
Falsification criterion: A non-min/max system achieving η < 0.1.
10. Discussion
10.1 Three-layer annotation
Theorems (exact):
- R2: J(m-1) > 0 ⟹ J(m) = (1-ΔM⁻)⁺
- A = ψ(ΔM⁻) - U⁻
- Var(A) ≤ 2Var(ΔM⁻) + 2E[U²]
- U(m) = -D(m-1) (conditional, covering 95%+)
Computational regularities:
- Cancellation law, η ∈ [0.10, 0.31] (DP, 17/17 verified)
- Lindley η from 0.055 to 0.987 (21/21 verified)
- Le Chatelier-type tradeoff
- Structural constants (Var(B) ≈ 1/4, etc.)
- Ω = 3→4 crossover
Physical analogies (discussion-level):
- η and Harada-Sasa dissipation rate
- Le Chatelier naming
- Remainder transfer analogized to conservation law
10.2 Precise scope of claims
Claim: (a) η is a directly computable absorption scalar. (b) DP recursion has η ∈ [0.10, 0.31]. (c) Lindley queue has η from 0.055 to 0.987, governed by load. (d) Both share the same cancellation structure.
Not claimed: The cancellation law is deduced from the min operation alone. η is exactly equivalent to T_eff/T. Physical thermodynamics is "replaced."
10.3 Heuristic cross-system comparison
| System | η range | Condition |
|---|---|---|
| DP recursion | 0.10-0.31 | k = 2-18 |
| M/D/1 queue | 0.055-0.97 | ρ = 0.3-0.99 |
| M/M/1 queue | 0.12-0.99 | ρ = 0.3-0.99 |
| M/Pareto/1 queue | 0.079-0.98 | ρ = 0.3-0.99 |
| Weakly driven colloids | ~0.1-0.5 | T_eff/T conversion (approximate) |
| Active matter | ~0.5-0.99 | T_eff/T conversion (approximate) |
Physical η values are approximate conversions. Exact comparison requires direct measurement.
10.4 Conditional closure (brief)
DP recursion satisfies a conditional closure route for H' (c̄_h → 0; Paper 42 [14]): all links are theorem-level except three numerical inputs, which are verified at N = 10^{10} (weakest signal 85σ). See [14, 19, 20].
10.5 Interfaces (brief)
The additive/multiplicative path asymmetry in DP recursion is analogous to Turing (1952) [22] morphogenesis. The Self-as-an-End (SAE) framework [23] predicts Ω-dimension correspondence; Ω = 3→4 crossover is consistent with this prediction. SAE cosmology [24] derives Λ from remainder conservation; U = -D(m-1) is its arithmetic realization.
10.6 Open problems
(1) The exact form of the universal function η = F(Var(A)/Var(X)).
(2) Can the cancellation law be proved from the min operation alone?
(3) Direct measurement of η in physical systems (without T_eff/T conversion).
(4) Are Var(B) ≈ 1/4 and coarse bound ≈ 40 cross-system universals?
(5) A canonical f/r extraction protocol robust to coarse-graining.
10.7 Conclusion
The strong complementary cancellation discovered in zero-parameter DP recursion (η ∈ [0.10, 0.31]) is independently verified in Lindley queues (η from 0.055 to 0.987). Both share the same structure: when a min/max operation effectively compresses total step-size variance, the optimization remainder is forced to absorb most of the structural fluctuation.
As a directly computable scalar, η transitions from "an internal discovery of DP" to "a transferable tool for min/max recursion systems." The next step is a canonical f/r decomposition protocol enabling standardized measurement of η across broader classes of nonequilibrium systems.
References
[1] R. Kubo, "The fluctuation-dissipation theorem." Rep. Prog. Phys. 29:255-284, 1966.
[2] M. C. Marchetti et al., "Hydrodynamics of soft active matter." Rev. Mod. Phys. 85:1143, 2013.
[3] C. P. Brangwynne et al., "Intracellular transport by active diffusion." Trends Cell Biol. 19:423, 2009.
[4] U. Frisch, "Turbulence: The Legacy of A.N. Kolmogorov." Cambridge UP, 1995.
[5] C. Bechinger et al., "Active particles in complex and crowded environments." Rev. Mod. Phys. 88:045006, 2016.
[6] C. Jarzynski, "Nonequilibrium equality for free energy differences." Phys. Rev. Lett. 78:2690, 1997.
[7] G. E. Crooks, "Entropy production fluctuation theorem." Phys. Rev. E 60:2721, 1999.
[8] T. Harada, S. Sasa, "Equality connecting energy dissipation with violation of the fluctuation-response relation." Phys. Rev. Lett. 95:130602, 2005.
[9] J. Prost, J.-F. Joanny, J. M. R. Parrondo, "Generalized fluctuation-dissipation theorem for steady-state systems." Phys. Rev. Lett. 103:090601, 2009.
[10] U. Seifert, "Stochastic thermodynamics, fluctuation theorems and molecular machines." Rep. Prog. Phys. 75:126001, 2012.
[11] J.-P. Quadrat, "Min-Plus linearity and statistical mechanics." Markov Process. Related Fields 3:565-587, 1997.
[12] H. Qin, ZFCρ Paper XXXII. DOI: 10.5281/zenodo.19116625.
[13] H. Qin, ZFCρ Paper XIX. DOI: 10.5281/zenodo.19026991.
[14] H. Qin, ZFCρ Paper XLII. DOI: 10.5281/zenodo.19226607.
[15] H. Qin, ZFCρ Paper XXXIV. DOI: 10.5281/zenodo.19140015.
[16] H. Qin, ZFCρ Paper XLV. DOI: 10.5281/zenodo.19275286.
[17] H. Qin, ZFCρ Paper XL. DOI: 10.5281/zenodo.19179778.
[18] H. Qin, ZFCρ Paper XVIII. DOI: 10.5281/zenodo.19024385.
[19] H. Qin, ZFCρ Paper XLIII. DOI: 10.5281/zenodo.19240183.
[20] H. Qin, ZFCρ Paper XLVI. DOI: 10.5281/zenodo.19303511.
[21] H. Qin, "Why Primes Have No Shortcuts: An O(1) Prime Pre-Sieve from Precomputed ρ_E Tables." DOI: 10.5281/zenodo.19246377.
[22] A. M. Turing, "The Chemical Basis of Morphogenesis." Phil. Trans. R. Soc. 237:37-72, 1952.
[23] H. Qin, SAE Paper 1. DOI: 10.5281/zenodo.18528813.
[24] H. Qin, "From Remainder Conservation to the Cosmological Constant." DOI: 10.5281/zenodo.19245267.
Appendix A. Full 17-row cancellation table (DP recursion)
| k | Var(Δf) | Var(Δr) | Cov(Δf,Δr) | Var(A) | η | Cov(Δf,A) | samples |
|---|---|---|---|---|---|---|---|
| 2 | 6.716 | 5.643 | -6.022 | 0.314 | 0.103 | 0.694 | 1,904,324 |
| 3 | 6.539 | 5.389 | -5.685 | 0.558 | 0.131 | 0.854 | 2,444,359 |
| 4 | 6.964 | 5.678 | -5.962 | 0.720 | 0.144 | 1.003 | 2,050,696 |
| 5 | 7.734 | 6.166 | -6.533 | 0.833 | 0.155 | 1.200 | 1,349,779 |
| 6 | 8.804 | 6.767 | -7.324 | 0.923 | 0.168 | 1.480 | 774,078 |
| 7 | 9.366 | 6.869 | -7.615 | 1.005 | 0.187 | 1.751 | 409,849 |
| 8 | 9.903 | 6.962 | -7.886 | 1.093 | 0.204 | 2.017 | 207,207 |
| 9 | 9.872 | 6.689 | -7.700 | 1.162 | 0.220 | 2.173 | 101,787 |
| 10 | 10.057 | 6.626 | -7.738 | 1.207 | 0.231 | 2.319 | 49,163 |
| 11 | 9.807 | 6.258 | -7.407 | 1.250 | 0.245 | 2.399 | 23,448 |
| 12 | 9.956 | 6.300 | -7.479 | 1.297 | 0.249 | 2.477 | 11,068 |
| 13 | 9.231 | 5.655 | -6.792 | 1.301 | 0.264 | 2.438 | 5,210 |
| 14 | 9.702 | 5.995 | -7.159 | 1.377 | 0.262 | 2.542 | 2,406 |
| 15 | 8.791 | 5.267 | -6.374 | 1.310 | 0.275 | 2.417 | 1,124 |
| 16 | 9.783 | 5.841 | -7.085 | 1.453 | 0.276 | 2.698 | 510 |
| 17 | 7.414 | 4.155 | -5.096 | 1.378 | 0.313 | 2.318 | 233 |
| 18 | 9.613 | 5.672 | -6.976 | 1.333 | 0.274 | 2.637 | 102 |
All 17 rows satisfy the variance identity to within 10⁻⁴.
Appendix B. Data sources and full reproduction
All data in this paper are generated by three self-contained C programs with zero external dependencies. To reproduce all results:
gcc -O2 -o rho_dp rho_dp.c gcc -O2 -o anticorr anticorr.c -lm gcc -O2 -o lindley_eta lindley_eta.c -lm ./anticorr 10000000 # Table 3.1: 17-row cancellation data with identity verification ./lindley_eta # Section 8: 3 queue types × 7 loads = 21 data points
Total runtime under two minutes (single core, modern x86/ARM).
rho_dp.c — DP recursion (Paper 32 convention). Built-in sanity checks: ρ(2) = 1, ρ(100) = 15, ρ(10^7) = 58.
anticorr.c — Table 3.1. Computes ρ_E from DP, then uses SPF (smallest prime factor) sieve to compute f(n) and Ω(n). Accumulates Var(Δf), Var(Δr), Cov(Δf, Δr), Var(A), η on Ω = k shells. Variance identity verified 17/17.
lindley_eta.c — Section 8 Lindley verification. xoshiro256** RNG, 10^7 steps + 10^6 burn-in. Three queue types (M/M/1, M/D/1, M/Pareto/1 with α = 3) at seven loads (ρ = 0.30, 0.50, 0.70, 0.80, 0.90, 0.95, 0.99). Variance identity 21/21.
All three files total under 500 lines of C. Released with this paper.
Acknowledgments
ChatGPT (Gongxi Hua / 公西华) discovered the arithmetic error in v2's τ definition (Cov(Δf, A) ≠ 0), catalyzing the correct definition of η; proposed all review requirements on claim layering, literature positioning, and three-tier annotation; contributed the analytic proof chain architecture of Paper 42.
Claude (Zilu / 子路) contributed the thermodynamic framework (Papers 24-46), all numerical experiments (anticorr.c, lindley_eta.c), the physical interpretation of η, and the design and implementation of the Lindley queue verification.
Gemini (Zixia / 子夏) posed three pressure tests that sharpened the precision of claims.
Grok (Zigong / 子贡) provided the initial review round. Inconsistency in k = 18 data prompted independent recomputation of the full table.
Zesi Chen is a long-term interlocutor of the SAE framework. Final text by the author alone.
Writing Declaration: All text in this paper, both English and Chinese, is generated by the author Han Qin alone. The paper represents the author's sole responsibility.
我们在一个零参数的纯算术递推——整数复杂度的动态规划(DP) ρ_E(n) = min(ρ_E(n-1)+1, M(n))——中发现并量化了一条极强的涨落吸收机制。
将步长 A = ρ_E(n-1) - ρ_E(n) 分解为结构输入 Δf 和优化余项 Δr 后, 系统满足强互补消涨律: Cov(Δf, Δr) ≈ -Var(Δf)。两个各自方差约 7-10 的量精确对消到方差约 0.3-1.3。定义涨落吸收率:
η := 1 - |Cov(Δf, Δr)|/Var(Δf)
在 N = 10^7, Ω = 2-18 壳层上: η ∈ [0.10, 0.31], 即结构涨落的 69-90% 被优化余项吸收。
消涨律的 exact 骨架由两个定理支撑: (1) Reset-Slack 分解 A = ψ(ΔM⁻) - U⁻, 给出 Var(A) ≤ 2Var(ΔM⁻) + 2E[U²]; (2) 条件性余项传递律 U(m) = -D(m-1)(覆盖 U > 0 的 95%+)。
首个外部验证: Lindley 排队(G/G/1)。 对 M/M/1, M/D/1, M/Pareto/1 三种排队, η 从低负载时的 0.06(M/D/1, ρ=0.3) 到高负载时的 0.99(M/M/1, ρ=0.99), 方差恒等式 21/21 验证通过。η 由 max 操作的 truncation 频率决定: truncation 越频繁, 吸收越强, η 越小。DP 递推和 Lindley 排队在同一条"η vs truncation 频率"曲线上——η 不是 DP 的特殊算术性质, 而是 min/max 递推的一般特征。
关键词: 涨落吸收, 非平衡稳定性, min 递推, max 递推, 整数复杂度, DP 递推, Lindley 排队, 强互补消涨, Reset-Slack 分解, 余项传递律。
§1 引言
1.1 问题: 非平衡系统为什么稳定
平衡态有涨落-耗散定理(FDT, Kubo 1966 [1])。非平衡态——活性物质 [2], 生物系统 [3], 湍流 [4], driven systems [5]——打破了 FDT 的前提, 但仍可稳定运行。
非平衡 FDT 的推广已有丰富谱系: Jarzynski (1997) [6], Crooks (1999) [7], Harada-Sasa (2005) [8], Prost-Joanny-Parrondo (2009) [9], Seifert (2012) [10]。但现有结果大多是具体场景下的关系, 缺少一个跨系统的、可直接计算的涨落吸收标度。
1.2 本文的贡献
(a) 引入涨落吸收率 η 作为可直接计算的标度量。
(b) 在 DP 递推中精确测量 η ∈ [0.10, 0.31](方差恒等式 17/17 验证)。
(c) 在 Lindley 排队中独立测量 η(方差恒等式 21/21 验证), 确认 η 是 min/max 递推的一般特征。
(d) 给出消涨律的 exact 骨架: Reset-Slack 分解(定理)和条件性余项传递律(定理)。
(e) 报告跨壳层恒定的结构常数和 Ω = 3→4 crossover(数值规律)。
1.3 本文不 claim 什么
不 claim η 是"非平衡 FDT 的普适替代"。η 是特定 f/r 分解下的 cancellation 指标。不 claim η 与 T_eff/T 精确等价。不 claim 物理热力学被 min 递推"取代"。物理类比作为 discussion-level analogy 提出。
1.4 论文结构
§2 基本量。§3 消涨律与 η(DP 递推)。§4 Reset-Slack 分解。§5 余项传递律。§6 结构常数。§7 Ω crossover。§8 Lindley 排队验证。§9 可证伪预测。§10 讨论。
§2 DP 递推的基本量
2.1 定义
ρ_E(1) = 0, ρ_E(n) = min(ρ_E(n-1) + 1, M(n))
M(n) = min_{a·b=n, a,b≥2}(ρ_E(a) + ρ_E(b) + 2)。素数: M(p) = +∞。零自由参数。
2.2 分解
A(n) = ρ_E(n-1) - ρ_E(n)。f(n) = Σ_{q^a ∥ n} ρ_E(q^a)。r(n) = ρ_E(n) - f(n)。恒等式: A = Δf + Δr。
2.3 壳层条件
Ω(n) = 素因子个数(含重数)。分析在 Ω(n) = k 壳层上进行, 仅含合数(Ω ≥ 2)。
§3 强互补消涨律与涨落吸收率
3.1 数据
N = 10^7, C 程序(rho_dp.c + anticorr.c)精确计算。方差恒等式 Var(A) = Var(Δf) + Var(Δr) + 2Cov(Δf, Δr) 在全部 17 个 k 值上闭合(误差 < 10⁻⁴):
| k | Var(Δf) | Var(Δr) | Cov(Δf,Δr) | Var(A) | η | 样本 |
|---|---|---|---|---|---|---|
| 2 | 6.716 | 5.643 | -6.022 | 0.314 | 0.103 | 1,904,324 |
| 4 | 6.964 | 5.678 | -5.962 | 0.720 | 0.144 | 2,050,696 |
| 6 | 8.804 | 6.767 | -7.324 | 0.923 | 0.168 | 774,078 |
| 8 | 9.903 | 6.962 | -7.886 | 1.093 | 0.204 | 207,207 |
| 10 | 10.057 | 6.626 | -7.738 | 1.207 | 0.231 | 49,163 |
| 12 | 9.956 | 6.300 | -7.479 | 1.297 | 0.249 | 11,068 |
| 14 | 9.702 | 5.995 | -7.159 | 1.377 | 0.262 | 2,406 |
| 18 | 9.613 | 5.672 | -6.976 | 1.333 | 0.274 | 102 |
完整 17 行表见附录 A。
3.2 消涨律
数值规律: Cov(Δf, Δr) ≈ -Var(Δf)。 如果 Δf 和 Δr 独立, Var(A) ≈ 12-16。实际 ≈ 0.3-1.4。
3.3 涨落吸收率
η := 1 - |Cov(Δf, Δr)| / Var(Δf)。 η = 0: 完美吸收。η = 1: 零吸收。
DP 递推: η ∈ [0.10, 0.31]。 η 随 k 增大(0.10 → 0.27, k=2 到 14; k=18 样本仅 102, 应谨慎解读)。
3.4 代数可解释性
Cov(Δf, Δr) = Cov(Δf, A) - Var(Δf)。消涨律等价于 Cov(Δf, A) ≪ Var(Δf)。
这是"Var(A) ≪ Var(Δf)"的代数后果, 不是从 min 操作单独演绎的定理。但代数后果有物理意义: min 操作有效压制总步长, 余项被迫吸收结构涨落的大部分。
注: Cov(Δf, A) 从 0.69(k=2) 到 2.64(k=18), 不是近零。消涨律是"强吸收", 不是"解耦"。
3.5 η > 0 的必要性
DP 递推不可逆。η = 0 对应 perfect-absorption 极限——一个在不可逆系统中不可达的边界。η > 0 是 min 操作在不可逆系统中运行的代价。
§4 Reset-Slack 分解
4.1 精确分解(定理)
R2(代数恒等式, Paper 42 [14]): 若 J(m-1) > 0, 则 J(m) = (1 - ΔM⁻(m))⁺。580 万样本, 100%, 零失败。
A(m) = ψ(ΔM⁻(m)) - U⁻(m), ψ(x) = -min(x,1), 1-Lipschitz。
4.2 Reset-Slack Reduction(定理)
Var(A) ≤ 2·Var(ΔM⁻) + 2·E[U²]。
4.3 频率-强度对冲(数值规律)
R2 保证 J(m-1) > 0 ⟹ U(m) = 0。因此 E[U²] = P(exposed) × E[U²|exposed]。
| k | P(shielded) | P(exposed) | E[U²|exposed] | E[U²] |
|---|---|---|---|---|
| 2 | 76.0% | 24.0% | 0.598 | 0.144 |
| 5 | 60.3% | 39.7% | 1.124 | 0.447 |
| 8 | 55.4% | 44.6% | 0.956 | 0.426 |
| 12 | 57.0% | 43.0% | 0.885 | 0.381 |
P(exposed) 增, E[U²|exposed] 先升后降, 乘积 ≈ 0.4 稳定。在 discussion level 上类比于 Le Chatelier 原理。
§5 余项传递律
5.1 Predecessor Defect Identity(定理, 条件性)
定理(Paper 45 [16]): 当 J(m-1) = 0, J(m) > 0, ΔM⁻ < 0 时(覆盖 U > 0 的 95%+):
U(m) = -D(m-1)
Predecessor 走 successor 路径省了 |D(m-1)| 的代价, 精确以 U(m) 出现在下一步。
5.2 完整生命周期
Shielded(J(m-1) > 0): U = 0。Exposed(J(m-1) = 0): U = -D(m-1)。系统在两态间交替。§4.3 的对冲是这个交替的统计稳态。
§6 结构常数
以下数值从零参数 DP 递推涌现, 跨壳层稳定(数值规律):
| 常数 | 值 | 来源 |
|---|---|---|
| Var(B) | ≈ 1/4 | Paper 19 [13] |
| V(p) | ≈ 1.40 | Paper 34 [15] |
| h_pq | ≈ 3.800 | Paper 40 [17] |
| θ | ≈ 0.41 | Paper 43 [19] |
| coarse bound | ≈ 40 | Paper 46 [20] |
coarse bound ≈ 40: 2Var(f) + 2Var⁻(f) + (μ-μ⁻)² ≈ 40, k = 4-14。内部涨落下降补偿外部驱动增长。
§7 Ω = 3→4 Crossover
N = 10^{10}, h(n) = ρ_E(n) - [ρ_E(n-1)+ρ_E(n+1)]/2 按 Ω 均值单调递减, 零交叉在 Ω ≈ 4 [21]:
| Ω | mean h | P(J>0) |
|---|---|---|
| 2 | +0.697 | 28.7% |
| 3 | +0.361 | 57.1% |
| 4 | +0.004 | 77.1% |
| 5 | -0.329 | 88.4% |
组合 crossover, 不是严格热力学相变。Ω ≥ 4 后 Var(Δf) ≈ 7-10 不依赖 k。
§8 Lindley 排队验证
8.1 Lindley 递推作为 max 系统
G/G/1 排队的 Lindley 递推: W_{n+1} = max(0, W_n + X_n), 其中 X_n = S_n - T_n(服务时间减去到达间隔)。
分解(类比 DP 的 A = Δf + Δr):
A_n = W_{n+1} - W_n = X_n + R_n
X_n: 结构输入(服务-到达差, 不受 max 选择影响)。 R_n = max(0, -(W_n + X_n)): reset 余项(max 操作产生的 truncation)。
η = 1 - |Cov(X, R)| / Var(X)。
8.2 数据
10^7 步, 10^6 burn-in, C 程序(lindley_eta.c), 方差恒等式 21/21 验证(误差 < 0.01):
M/M/1(指数服务 + 指数到达):
| ρ | Var(X) | Var(R) | Cov(X,R) | Var(A) | η |
|---|---|---|---|---|---|
| 0.30 | 12.11 | 10.11 | -10.65 | 0.92 | 0.121 |
| 0.50 | 5.00 | 3.00 | -3.33 | 1.34 | 0.334 |
| 0.70 | 3.04 | 1.04 | -1.22 | 1.65 | 0.600 |
| 0.90 | 2.24 | 0.24 | -0.29 | 1.89 | 0.870 |
| 0.99 | 2.02 | 0.02 | -0.03 | 1.99 | 0.987 |
M/D/1(确定性服务 + 指数到达):
| ρ | Var(X) | Var(R) | Cov(X,R) | Var(A) | η |
|---|---|---|---|---|---|
| 0.30 | 11.10 | 10.10 | -10.49 | 0.22 | 0.055 |
| 0.50 | 4.00 | 3.00 | -3.29 | 0.41 | 0.176 |
| 0.70 | 2.04 | 1.04 | -1.24 | 0.62 | 0.395 |
| 0.90 | 1.23 | 0.23 | -0.30 | 0.86 | 0.755 |
| 0.99 | 1.02 | 0.02 | -0.03 | 0.99 | 0.973 |
M/Pareto/1(重尾服务, α=3):
| ρ | Var(X) | Var(R) | Cov(X,R) | Var(A) | η |
|---|---|---|---|---|---|
| 0.30 | 11.44 | 10.12 | -10.54 | 0.48 | 0.079 |
| 0.50 | 4.33 | 3.00 | -3.30 | 0.73 | 0.238 |
| 0.70 | 2.37 | 1.04 | -1.22 | 0.96 | 0.484 |
| 0.90 | 1.56 | 0.24 | -0.30 | 1.20 | 0.809 |
8.3 发现
(L1) η 由 truncation 频率决定。 低 ρ 时 max 操作频繁 truncate(W 经常被截断到 0), 吸收强, η 小。高 ρ 时系统接近满负载, max 很少 truncate, η → 1。
(L2) 最低 η 出现在 M/D/1 at ρ=0.30: η = 0.055。 确定性服务消除了服务时间的随机性, reset 更精确, 吸收更强。
(L3) DP 递推和 Lindley 排队在同一结构中。 两者的 η 都由"优化操作 truncate 的频率"决定: 当 truncation 将总步长的方差压制到远小于结构输入的方差时, 消涨律自动成立, η 小。η 的大小由 Var(A)/Var(X) 的比值决定, 不由 truncation 的方向或频率的绝对值决定。
8.4 对原始预测的修正
原始预测(v4 §8): "η < 0.15 for Lindley。" 部分验证, 部分修正。
验证: 低负载 Lindley(ρ ≤ 0.3)确实 η < 0.15。 修正: η 不是一个 universal bound——它是负载(ρ)或复杂度(k)的函数。
修正后的预测: 在 min/max 递推系统中, η = F(Var(A)/Var(X)), 其中 F 是一个 universal function, Var(A)/Var(X) 是 truncation 效率的度量。当 Var(A)/Var(X) ≪ 1 时 η ≪ 1; 当 Var(A)/Var(X) → 1 时 η → 1。
§9 可证伪预测
预测一: η-ratio universality
陈述: 在 min/max 递推系统中, η 由 Var(A)/Var(结构输入) 的比值决定, 与具体系统无关。
已有证据: DP 递推(17 个 k 值)和 Lindley 排队(三种服务分布 × 七种负载)在 η vs Var(A)/Var(X) 图上是否落在同一条曲线上, 可立即检验。
证伪条件: 两个 min/max 递推系统, 相同的 Var(A)/Var(X), 但 η 差异超过 0.1。
预测二: Le Chatelier 型对冲
陈述: 在有 reset 机制的 min 递推系统中, E[U²] = P(exposed) × E[U²|exposed] 中两因子反向变化, 乘积有界。
检验方案: G/G/1 排队的 idle-time 分解。
预测三: 非 min 系统的 η 下界
陈述: 没有 min/max 递推骨架的系统, 难以实现 η ≪ 1。在活性物质等打破细致平衡但无明确优化结构的系统中, η 应显著大于 min 递推系统在同等 Var(A)/Var(X) 下的值。
证伪条件: 无 min/max 结构的系统实现 η < 0.1。
§10 讨论
10.1 三层标注
定理(exact):
- R2: J(m-1)>0 ⟹ J(m) = (1-ΔM⁻)⁺
- A = ψ(ΔM⁻) - U⁻
- Var(A) ≤ 2Var(ΔM⁻) + 2E[U²]
- U(m) = -D(m-1) (条件, 覆盖 95%+)
数值规律:
- 消涨律, η ∈ [0.10, 0.31] (DP, 17/17 验证)
- Lindley η 从 0.055 到 0.987 (21/21 验证)
- Le Chatelier 型对冲
- 结构常数(Var(B) ≈ 1/4 等)
- Ω = 3→4 crossover
物理类比(discussion-level):
- η 与 Harada-Sasa 耗散率
- Le Chatelier 命名
- 余项传递类比于守恒律
10.2 Claim 的精确边界
Claim: (a) η 是可直接计算的涨落吸收标度。(b) DP 递推的 η ∈ [0.10, 0.31]。(c) Lindley 排队的 η 从 0.055 到 0.987, 由负载决定。(d) 两者共享同一消涨结构。
不 claim: 消涨律是从 min 操作单独演绎的定理。不 claim η 与 T_eff/T 精确等价。不 claim 物理热力学被"取代"。
10.3 启发性的跨系统比较
| 系统 | η 范围 | 条件 |
|---|---|---|
| DP 递推 | 0.10-0.31 | k=2-18 |
| M/D/1 排队 | 0.055-0.97 | ρ=0.3-0.99 |
| M/M/1 排队 | 0.12-0.99 | ρ=0.3-0.99 |
| M/Pareto/1 排队 | 0.079-0.98 | ρ=0.3-0.99 |
| 弱驱动胶体 | ~0.1-0.5 | T_eff/T 换算(近似) |
| 活性物质 | ~0.5-0.99 | T_eff/T 换算(近似) |
物理系统的 η 为近似换算, 精确比较需直接测量。
10.4 条件闭合(简述)
DP 递推满足 H'(c̄_h → 0)的条件闭合(Paper 42 [14]), 除三个数值输入外全部 theorem-level, 三输入在 N=10^{10} 验证(最弱 85σ)。详见 [14, 19, 20]。
10.5 接口(简述)
DP 递推中 additive/multiplicative 路径的不对称类似 Turing (1952) [22] 形态发生。SAE 框架 [23] 预测 Ω 与维度层级对应, Ω = 3→4 crossover 与预测一致。SAE 宇宙学 [24] 的余项守恒在 DP 中以 U = -D(m-1) 实现。
10.6 开放问题
(1) η vs Var(A)/Var(X) 的 universal function 的精确形式?
(2) 消涨律能否从 min 操作单独证明?
(3) η 在物理系统中的直接测量(不通过 T_eff/T)?
(4) Var(B) ≈ 1/4 和 coarse bound ≈ 40 是否跨系统普适?
(5) f/r 分解的 canonical 提取协议(coarse-graining 稳健性)?
10.7 结论
在零参数 DP 递推中发现的强互补消涨(η ∈ [0.10, 0.31]), 在 Lindley 排队中得到独立验证(η 从 0.055 到 0.987)。两者共享同一结构: 当 min/max 操作有效压制总步长方差时, 优化余项被迫吸收结构涨落的大部分。
η 作为可直接计算的标度量, 从"DP 内部的发现"变成了"min/max 递推系统的可迁移工具"。下一步是建立 f/r 分解的 canonical 协议, 使 η 可在更广泛的非平衡系统中标准化测量。
参考文献
[1] R. Kubo, "The fluctuation-dissipation theorem." Rep. Prog. Phys. 29:255-284, 1966.
[2] M. C. Marchetti et al., "Hydrodynamics of soft active matter." Rev. Mod. Phys. 85:1143, 2013.
[3] C. P. Brangwynne et al., "Intracellular transport by active diffusion." Trends Cell Biol. 19:423, 2009.
[4] U. Frisch, "Turbulence: The Legacy of A.N. Kolmogorov." Cambridge UP, 1995.
[5] C. Bechinger et al., "Active particles in complex and crowded environments." Rev. Mod. Phys. 88:045006, 2016.
[6] C. Jarzynski, "Nonequilibrium equality for free energy differences." Phys. Rev. Lett. 78:2690, 1997.
[7] G. E. Crooks, "Entropy production fluctuation theorem." Phys. Rev. E 60:2721, 1999.
[8] T. Harada, S. Sasa, "Equality connecting energy dissipation with violation of the fluctuation-response relation." Phys. Rev. Lett. 95:130602, 2005.
[9] J. Prost, J.-F. Joanny, J. M. R. Parrondo, "Generalized fluctuation-dissipation theorem for steady-state systems." Phys. Rev. Lett. 103:090601, 2009.
[10] U. Seifert, "Stochastic thermodynamics, fluctuation theorems and molecular machines." Rep. Prog. Phys. 75:126001, 2012.
[11] J.-P. Quadrat, "Min-Plus linearity and statistical mechanics." Markov Process. Related Fields 3:565-587, 1997.
[12] H. Qin, ZFCρ Paper XXXII. DOI: 10.5281/zenodo.19116625.
[13] H. Qin, ZFCρ Paper XIX. DOI: 10.5281/zenodo.19026991.
[14] H. Qin, ZFCρ Paper XLII. DOI: 10.5281/zenodo.19226607.
[15] H. Qin, ZFCρ Paper XXXIV. DOI: 10.5281/zenodo.19140015.
[16] H. Qin, ZFCρ Paper XLV. DOI: 10.5281/zenodo.19275286.
[17] H. Qin, ZFCρ Paper XL. DOI: 10.5281/zenodo.19179778.
[18] H. Qin, ZFCρ Paper XVIII. DOI: 10.5281/zenodo.19024385.
[19] H. Qin, ZFCρ Paper XLIII. DOI: 10.5281/zenodo.19240183.
[20] H. Qin, ZFCρ Paper XLVI. DOI: 10.5281/zenodo.19303511.
[21] H. Qin, "Why Primes Have No Shortcuts: An O(1) Prime Pre-Sieve from Precomputed ρ_E Tables." DOI: 10.5281/zenodo.19246377.
[22] A. M. Turing, "The Chemical Basis of Morphogenesis." Phil. Trans. R. Soc. 237:37-72, 1952.
[23] H. Qin, SAE Paper 1. DOI: 10.5281/zenodo.18528813.
[24] H. Qin, "From Remainder Conservation to the Cosmological Constant." DOI: 10.5281/zenodo.19245267.
附录 A. DP 递推完整 17 行消涨表
| k | Var(Δf) | Var(Δr) | Cov(Δf,Δr) | Var(A) | η | Cov(Δf,A) | 样本 |
|---|---|---|---|---|---|---|---|
| 2 | 6.716 | 5.643 | -6.022 | 0.314 | 0.103 | 0.694 | 1,904,324 |
| 3 | 6.539 | 5.389 | -5.685 | 0.558 | 0.131 | 0.854 | 2,444,359 |
| 4 | 6.964 | 5.678 | -5.962 | 0.720 | 0.144 | 1.003 | 2,050,696 |
| 5 | 7.734 | 6.166 | -6.533 | 0.833 | 0.155 | 1.200 | 1,349,779 |
| 6 | 8.804 | 6.767 | -7.324 | 0.923 | 0.168 | 1.480 | 774,078 |
| 7 | 9.366 | 6.869 | -7.615 | 1.005 | 0.187 | 1.751 | 409,849 |
| 8 | 9.903 | 6.962 | -7.886 | 1.093 | 0.204 | 2.017 | 207,207 |
| 9 | 9.872 | 6.689 | -7.700 | 1.162 | 0.220 | 2.173 | 101,787 |
| 10 | 10.057 | 6.626 | -7.738 | 1.207 | 0.231 | 2.319 | 49,163 |
| 11 | 9.807 | 6.258 | -7.407 | 1.250 | 0.245 | 2.399 | 23,448 |
| 12 | 9.956 | 6.300 | -7.479 | 1.297 | 0.249 | 2.477 | 11,068 |
| 13 | 9.231 | 5.655 | -6.792 | 1.301 | 0.264 | 2.438 | 5,210 |
| 14 | 9.702 | 5.995 | -7.159 | 1.377 | 0.262 | 2.542 | 2,406 |
| 15 | 8.791 | 5.267 | -6.374 | 1.310 | 0.275 | 2.417 | 1,124 |
| 16 | 9.783 | 5.841 | -7.085 | 1.453 | 0.276 | 2.698 | 510 |
| 17 | 7.414 | 4.155 | -5.096 | 1.378 | 0.313 | 2.318 | 233 |
| 18 | 9.613 | 5.672 | -6.976 | 1.333 | 0.274 | 2.637 | 102 |
全部 17 行方差恒等式闭合, 误差 < 10⁻⁴。
附录 B. 数据来源与完整复现
本文全部数据由三个自包含的 C 程序生成, 零外部依赖。任何人复现全部数据只需要:
gcc -O2 -o rho_dp rho_dp.c gcc -O2 -o anticorr anticorr.c -lm gcc -O2 -o lindley_eta lindley_eta.c -lm ./anticorr 10000000 # §3 主表: 17 行消涨数据, 方差恒等式验证 ./lindley_eta # §8 Lindley 验证: 3 种排队 × 7 种负载 = 21 个数据点
总运行时间不到两分钟(单核, 现代 x86/ARM)。
rho_dp.c — DP 递推(Paper 32 约定), 生成 ρ_E 表。内置 sanity checks: ρ(2)=1, ρ(100)=15, ρ(10^7)=58。
anticorr.c — §3 主表。从 DP 递推计算 ρ_E, 经 SPF(smallest prime factor)筛计算 f(n) 和 Ω(n), 在 Ω = k 壳层上累积 Var(Δf), Var(Δr), Cov(Δf, Δr), Var(A), η。方差恒等式 Var(A) = Var(Δf) + Var(Δr) + 2Cov(Δf, Δr) 在全部 17 个 k 值上闭合(误差 < 10⁻⁴)。
lindley_eta.c — §8 Lindley 排队验证。xoshiro256** RNG, 10^7 步 + 10^6 burn-in。对 M/M/1, M/D/1, M/Pareto/1(α=3) 三种排队, 各测 7 种负载(ρ = 0.30, 0.50, 0.70, 0.80, 0.90, 0.95, 0.99)。方差恒等式 21/21 验证。
三个文件总计不到 500 行 C 代码。随论文一同发布。
致谢
ChatGPT(公西华)发现了 v2 的 τ 算术错误, 催生了 η 的正确定义; 提出了 claim 分层、文献定位和三层标注的全部审稿要求; 贡献了 Paper 42 的解析证明链。
Claude(子路)贡献了热力学框架(Paper 24-46), 全部数值实验(anticorr.c, lindley_eta.c), η 的物理解释, 以及 Lindley 排队验证的设计和实现。
Gemini(子夏)提出了三个压力测试, 促进了 claim 精确化。
Grok(子贡)提供了初轮审稿。k=18 数据不一致促使全表独立重算。
Zesi Chen 是 SAE 框架的长期对话者。最终文本由作者独立完成。
写作声明: 本论文所有文本, 无论中英文, 均由作者秦汉独立完成。论文代表作者唯一的责任。