A procedure is developed and tested to recover the distribution of
connectivity of an a priori unknown network, by sampling the dynamics of an
ensemble made of reactive walkers. The relative weight between reaction and
relocation is gauged by a scalar control parameter, which can be adjusted at
will. Different equilibria are attained by the system, following the externally
imposed modulation, and reflecting the interplay between reaction and diffusion
terms. The information gathered on the observation node is used to predict the
stationary density as displayed by the system, via a direct implementation of
the celebrated Heterogeneous Mean Field (HMF) approximation. This knowledge
translates into a linear problem which can be solved to return the entries of
the sought distribution. A variant of the model is then considered which
consists in assuming a localized source where the reactive constituents are
injected, at a rate that can be adjusted as a stepwise function of time. The
linear problem obtained when operating in this setting allows one to recover a
fair estimate of the underlying system size. Numerical experiments are carried
so as to challenge the predictive ability of the theory