This study proposes a mathematical model of clustered multi-temperature joint distribution in fuzzy environment. In this model, a Z-shaped function is used to depict customer satisfaction. For the imprecise model, triangular fuzzy numbers are used to represent travel times. By redefining the movement procedure of fireflies, two versions of discrete firefly algorithms are developed, which differ in the population initialization strategy. Finally, experiments are carried out and computational results are reported to confirm the effectiveness of the proposed algorithms.