#include <cstdio>
#include <cstdint>
#include <cassert>
#include <cctype>
#include <algorithm>
#include <string>
#include <vector>

#define fprintf(...) void(0)

intmax_t gcd(intmax_t m, intmax_t n) {
  while (n) std::swap(m%=n, n);
  return m;
}

class Polynomial {
  std::vector<intmax_t> coef;  // ascending order

public:
  Polynomial(intmax_t n=1): coef{n} {}  // constant
  Polynomial(const std::vector<intmax_t> &coef): coef(coef) {}

  Polynomial &operator +=(const Polynomial &oth) {
    if (coef.size() < oth.coef.size())
      coef.resize(oth.coef.size());

    for (size_t i=0; i<oth.coef.size(); ++i)
      coef[i] += oth[i];

    return *this;
  }

  Polynomial &operator -=(const Polynomial &oth) {
    if (coef.size() < oth.coef.size())
      coef.resize(oth.coef.size());

    for (size_t i=0; i<oth.coef.size(); ++i)
      coef[i] -= oth[i];

    return *this;
  }

  Polynomial &operator *=(const Polynomial &oth) {
    std::vector<intmax_t> res(coef.size()+oth.coef.size()-1);
    for (size_t i=0; i<coef.size(); ++i)
      for (size_t j=0; j<oth.coef.size(); ++j)
        res[i+j] += coef[i]*oth[j];

    coef.swap(res);
    return *this;
  }

  Polynomial &operator %=(Polynomial oth) {
    if (coef.size() < oth.coef.size()) return *this;

    for (size_t i=coef.size(); i-->=oth.coef.size();) {
      intmax_t a=coef[i], b=oth.coef.back();
      assert(b != 0);
      if (a == 0) continue;
      intmax_t g=gcd(a, b);
      intmax_t li=b/g, lj=a/g;

      if (li != 1)
        for (intmax_t &ai: coef) ai *= li;
      if (lj != 1)
        for (intmax_t &aj: oth.coef) aj *= lj;

      for (size_t jl=i, jr=oth.coef.size(); jr--; --jl) {
        coef[jl] -= oth[jr];
      }
    }

    coef.resize(oth.coef.size()-1);
    normalize();
    return *this;
  }

  Polynomial &normalize() {
    while (coef.size() > 1 && coef.back() == 0)
      coef.pop_back();

    intmax_t g=abs(coef[0]);
    for (size_t i=1; i<coef.size(); ++i)
      g = gcd(abs(g), abs(coef[i]));

    if (coef.back() < 0) g = -g;
    if (g != 0)
      for (size_t i=0; i<coef.size(); ++i)
        coef[i] /= g;

    return *this;
  }

  Polynomial &operator ^=(int iexp) {
    Polynomial res({1});
    for (Polynomial dbl(coef); iexp; iexp>>=1) {
      if (iexp & 1) res *= dbl;
      dbl *= dbl;
    }
    coef.swap(res.coef);
    return *this;
  }

  Polynomial operator +() const {
    return *this;
  }

  Polynomial operator -() const {
    std::vector<intmax_t> neg=coef;
    for (intmax_t &n: neg) n = -n;
    return Polynomial(neg);
  }

  intmax_t &operator [](intmax_t n) {
    return coef[n];
  }

  intmax_t operator [](intmax_t n) const {
    return coef[n];
  }

  intmax_t operator ()(intmax_t x) const {
    intmax_t res=0;
    for (size_t i=coef.size(); i--;)
      res = res*x + coef[i];

    return res;
  }

  operator bool() const {
    return coef.size() > 1 || coef[0] != 0;
  }

  void swap(Polynomial &oth) {
    coef.swap(oth.coef);
  }

  std::string to_str() const {
    std::string res;
    for (size_t i=coef.size(); i--;) {
      intmax_t n=coef[i];
      if (n == 0) continue;
      if (n > 0) {
        if (!res.empty()) res += '+';
      } else {
        res += '-';
        n = -n;
      }
      if (n != 1 || i == 0)
        res += std::to_string(n);

      if (i > 0) {
        res += 'x';
        if (i > 1) res += "^" + std::to_string(i);
      }
    }
    if (res.empty()) res = "0";
    return res;
  }
};

std::string expl_times(const std::string &s) {
  std::string res(1, s[0]);
  for (size_t i=1; i<s.length(); ++i) {
    if (s[i] == 'x' || s[i] == '(')
      if (isdigit(s[i-1]) || s[i-1] == 'x' || s[i-1] == ')')
        res += '*';

    res += s[i];
  }
  return res;
}

Polynomial parse(const std::string &s, size_t &i, size_t preced=0) {
  static const std::vector<std::string> ops={"+-", "*", "^"};

  if (preced == ops.size()) {
    if (s[i] == '(') {
      Polynomial res=parse(s, ++i);
      assert(s[i] == ')');
      ++i;
      return res;
    }
    if (s[i] == '-') {
      return -parse(s, ++i, preced);
    }
    if (isdigit(s[i])) {
      intmax_t res=s[i]-'0';
      while (++i < s.length() && isdigit(s[i]))
        res = res*10+s[i]-'0';

      return res;
    }
    if (s[i] == 'x') {
      ++i;
      return Polynomial({0, 1});
    }
    assert(false);
  }

  Polynomial res=parse(s, i, preced+1);
  while (i < s.length()) {
    char op=s[i];
    if (!std::count(ops[preced].begin(), ops[preced].end(), op)) break;
    Polynomial tmp=parse(s, ++i, preced+1);

    if (op == '+') res += tmp;
    if (op == '-') res -= tmp;
    if (op == '*') res *= tmp;
    if (op == '^') res ^= tmp[0];
  }
  return res;
}

int testcase_ends() {
  char buf[16384];
  scanf("%s", buf);
  std::string lhs=expl_times(buf);
  if (lhs == ".") return 1;
  scanf("%s", buf);
  std::string rhs=expl_times(buf);

  size_t i=0;
  Polynomial lp=parse(lhs, i);
  size_t j=0;
  Polynomial rp=parse(rhs, j);

  while (rp)
    (lp %= rp).swap(rp);

  lp.normalize();
  printf("%s\n", lp.to_str().c_str());
  return 0;
}

int main() {
  while (!testcase_ends()) {}
}



Measured memory_bytes: 2728