Rekurencyjna funkcja skumulowana

Muszę napisać skumulowaną funkcję sumowania w R, ale uderzyłem w ścianę z cegieł. Funkcja ma następującą strukturę:

a*x1
a*x2 + a^2*x1
a*x3 + a^2*x2 + a^3*x1
a*x4 + a^2*x3 + a^3*x2 + a^4*x1

I tak dalej. Wydaje się, że cumsum nie działa dla tego typu funkcji. Czy jest jakiś sposób, w jaki można to zaimplementować w R?