For a binary integer program (IP) max c T x, Ax ≤ b, x ∈ {0, 1} n , where A ∈ R m×n and c ∈ R n have independent Gaussian entries and the right-hand side b ∈ R m satisfies that its negative coordinates have ℓ2 norm at most n/10, we prove that the gap between the value of the linear programming relaxation and the IP is upper bounded by poly(m)(log n) 2 /n with probability at least 1 − 1/n 7 − 1/2 Ω(m) . Our results give a Gaussian analogue of the classical integrality gap result of Dyer and Frieze (Math. of O.R., 1989) in the case of random packing IPs. In constrast to the packing case, our integrality gap depends only polynomially on m instead of exponentially. By recent breakthrough work of Dey, Dubey and Molinaro (SODA, 2021), the bound on the integrality gap immediately implies that branch and bound requires n poly(m) time on random Gaussian IPs with good probability, which is polynomial when the number of constraints m is fixed.