Problem: There is an infinite board of square tiles. Initially all tiles are white.

Vova has a red marker and a blue marker. Red marker can color a tiles. Blue marker can color b tiles. If some tile isn't white then you can't use marker of any color on it. Each marker must be drained completely, so at the end there should be exactly a red tiles and exactly b blue tiles across the board.

Vova wants to color such a set of tiles that:

  * they would form a rectangle, consisting of exactly a+b colored tiles; 
  * all tiles of at least one color would also form a rectangle. 



Here are some examples of correct colorings:

<image>

Here are some examples of incorrect colorings:

<image>

Among all correct colorings Vova wants to choose the one with the minimal perimeter. What is the minimal perimeter Vova can obtain?

It is guaranteed that there exists at least one correct coloring.

Input

A single line contains two integers a and b (1 ≤ a, b ≤ 10^{14}) — the number of tiles red marker should color and the number of tiles blue marker should color, respectively.

Output

Print a single integer — the minimal perimeter of a colored rectangle Vova can obtain by coloring exactly a tiles red and exactly b tiles blue.

It is guaranteed that there exists at least one correct coloring.

Examples

Input

4 4


Output

12


Input

3 9


Output

14


Input

9 3


Output

14


Input

3 6


Output

12


Input

506 2708


Output

3218

Note

The first four examples correspond to the first picture of the statement.

Note that for there exist multiple correct colorings for all of the examples.

In the first example you can also make a rectangle with sides 1 and 8, though its perimeter will be 18 which is greater than 8.

In the second example you can make the same resulting rectangle with sides 3 and 4, but red tiles will form the rectangle with sides 1 and 3 and blue tiles will form the rectangle with sides 3 and 3.

Solution (in CPP):

Solution:
#include <bits/stdc++.h>
using namespace std;
int mod = 1000000007;
const int inf = 1034567891;
const long long LL_INF = 1234567890123456789ll;
template <typename Arg1>
void __f(const char* name, Arg1&& arg1) {
  cout << name << " : " << arg1 << '\n';
}
template <typename Arg1, typename... Args>
void __f(const char* names, Arg1&& arg1, Args&&... args) {
  const char* comma = strchr(names + 1, ',');
  cout.write(names, comma - names) << " : " << arg1 << " | ";
  __f(comma + 1, args...);
}
template <typename T>
T GCD(T a, T b) {
  long long t;
  while (a) {
    t = a;
    a = b % a;
    b = t;
  }
  return b;
}
template <typename T>
string toString(T a) {
  return to_string(a);
}
template <typename T>
void toInt(string s, T& x) {
  stringstream str(s);
  str >> x;
}
inline int add(int x, int y) {
  x += y;
  if (x >= mod) x -= mod;
  return x;
}
inline int sub(int x, int y) {
  x -= y;
  if (x < 0) x += mod;
  return x;
}
inline int mul(int x, int y) { return (x * 1ll * y) % mod; }
inline int powr(int a, long long b) {
  int x = 1 % mod;
  while (b) {
    if (b & 1) x = mul(x, a);
    a = mul(a, a);
    b >>= 1;
  }
  return x;
}
inline int inv(int a) { return powr(a, mod - 2); }
int main() {
  ios_base::sync_with_stdio(false);
  cin.tie(NULL);
  long long a, b;
  cin >> a >> b;
  vector<long long> divisorsA, divisorsB;
  for (long long i = 1; i * i <= a; i++) {
    if (a % i == 0) {
      divisorsA.push_back(i);
      if (i != a / i) divisorsA.push_back(a / i);
    }
  }
  int mA = divisorsA.size();
  for (long long i = 1; i * i <= b; i++) {
    if (b % i == 0) {
      divisorsB.push_back(i);
      if (i != b / i) divisorsB.push_back(b / i);
    }
  }
  int mB = divisorsB.size();
  sort(divisorsA.begin(), divisorsA.end());
  sort(divisorsB.begin(), divisorsB.end());
  long long n = a + b;
  long long ans = 2 * (a + b + 1);
  for (long long i = 1; i * i <= n; i++) {
    if (n % i == 0) {
      long long x = i, y = n / i;
      int ind = lower_bound(divisorsA.begin(), divisorsA.end(), x) -
                divisorsA.begin();
      if (ind == mA) ind--;
      if (divisorsA[ind] > x) ind--;
      long long v = divisorsA[ind];
      if (a / v <= y) {
        ans = min(ans, 2 * (x + y));
      }
      ind = lower_bound(divisorsB.begin(), divisorsB.end(), x) -
            divisorsB.begin();
      if (ind == mB) ind--;
      if (divisorsB[ind] > x) ind--;
      v = divisorsB[ind];
      if (b / v <= y) {
        ans = min(ans, 2 * (x + y));
      }
    }
  }
  cout << ans << '\n';
  return 0;
}