#include <bits/stdc++.h>
using namespace std;
using ll = long long;
using vll = vector<ll>;
using ull = unsigned long long;
using int2 = pair<ll, ll>;
using int3 = tuple<ll, ll, ll>;
template <typename A>
string to_string(const A &v);
string to_string(const string &s) { return '"' + s + '"'; }
string to_string(const char *s) { return to_string((string)s); }
string to_string(char c) { return string(1, c); }
string to_string(bool b) { return (b ? "true" : "false"); }
template <typename A, typename B>
string to_string(const pair<A, B> &p) {
  return "(" + to_string(p.first) + ", " + to_string(p.second) + ")";
}
template <typename A, typename B, typename C>
string to_string(const tuple<A, B, C> &p) {
  return "(" + to_string(std::get<0>(p)) + ", " + to_string(std::get<1>(p)) +
         ", " + to_string(std::get<2>(p)) + ")";
}
template <typename A>
string to_string(const A &v) {
  bool first = true;
  string res = "{";
  for (const auto &x : v) {
    if (!first) {
      res += ", ";
    }
    first = false;
    res += to_string(x);
  }
  res += "}";
  return res;
}
template <typename T>
void sort(vector<T> &val) {
  sort(val.begin(), val.end());
}
template <typename T>
void reverse(vector<T> &val) {
  reverse(val.begin(), val.end());
}
template <typename T>
void sort_b(vector<T> &val) {
  sort(val.begin(), val.end(), greater<T>());
}
template <typename T>
T pair_sum(pair<T, T> &val) {
  return val.first + val.second;
}
template <typename T>
T min(T a, T b, T c) {
  return min(min(a, b), c);
}
template <typename T>
T max(T a, T b, T c) {
  return max(max(a, b), c);
}
template <typename T>
T mul_mod(T a, T b, T mod) {
  T res = 0;
  a %= mod;
  while (b) {
    if (b & 1) {
      res = (res + a) % mod;
    }
    a = (2 * a) % mod;
    b >>= 1;
  }
  return res;
}
void debug_out() { cerr << endl; }
template <typename Head, typename... Tail>
void debug_out(Head H, Tail... T) {
  cerr << " " << to_string(H);
  debug_out(T...);
}
class mytask {
 public:
  mytask() {}
  int solve() {
    ll t;
    cin >> t;
    while (t--) {
      ll p, q;
      cin >> p >> q;
      if (p % q != 0) {
        cout << p << "\n";
        continue;
      }
      map<ll, ll> del_q;
      for (ll i = 2; i * i <= q; i++) {
        while (q % i == 0) {
          del_q[i]++;
          q /= i;
        }
      }
      if (q > 1) del_q[q]++;
      ;
      ;
      ll res = 1;
      for (auto [k, v] : del_q) {
        ll tmp_p = p;
        ll tmp_d = pow(k, v);
        while (tmp_p % tmp_d == 0) {
          tmp_p /= k;
        }
        res = max(res, tmp_p);
      }
      cout << res << "\n";
    }
    return 0;
  }
};
int main() {
  auto task = mytask();
  task.solve();
  return 0;
}
