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.