From the Gauss-Bonnet-Chern theorem, the Euler characteristic of NUT-Kerr-Newman black hole is calculated to be some discrete numbers from 0 to 2. We find that the Bekenstein-Hawking entropy is the largest entropy in topology by taking into account of the relationship between the entropy and the Euler characteristic. The NUT-Kerr-Newman black hole evolves from the torus-like topological structure to the spherical structure with the changes of mass, angular momentum, electric and NUT charges. In this process, the Euler characteristic and the entropy are changed discontinuously, which give the topological aspect of the first-order phase transition of NUT-Kerr-Newman black hole. The corresponding latent heat of the topological phase transition is also obtained. The estimated latent heat of the black hole evolving from the star just lies in the range of the energy of gamma ray bursts.