This paper introduces three different simulation algorithms of the shifted Poisson distribution. The first algorithm is the inverse transform method, the second is the rejection sampling, and the third is gamma-Poisson hierarchy sampling. Three algorithms have different regions of parameters at which they are efficient. We numerically compare those algorithms with different sets of parameters. As an application, we give a simulation method of the constant elasticity of variance model.