After this question popped up again, it seemed to me to scream out for a use of the Fourier Transform (FT). I have decided to post this as an answer, since this approach is transparent and provides a systematic method to obtain all (distributional) solutions. A simple formal manipulation (to be justified below) shows that the FT $g$ of a solution satisfies the equation $$e^{\{-iy+1+iy\}}g(y)=0.$$ In classical terms, this only provides the zero solution but in the sense of distributions it has many non-trivial ones, in fact, suitable combinations of $\delta$-functions with singularities at the zeros of the bracket function in brackets. This easily translates into a precise description of all possible solutions of the original equation. In order to justify this, a few words on the FT for distributions are, perhaps, in order. The problem of extending the FT to distributions, which was motivated by applications in physics, was studied in detail at the middle of the last century. L. Schwartz considered the case of tempered distributions in his seminal monograph. But this can be done for many other situations. The basic starting point for such an extension is an initial setting where it is an isomorphism between two l.c. spaces of test functions (that is, functions with good smoothness and/or growth properties). One then uses transposes to translate this into an isomorphism between corresponding spaces of distributions. Schwartz used the (symmetric) case of the smooth functions of rapid decrease but there are many non symmetric cases which are of interest. We require the FT’s for arbitrary distributions on the line and so used the fact that it is an isomorphism between the smooth functions on the line with compact support and a suitable space of entire functions of exponential growth (known as the Paley-Schwartz theorem--details are in Strichartz). Dualising, we get a version of the FT which works for ALL distributions--it takes its values in a space of analytic functionals (a superspace of the dual space of the Fr´echet space of entire functions). This can be described explicitly, but for our purposes it suffices to know that it contains the delta functions on the complex plane, together with suitable combinations.