On the powers of the descent set statistic

We study the sum of the $r$th powers of the descent set statistic and how many small prime factors occur in these numbers. Our results depend upon the base $p$ expansion of $n$ and $r$.