  Recent experimental progress has made it possible to detect in real-time
single electrons tunneling through Coulomb blockade nanostructures, thereby
allowing for precise measurements of the statistical distribution of the number
of transferred charges, the so-called full counting statistics. These
experimental advances call for a solid theoretical platform for equally
accurate calculations of distribution functions and their cumulants. Here we
develop a general framework for calculating zero-frequency current cumulants of
arbitrary orders for transport through nanostructures with strong Coulomb
interactions. Our recursive method can treat systems with many states as well
as non-Markovian dynamics. We illustrate our approach with three examples of
current experimental relevance: bunching transport through a two-level quantum
dot, transport through a nano-electromechanical system with dynamical
Franck-Condon blockade, and transport through coherently coupled quantum dots
embedded in a dissipative environment. We discuss properties of high-order
cumulants as well as possible subtleties associated with non-Markovian
dynamics.
