Variational Integrators and Graph-Based Solvers for Multibody Dynamics in Maximal Coordinates