Inertia plays a crucial role in rigid body dynamics, and the associated mass matrix is of various forms in representation. The influences of accuracy of different representation, however, have not drawn enough attention in previous researches about numerical simulation of rigid body dynamics. In the paper, the inertia representation is intensively investigated for rigid body dynamics and a modified formulation is derived through splitting the kinetic energy into two parts: a square term of velocity and a quadratic form in the derivatives with quadratic coefficients in generalized displacement, of which the proportion is controlled by a scaling parameter. Although the kinetic energy with different scaling parameters is theoretically equivalent in dynamics, error estimation demonstrates that accuracy of numerical scheme crucially depends on the particular value of scaling parameter if only rotational coordinates are expressed in pseudo vectors. This attractive feature distinguishes the modified formulation from others in numerical significance. According to the modified representation of inertia, a variational integrator is derived for rigid body dynamics in pseudo vectors. Numerical results demonstrate that the variational integrator, of which the scaling parameter is selected as the arithmetic mean of three principal moments of inertia tensor, is of impressively higher accuracy in simulation, especially compared with the integrations derived with the original formulation of mass matrix.