🩹 stats: fix [g]zscore shape-typing on numpy<2.1

Author: jorenham

scipy-stubs/stats/_stats_py.pyi

@@ -689,14 +689,11 @@ def sem(
 # NOTE: keep in sync with `gzscore`
 @overload  # +integer, known shape
 def zscore(
-    a: onp.CanArray[_ShapeT, np.dtype[npc.integer | np.bool_]],
-    axis: int | None = 0,
-    ddof: int = 0,
-    nan_policy: NanPolicy = "propagate",
+    a: onp.CanArrayND[npc.integer | np.bool_, _ShapeT], axis: int | None = 0, ddof: int = 0, nan_policy: NanPolicy = "propagate"
 ) -> onp.ArrayND[np.float64, _ShapeT]: ...
 @overload  # known inexact dtype, known shape
 def zscore(
-    a: onp.CanArray[_ShapeT, np.dtype[_InexactT]], axis: int | None = 0, ddof: int = 0, nan_policy: NanPolicy = "propagate"
+    a: onp.CanArrayND[_InexactT, _ShapeT], axis: int | None = 0, ddof: int = 0, nan_policy: NanPolicy = "propagate"
 ) -> onp.ArrayND[_InexactT, _ShapeT]: ...
 @overload  # float 1d
 def zscore(
@@ -734,15 +731,15 @@ def zscore(
 # NOTE: keep in sync with `zscore`
 @overload  # +integer, known shape
 def gzscore(
-    a: onp.CanArray[_ShapeT, np.dtype[npc.integer | np.bool_]],
+    a: onp.CanArrayND[npc.integer | np.bool_, _ShapeT],
     *,
     axis: int | None = 0,
     ddof: int = 0,
     nan_policy: NanPolicy = "propagate",
 ) -> onp.ArrayND[np.float64, _ShapeT]: ...
 @overload  # known inexact dtype, known shape
 def gzscore(
-    a: onp.CanArray[_ShapeT, np.dtype[_InexactT]], *, axis: int | None = 0, ddof: int = 0, nan_policy: NanPolicy = "propagate"
+    a: onp.CanArrayND[_InexactT, _ShapeT], *, axis: int | None = 0, ddof: int = 0, nan_policy: NanPolicy = "propagate"
 ) -> onp.ArrayND[_InexactT, _ShapeT]: ...
 @overload  # float 1d
 def gzscore(