  Diffusion on graphs is ubiquitous with numerous high-impact applications. In
these applications, complete diffusion histories play an essential role in
terms of identifying dynamical patterns, reflecting on precaution actions, and
forecasting intervention effects. Despite their importance, complete diffusion
histories are rarely available and are highly challenging to reconstruct due to
ill-posedness, explosive search space, and scarcity of training data. To date,
few methods exist for diffusion history reconstruction. They are exclusively
based on the maximum likelihood estimation (MLE) formulation and require to
know true diffusion parameters. In this paper, we study an even harder problem,
namely reconstructing Diffusion history from A single SnapsHot} (DASH), where
we seek to reconstruct the history from only the final snapshot without knowing
true diffusion parameters. We start with theoretical analyses that reveal a
fundamental limitation of the MLE formulation. We prove: (a) estimation error
of diffusion parameters is unavoidable due to NP-hardness of diffusion
parameter estimation, and (b) the MLE formulation is sensitive to estimation
error of diffusion parameters. To overcome the inherent limitation of the MLE
formulation, we propose a novel barycenter formulation: finding the barycenter
of the posterior distribution of histories, which is provably stable against
the estimation error of diffusion parameters. We further develop an effective
solver named DIffusion hiTting Times with Optimal proposal (DITTO) by reducing
the problem to estimating posterior expected hitting times via the
Metropolis--Hastings Markov chain Monte Carlo method (M--H MCMC) and employing
an unsupervised graph neural network to learn an optimal proposal to accelerate
the convergence of M--H MCMC. We conduct extensive experiments to demonstrate
the efficacy of the proposed method.
