  We present analytical results for the emerging structure of networks that
evolve via a combination of growth (by node addition and random attachment) and
contraction (by random node deletion). To this end we consider a network model
in which at each time step a node addition and random attachment step takes
place with probability $P_{add}$ and a random node deletion step takes place
with probability $P_{del}=1-P_{add}$. The balance between the growth and
contraction processes is captured by the parameter $\eta=P_{add}-P_{del}$. The
case of pure network growth is described by $\eta=1$. In case that $0<\eta<1$
the rate of node addition exceeds the rate of node deletion and the overall
process is of network growth. In the opposite case, where $-1<\eta<0$, the
overall process is of network contraction, while in the special case of
$\eta=0$ the expected size of the network remains fixed, apart from
fluctuations. Using the master equation we obtain a closed form expression for
the time dependent degree distribution $P_t(k)$. The degree distribution
$P_t(k)$ includes a term that depends on the initial degree distribution
$P_0(k)$, which decays as time evolves, and an asymptotic distribution
$P_{st}(k)$. In the case of pure network growth ($\eta=1$) the asymptotic
distribution $P_{st}(k)$ follows an exponential distribution, while for
$-1<\eta<1$ it consists of a sum of Poisson-like terms and exhibits a
Poisson-like tail. In the case of overall network growth ($0 < \eta < 1$) the
degree distribution $P_t(k)$ eventually converges to $P_{st}(k)$. In the case
of overall network contraction ($-1 < \eta < 0$) we identify two different
regimes. For $-1/3 < \eta < 0$ the degree distribution $P_t(k)$ quickly
converges towards $P_{st}(k)$. In contrast, for $-1 < \eta < -1/3$ the
convergence of $P_t(k)$ is initially very slow and it gets closer to
$P_{st}(k)$ only shortly before the network vanishes.
