Simulation-based population inference of LISA's Galactic binaries: Bypassing the global fit