We derive the multi-fractal scaling of probability distributions of multi-particle configurations for the binary reaction-diffusion system A + A → ∅ in d ≤ 2 and for the ternary system 3A → ∅ in d = 1. For the binary reaction we find that the probability Pt(N, ∆V ) of finding N particles in a fixed volume element ∆V at time t decays in the limit of large time as (For the ternary reaction in one dimension we find that Pt(N, ∆V ) ∼ (. The principal tool of our study is the dynamical renormalization group. We compare predictions of ε-expansions for Pt(N, ∆V ) for binary reaction in one dimension against exact known results. We conclude that the ε-corrections of order two and higher are absent in the above answer for Pt(N, ∆V ) for N = 1, 2, 3, 4. Furthermore we conjecture the absence of ε 2 -corrections for all values of N .