Definir el avance operador diferencia $$\Delta f(x) = f(x+1) - f(x)$$ Yo creo que si $f(x)$ es un polinomio con coeficientes enteros, $\Delta^k f(x)$ es divisible por k!. Por la linealidad es suficiente con considerar un solo monomio $f(x) = x^n$. He comprobado esto para valores pequeños de a$n$$k$, y creen que una simple prueba de que debe de existir, pero soy incapaz de encontrarla.
En particular, la fuerza bruta da $$\Delta^k x^n = \sum_j x^{n-j} \left[ \binom{n}{j} \sum_i (-1)^{k-i} \binom{k}{i} i^j \right]$$ pero los términos entre paréntesis parecen no tener solución de forma cerrada (ver (20)-(25) de http://mathworld.wolfram.com/BinomialSums.html).
Motivación: tengo un desconocido entero coeficiente polinomio de grado $n$ de la muestra a $x = 0, 1, \ldots, n$, y quieren demostrar que todos los resultados intermedios en el clásico dividido diferencia del algoritmo son enteros.