On a stable variational approximation of the cut locus, and a non-local isoperimetric problem