Magnetar is a kind of pulsar powered by magnetic field energy. The study of magnetars is an important hotspot in the field of pulsars. In this paper, according to the work of Zhu Cui, et al. (Zhu C, Gao Z F, Li X D, Wang N, Yuan J P, Peng Q H 2016 Mod. Phys. Lett. A 311650070 ), we reinvestigate the Landau-level stability of electrons in a superhigh magnetic field (SMF), $B\gg B_{\rm cr}$ ( B _cris a quantum critical magnetic field with a value of 4.414×10 ¹³G), and its influence on the pressure of electrons in magnetar. First, we briefly review the pressure of electrons in neutron star (NS) with a weak-magnetic field limit ( $ B\ll B $ _cr). Then, we introduce an electron Landau level stability coefficient g _νand a Dirac- δfunction to deduce a modified pressure formula for the degenerate and relativistic electrons in an SMF in an application range of matter density ρ≥ 10 ⁷g·cm ^–3and B _cr $ \ll $ B< 10 ¹⁷G. By modifying the phase space of relativistic electrons, the SMF can enhance the electron number density n _e, and reduce the maximum of electron Landau level number ν _max, which results in a redistribution of electrons. As Bincreases, more and more electrons will occupy higher Landau levels, and the electron Landau level stability coefficient g _νwill decrease with the augment of Landau energy-level number ν. By modifying the phase space of relativistic electrons, the electron number density n _eincreases with the MF strength increasing, leading the electron pressure P _eto increase. Utilizing the modified expression of electron pressure, we discuss the phenomena of Fermion spin polarization and electron magnetization in the SMF, and the modification of the equation of state by the SMF. We calculate the baryon number density, magnetization pressure, and the difference between pressures in the direction parallel to and perpendicular to the magnetic field in the frame of the relativistic mean field model. Moreover, we find that the pressure anisotropy due to the strong magnetic field is very small and can be ignored in the present model. We compare our results with the results from other similar studies, and examine their similarities and dissimilarities. The similarities include 1) the abnormal magnetic moments of electrons and the interaction between them are ignored; 2) the electron pressure relate to magnetic field intensity B, electron number density n _eand electron Fermi energy $E_{{\rm{F}}}^{{\rm{e}}}$ , and the latter two are complex functions containing B; 3) with n _eand $E_{{\rm{F}}}^{{\rm{e}}}$ fixed, P _eincreases with Brising; 4) as Bincreases, the pressure-density curves fitted by the results from other similar studies have irregular protrusions or fluctuations, which are caused by the transformation of electron energy state from partial filling to complete filling at the ν-level or the transition of electrons from the νto the ( ν+1)-level. This phenomenon is believed to relate to the behavior of electrons near the Fermi surface in a strong magnetic field, which essentially reflects the Landau level instability. Finally, the future research direction is prospected. The present results provide a reference for future studies of the equation of state and emission mechanism of high- Bpulsar, magnetar and strongly magnetized white dwarf.