It's worth noting this special case: if $G$ is isomorphic to $(\mathbb{Z}/p)^r$ then it can be regarded as a vector space of dimension $r$ over the field $\mathbb{Z}/p$, and the subgroups are just the subspaces. The number of linearly independent lists of length $n$ is
$(p^r-1)(p^r-p)\dotsb(p^r-p^{n-1})$
(choose a nonzero vector $v_1$, then a vector $v_2$ not in the one-dimensional space spanned by $v_1$, then $v_2$ not in the 2-dimensional space spanned by $v_1$ and $v_2$, and so on). For any subspace $V$ of dimension $n$, the number of bases is
$(p^n-1)(p^n-p)\dotsb(p^n-p^{n-1})$
by the same argument. It follows that the number $N_n$ of subspaces of dimension $n$ is
$ N_n = \frac{(p^r-1)\dotsb(p^r-p^{n-1})}{(p^n-1)\dotsb(p^n-p^{n-1})} = \frac{(p^r-1)(p^{r-1}-1)\dotsb(p-1)}{(p^n-1)(p^{n-1}-1)\dotsb(p-1)} $$ N_n = \frac{(p^r-1)\dotsb(p^r-p^{n-1})}{(p^n-1)\dotsb(p^n-p^{n-1})} = \frac{(p^r-1)(p^{r-1}-1)\dotsb(p^{r-n+1}-1)}{(p^n-1)(p^{n-1}-1)\dotsb(p-1)} $