Given multiple polygons that overlap in multiple ways, I would like to export from these features all polygons that don't overlap with others, iteratively. 

The product would be a number of features with no overlap that, when summed together, make up the original. 

The products could then be used as input to Zonal Statistics, and this would be much faster than iterating Zonal Statistics over each polygon. 

I have been trying to code this in ArcPy without success. 

Does code to do this already exist?