It has worked in Pro since I can remember... the arcpy's differ though
import arcpy
coords = [(2524240.5113, 2562562.7555), (2524240.8774, 2562712.7555), (2524290.3774, 2562712.6347), (2524290.0113, 2562562.6347), (2524240.5113, 2562562.7555) ]
part = arcpy.Array(arcpy.Point(*coord) for coord in coords)
part
<Array [<Point (2524240.5113, 2562562.7555,
part[:2]
<Array [<Point (2524240.5113, 2562562.7555,
pnts = [p for p in part]
pnts
[<Point (2524240.5113, 2562562.7555,
arr = np.asarray([(p.X, p.Y) for p in part])
arr
array([[2524240.5113, 2562562.7555],
[2524240.8774, 2562712.7555],
[2524290.3774, 2562712.6347],
[2524290.0113, 2562562.6347],
[2524240.5113, 2562562.7555]])
np.mean(arr, axis=0)
array([2524260.45774, 2562622.70718])