In case you're wondering how to get differentials to act like operators in Mathematica, I stumbled across a package Carl Woll made to solve this issue in this question. There's a a more recent version of the package in his github than is linked in his comment. Here's how to write it using the package:
$\left(f\left(x\right)\frac{\partial}{\partial x}\right)^nf\left(x\right)$=Simplify[((Subscript[operator[DifferentialOperator[]], x] f[x])^n)[f[x]]]
What I want is a summation which can give the same result. I tried telling Mathematica to assume n was a positive integer but it didn't simplify any more. Can anyone help me?
Example for n=3:
$f(x)(f'(x)^3+4f(x)f'(x)f''(x)+f(x)^2f'''(x))$