We study the recovery of a finitely supported distribution, a complex linear combination of Dirac measures, from intensity measurements. The distribution µ = s j=1 c j δt j is given by a coefficient vector c ∈ C s and its support {t 1 , t 2 , . . . , ts} is contained in [0, Λ] for some Λ > 0. The intensity measurements evaluate (squared) magnitudes of a set of linear functionals applied to µ, obtained by sampling μ, the Fourier transform of µ, or by evaluating differences between modulated samples. Following a strategy by Alexeev et al., the structure of the linear functionals, and hence of the non-linear magnitude measurement, is encoded with a graph, where the vertices represent point evaluations of μ at {v 1 , v 2 , . . . , vn} ⊂ [−Ω, Ω] and each edge represents a (modulated) difference between vertices incident with it. We show that a Ramanujan graph with degree d ≥ 3 and n > 6(1+6/ ln(s/ΛΩ))svertices provides M = (d + 1)n magnitudes that are sufficient for identifying the complex measure up to an overall unimodular multiplicative constant. At the cost of including an additional oversampling step and with an additional requirement that n − 1 is prime, we construct an explicit recovery algorithm that is based on the Prony method.