Las ecuaciones funcionales lineales pueden resolverse con la transformada de Fourier. Sea $\lambda_k$ sean las raíces de la ecuación $e^\lambda-1=\lambda$ . Hay infinitas raíces de este tipo. Entonces
$$f(z)=\sum_k a_ke^{\lambda_k z}$$
es una solución. Aquí la suma puede ser finita, y $a_k$ arbitraria, o la suma puede ser infinita, y $a_k$ tienden a cero con tal rapidez que la serie converge.
La solución trivial está cubierta, si interpretas la respuesta correctamente. La ecuación $e^\lambda-1=\lambda$ tiene una DOBLE raíz en $0$ . Para el caso de una raíz doble se incluye no sólo el término $e^{\lambda z}$ sino también el término $ze^{\lambda z}$ . Del mismo modo, para las raíces múltiples $z^ke^{\lambda z}$ . Así que la raíz $\lambda=0$ cubre exactamente la solución $az+b$ .
EDIT: En realidad todas las soluciones enteras se pueden representar de esta manera, para la prueba me remito a Gelfond, Calculus of finite differences, MR0342890, Chap. 5 Sect 7, Thm II and Corollaries. Él da la ecuación $f'(z)=f(z-1)$ como ejemplo, pero su ecuación se trata de forma similar.
EDIT2: El problema se puede generalizar como sigue: Sea $\omega$ sea una distribución con soporte limitado. Ecuación $f\star w=0$ , donde $\star$ es la convolución, se llama ecuación de convolución y sus soluciones se denominan media-periódica funciones (fonctions moyenne-periodiques). El caso que nos ocupa es $\omega(x)=\delta(x+1)-\delta(x)-\delta'(x)$ donde $\delta$ es la función delta. La teoría de las funciones periódicas medias fue creada por Delsarte y Schwartz en 1940. Schwartz tiene un resultado general de que todas las funciones periódicas medias pueden obtenerse como límite de sumas exponenciales, como en este problema.