Implement gcd in sdk.

sdk/lib/_internal/compiler/js_lib/js_number.dart

 // Copyright (c) 2012, the Dart project authors.  Please see the AUTHORS file
 // for details. All rights reserved. Use of this source code is governed by a
 // BSD-style license that can be found in the LICENSE file.
 
 part of _interceptors;
 
 /**
  * The super interceptor class for [JSInt] and [JSDouble]. The compiler
  * recognizes this class as an interceptor, and changes references to
  * [:this:] to actually use the receiver of the method, which is
@@ skipping 392 matching lines @@
     while (e > 0) {
       if (e.isOdd) {
         r = (r * b) % m;
       }
       e ~/= 2;
       b = (b * b) % m;
     }
     return r;
   }
 
-  // Returns 1/this % m, with m > 0.
-  int modInverse(int m) {
-    if (m is! int) throw argumentErrorValue(m);
-    if (m <= 0) throw new RangeError(m);
-    if (m == 1) return 0;
-    int t = this;
-    if ((t < 0) || (t >= m)) t %= m;
-    if (t == 1) return 1;
-    final bool ac = m.isEven;
-    if ((t == 0) || (ac && t.isEven)) throw new RangeError("Not coprime");
-    int u = m;
-    int v = t;
+  // If inv is false, returns gcd(x, y).
+  // If inv is true and gcd(x, y) = 1, returns d, so that c*x + d*y = 1.
+  // If inv is true and gcd(x, y) != 1, throws RangeError("Not coprime").
+  static int _binaryGcd(int x, int y, bool inv) {
+    int s = 1;
+    if (!inv) {
+      while (x.isEven && y.isEven) {
+        x ~/= 2;
+        y ~/= 2;
+        s *= 2;
+      }
+      if (y.isOdd) {
+        var t = x;
+        x = y;
+        y = t;
+      }
+    }
+    final bool ac = x.isEven;
+    int u = x;
+    int v = y;
     int a = 1,
         b = 0,
         c = 0,
         d = 1;
     do {
       while (u.isEven) {
         u ~/= 2;
         if (ac) {
           if (!a.isEven || !b.isEven) {
-            a += t;
-            b -= m;
+            a += y;
+            b -= x;
           }
           a ~/= 2;
         } else if (!b.isEven) {
-          b -= m;
+          b -= x;
         }
         b ~/= 2;
       }
       while (v.isEven) {
         v ~/= 2;
         if (ac) {
           if (!c.isEven || !d.isEven) {
-            c += t;
-            d -= m;
+            c += y;
+            d -= x;
           }
           c ~/= 2;
         } else if (!d.isEven) {
-          d -= m;
+          d -= x;
         }
         d ~/= 2;
       }
       if (u >= v) {
         u -= v;
         if (ac) a -= c;
         b -= d;
       } else {
         v -= u;
         if (ac) c -= a;
         d -= b;
       }
     } while (u != 0);
+    if (!inv) return s*v;
     if (v != 1) throw new RangeError("Not coprime");
     if (d < 0) {
-      d += m;
-      if (d < 0) d += m;
-    } else if (d > m) {
-      d -= m;
-      if (d > m) d -= m;
+      d += x;
+      if (d < 0) d += x;
+    } else if (d > x) {
+      d -= x;
+      if (d > x) d -= x;
     }
     return d;
   }
 
+  // Returns 1/this % m, with m > 0.
+  int modInverse(int m) {
+    if (m is! int) throw new ArgumentError(m);
+    if (m <= 0) throw new RangeError(m);
+    if (m == 1) return 0;
+    int t = this;
+    if ((t < 0) || (t >= m)) t %= m;
+    if (t == 1) return 1;
+    if ((t == 0) || (t.isEven && m.isEven)) throw new RangeError("Not coprime");
+    return _binaryGcd(m, t, true);
+  }
+
+  // Returns gcd of abs(this) and abs(other), with this != 0 and other !=0.
+  int gcd(int other) {
+    if (other is! int) throw new ArgumentError(other);
+    if ((this == 0) || (other == 0)) throw new RangeError(0);
+    int x = this.abs();
+    int y = other.abs();
+    if ((x == 1) || (y == 1)) return 1;
+    return _binaryGcd(x, y, false);
+  }
+
   // Assumes i is <= 32-bit and unsigned.
   static int _bitCount(int i) {
     // See "Hacker's Delight", section 5-1, "Counting 1-Bits".
 
     // The basic strategy is to use "divide and conquer" to
     // add pairs (then quads, etc.) of bits together to obtain
     // sub-counts.
     //
     // A straightforward approach would look like:
     //
@@ skipping 34 matching lines @@
 }
 
 class JSDouble extends JSNumber implements double {
   const JSDouble();
   Type get runtimeType => double;
 }
 
 class JSPositiveInt extends JSInt {}
 class JSUInt32 extends JSPositiveInt {}
 class JSUInt31 extends JSUInt32 {}