I recently discovered that the following identity is true:
$$ \boxed{\frac{\partial^{n+1}}{\partial x^{n+1}}\left(\frac{(xy-1)^n}{n!} \log \frac{1}{1-x}\right) \equiv \frac{y^{n+1}}{1-xy} \pmod{(y-1)^{n+1}}} $$
To add a bit more context, for any polynomial $f(x)$ of degree at most $n$, we could apply the linear functional $T(y^i) = f(i)$ to the formula above, and it will yield the result about the values of $f(x)$, specifically
$$ \frac{\partial^{n+1}}{\partial x^{n+1}}\left( \sum\limits_{i=0}^n \frac{(-1)^{n-i} f(i)}{i!(n-i)!} x^i \right) \left(\sum\limits_{j=1}^{\infty} \frac{x^j}{j} \right) = \sum\limits_{k=1}^\infty f(n+k) x^{k-1}. $$
It then allows, for a given $f(0), \dots, f(n)$ and integer $s$ to find $f(s),\dots, f(s+n)$ in a single convolution.
The result above can as well be obtained by somewhat mindless formula bash starting from the Lagrange interpolation formula, but I'd really like to know if there is any deeper meaning to the expression above in terms of $x$ and $y$.