• Let $f(n)$ be an arbitrary function with integer values.
• Let $a(n)$ be an integer sequence such that $$ \frac{1}{1-x}=\sum\limits_{n=0}^{\infty}a(n)x^n\prod\limits_{k=1}^{n+1}(1-f(k)x) $$
• Start with vector $\nu$ of length $n$ with elements $\nu_i = 1$ (that is, $\nu = \{1, 1, \dotsc, 1\}$) and for $i$ from $1$ to $n-1$ and for $j$ from $i+1$ to $n$ apply $\nu_{j} = f(i)\nu_{j-1} + \nu_{j}$.

I conjecture that after the whole transform we have vector $\nu$ with elements $\nu_i = a(i-1)$.

Here is the PARI/GP program to check it numerically:

f(n) = n
upto1(n) = my(v1); v1 = vector(n, i, 1); for(i=1, n-1, for(j=i+1, n, v1[j] = f(i)*v1[j-1] + v1[j])); v1
h(n, x) = my(v1); v1 = upto1(n+1); sum(i=0, n, v1[i+1]*x^i*prod(k=1, i+1, 1-f(k)*x)) + x*O(x^n)
test(n) = my(x = 'x); Vec(h(n, x)) == vector(n+1, i, 1)

Is there a way to prove it?