In this paper, we present a general Mobius inversion transform formula for hcp lattices. This formula can be applied to hcp lattices with a non-ideal c/a value and to obtain the pair potential between atoms in these lattices from the cohesive energy. Also, the three-body interaction among atoms in the lattices can be taken into account in the method. This method gives a useful means to obtain interatomic interactions in the interatomic force model. The method has been applied to zinc, and the pair potential obtained is used to calculate the phonon dispersion relations for some high-symmetry directions. It is found that, by properly considering a three-body interaction, one can acquire satisfactory results.