  We use Monte-Carlo simulations, combined with homogeneously determined age
and mass distributions based on multi-wavelength photometry, to constrain the
cluster formation history and the rate of bound cluster disruption in the Large
Magellanic Cloud (LMC) cluster system. We evolve synthetic star cluster systems
formed with a power-law initial cluster mass function (ICMF) of spectral index
$\alpha =-2$ assuming different cluster disruption time-scales. For each of
these disruption time-scales we derive the corresponding cluster formation rate
(CFR) required to reproduce the observed cluster age distribution. We then
compare, in a Poissonian $\chi^2$ sense, model mass distributions and model
two-dimensional distributions in log(mass) vs. log(age) space of the detected
surviving clusters to the observations. Because of the bright detection limit
($M_V^{\rm lim} \simeq -4.7$ mag) above which the observed cluster sample is
complete, one cannot constrain the characteristic disruption time-scale for a
$10^4$ M$_\odot$ cluster, $t_4^{\rm dis}$ (where the disruption time-scale
depends on cluster mass as $t_{\rm dis} = t_4^{\rm dis} (M_{\rm cl} / 10^4 {\rm
M}_\odot)^0.62$), to better than $t_4^{\rm dis} \ge 1$ Gyr. We conclude that
the CFR has increased from 0.3 clusters Myr$^{-1}$ 5 Gyr ago, to a present rate
of $(20-30)$ clusters Myr$^{-1}$. For older ages the derived CFR depends
sensitively on our assumption of the underlying CMF shape. If we assume a
universal Gaussian ICMF, then the CFR has increased steadily over a Hubble time
from $\sim 1$ cluster Gyr$^{-1}$ 15 Gyr ago to its present value. If the ICMF
has always been a power law with a slope close to $\alpha=-2$, the CFR exhibits
a minimum some 5 Gyr ago, which we tentatively identify with the well-known age
gap in the LMC's cluster age distribution.
