We develop a random-sampling method for finding accurate symplectic integrators which best match the exact trajectory of a one-dimensional harmonic oscillator. We recover several well-known algorithms. We demonstrate the usefulness of the random sampling method by finding and validating a new integrator, applying it to the classical many-body problem.