We recently proposed the binary interaction approximation (BIA) to N-body problems, which, in principle, excludes the interparticle force evaluation if the exact solutions are known for the corresponding two-body problems such as the Coulombic and gravitational interactions. In this article, a detailed introduction to the BIA is given, including the error analysis to give the expressions for the approximation error in the total angular momentum and the total energy of the entire system. It is shown that, although the energy conservation of the BIA scheme is worse than the 4th order Hermite integrator (HMT4) for similar elapsed, or the wall-clock times, the individual errors in position and in velocity are much better than HMT4. The energy error correction scheme to the BIA is also introduced that does not deteriorate the individual errors in position and in velocity. It is suggested that the BIA scheme is applicable to the tree method, the particle-mesh (PM), and the particle-particle-particle-mesh (PPPM) schemes simply by replacing the force evaluation and the conventional time integrator with the BIA scheme. (C) 2015 Elsevier B.V. All rights reserved.