In this paper, we shall utilize Nevanlinna value distribution theory to study the solvability of the difference equations of the form: f (z)(n) + p(z)(Delta(f)(c))(m) = r(z)e(q(z)) and f (z)(n) + p(z)e(q(z))(Delta(c)f)(m) = r(z), and we shall study the growth of their entire solutions. Moreover, we will give a number of examples to show that the results in this paper are the best possible in certain senses. This article extends earlier results by Liu et al. (Czechoslov. Math. J.61:565-576, 2011; Ann. Pol. Math. 102:129-142, 2011).