内側のsumと外側のsumを入れ替え、sum(k^a,k,1,n)をs(a)に置き換えると（再帰的に定義）、

s[p_]:=s[p]=Simplify[((n+1)^(p+1)-1-Sum[Binomial[p+1,a] s[a],{a,0,p-1}])/(p+1)]

となります。

関数Simplifyは、適当に括弧をはずし、適当に因数分解してくれるような感じの関数です。
何度も呼び出す可能性がある場合は、メモリに記憶させるとよく、途中に、“=s[p]”を入れておきます。

これで、s[100]を計算させると、

(-94598037819122125295227433069493721872702841533066936133385696204311395415197247711*n)/33330
+ 18674771685049011296614057325260991542019103825338734064995339343748060441015725*n^3
- (16293234618989521508515025064456465992824384487957638029599182473343901462949018943*n^5)/442
(中略)
(43232541100*n^91)/3 - 66698170*n^93 + 298760*n^95 - (2695*n^97)/2 + (25*n^99)/3 + n^100/2 + n^101/101

となりました。上の定義式が、希望されているものなのではないでしょうか？
以上、Mathematicaでの書式です。適当に記述をmaximaに変換してお試し下さい。