In this paper, we design a precise integration method with a variable time step for the general three-body problem that maintains all the conserved quantities. Our method is based on a logarithmic Hamiltonian leapfrog with chain vectors proposed by Mikkola & Tanikawa and features an energy-preserving parameter. Although the proposed method is merely second-order accurate, it can precisely trace some periodic orbits. This is not possible with generic geometric eighth-order numerical integrators and the logarithmic Hamiltonian leapfrog approach. Further, similar to logarithmic Hamiltonian leapfrog, our method is analytically shown to have Lagrangian solutions. Prior to the presented integration method, no integration method was known to preserve all the conserved quantities, in addition to presenting triangular Lagrangian solutions. Because our method is implicit, it requires an iteration method. Therefore, the proposed approach seems to be computationally intensive. However, our method is less computationally burdensome than a generic explicit eighth-order symplectic method.