We extend covariance velocity map ion imaging to four particles, establishing cumulant mapping and allowing for measurements that provide insights usually associated with coincidence detection, but at much higher count rates. Without correction, a fourfold covariance analysis is contaminated by the pairwise correlations of uncorrelated events, but we have addressed this with the calculation of a full cumulant, which subtracts pairwise correlations. We demonstrate the approach on the four-body breakup of formaldehyde following strong field multiple ionization in few-cycle laser pulses. We compare Coulomb explosion imaging for two different pulse durations (30 and 6 fs), highlighting the dynamics that can take place on ultrafast timescales. These results have important implications for Coulomb explosion imaging as a tool for studying ultrafast structural changes in molecules, a capability that is especially desirable for high-count-rate x-ray free-electron laser experiments.