We calculate the evolution of 80 - 140 Msun stars that undergo pair-instability pulsations during oxygen burning for various masses and metallicities.  We find that the presupernova mass ejection caused by pulsations is more extensive for more massive stars.  Some ejecta includes a large amount of C+O core materials.  We calculate the subsequent evolution through Fe core-collapse.  We then construct possible explosion models including circumstellar interaction and apply for a certain type of superluminous supernovae.  We also discuss the maximum mass of black holes in this mass range to compare with the gravitational wave observations.