Implement gcd in sdk.

runtime/lib/integers.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.
 
 // TODO(srdjan): fix limitations.
 // - shift amount must be a Smi.
 class _IntegerImplementation extends _Num {
   // The Dart class _Bigint extending _IntegerImplementation requires a
   // default constructor.
 
@@ skipping 272 matching lines @@
     while (e > 0) {
       if (e.isOdd) {
         r = (r * b) % m;
       }
       e >>= 1;
       b = (b * b) % m;
     }
     return r;
   }
 
-  // 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;
-    if (m is _Bigint) {
-      return _toBigint().modInverse(m);
+  // 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 = 0;
+    if (!inv) {
+      while (x.isEven && y.isEven) {
+        x >>= 1;
+        y >>= 1;
+        s++;
+      }
+      if (y.isOdd) {
+        var t = x;
+        x = y;
+        y = t;
+      }
     }
-    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;
+    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 >>= 1;
         if (ac) {
           if (!a.isEven || !b.isEven) {
-            a += t;
-            b -= m;
+            a += y;
+            b -= x;
           }
           a >>= 1;
         } else if (!b.isEven) {
-          b -= m;
+          b -= x;
         }
         b >>= 1;
       }
       while (v.isEven) {
         v >>= 1;
         if (ac) {
           if (!c.isEven || !d.isEven) {
-            c += t;
-            d -= m;
+            c += y;
+            d -= x;
           }
           c >>= 1;
         } else if (!d.isEven) {
-          d -= m;
+          d -= x;
         }
         d >>= 1;
       }
       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 v << s;
     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;
+    if (m is _Bigint) {
+      return _toBigint().modInverse(m);
+    }
+    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;
+    if (other is _Bigint) {
+      return _toBigint().gcd(other);
+    }
+    return _binaryGcd(x, y, false);
+  }
 }
 
 class _Smi extends _IntegerImplementation implements int {
   factory _Smi._uninstantiable() {
     throw new UnsupportedError(
         "_Smi can only be allocated by the VM");
   }
   int get _identityHashCode {
     return this;
   }
@@ skipping 201 matching lines @@
   // Shift by mint exceeds range that can be handled by the VM.
   int _shrFromInt(int other) {
     if (other < 0) {
       return -1;
     } else {
       return 0;
     }
   }
   int _shlFromInt(int other) native "Mint_shlFromInt";
 }