S(k,x):=k*SUM(p^(k-1),p,1,x-1)

;User=Simp(User)
VECTOR(S(k,x),k,0,3)=[0,x-1,x*(x-1),x*(2*x^2-3*x+1)/2]