Solution.
Observe that g(z)=e−2πiz is such a function. Given any other f(z) satisfying the same identities, set h(z)=g(z)f(z). Then
h(z+1)=g(z+1)f(z+1)=g(z)f(z)=h(z)h(z+i)=g(z+i)f(z+i)=e2πg(z)e2πf(z)=h(z).
In other words, h is doubly periodic. Notice that g(z) never vanishes, so h is also entire. In particular, h is continuous on [0,1]×[0,1], so by double periodicity, h is bounded everywhere. By Liouville's theorem, h is constant, so f(z)=Ce−2πiz for some C∈C, which completes the proof.