An efficient integration method for binaries in N-body simulations