Added ::GCD and ::LCM: <numeric> [c++17] contains gcd and lcm impleme…

…ntation (#20583)

src/operator/mshadow_op.h

@@ -37,6 +37,17 @@
 #include <cuda_fp16.h>
 #endif
 
+#define MXNET_HAS_GCD_LCM() 0
+#if __cplusplus >= 201703L
+#ifdef __has_gcd_lcm
+#if __has_gcd_lcm(<numeric>)
+#include <numeric>
+#undef MXNET_HAS_GCD_LCM
+#define MXNET_HAS_GCD_LCM() 1
+#endif
+#endif
+#endif
+
 namespace mxnet {
 namespace op {
 namespace mshadow_op {
@@ -1808,6 +1819,9 @@ struct gcd : public mxnet_op::tunable {
   template <typename DType>
   MSHADOW_XINLINE static typename enable_if<is_integral<DType>::value, DType>::type Map(DType a,
                                                                                         DType b) {
+#if MXNET_HAS_GCD_LCM()
+    return std::gcd(a, b);
+#else
     // minus cases.
     if (a < 0) {
       a = -a;
@@ -1839,6 +1853,7 @@ struct gcd : public mxnet_op::tunable {
       c = b;
     }
     return c;
+#endif
   }
 
   template <typename DType>
@@ -1853,6 +1868,9 @@ struct lcm : public mxnet_op::tunable {
   template <typename DType>
   MSHADOW_XINLINE static typename enable_if<is_integral<DType>::value, DType>::type Map(DType a,
                                                                                         DType b) {
+#if MXNET_HAS_GCD_LCM()
+    return std::lcm(a, b);
+#else
     // minus cases.
     if (a < 0) {
       a = -a;
@@ -1882,6 +1900,7 @@ struct lcm : public mxnet_op::tunable {
       c = tmp_a / b * tmp_b;
     }
     return c;
+#endif
   }
   template <typename DType>
   MSHADOW_XINLINE static typename enable_if<!is_integral<DType>::value, DType>::type Map(DType a,