On the powers of the descent set statistic

We study the sum of the rth 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.