  Using an artificial neutral network we explore the parameter space of
supergravity grand unified models consistent with the combined Fermilab E989
and Brookhaven E821 data on $(g-2)_\mu$. The analysis indicates that the region
favored by the data is the one generated by gluino-driven radiative breaking of
the electroweak symmetry. This region naturally leads to a split sparticle
spectrum with light sleptons and weakinos but heavy squarks, with the stau and
the chargino as the lightest charged particles. We show that if the entire
deviation from the standard model $(g-2)_{\mu}$ arises from supersymmetry, then
supersymmetry is discoverable at HL-LHC and HE-LHC via production and decay of
sleptons within the optimal integrated luminosity of HL-LHC and with a smaller
integrated luminosity at HE-LHC.
