#include <bits/stdc++.h>
using namespace std;
const int mod = 1e9 + 7;
const int inf = 0x4f4f4f4f4f;
int n, m, x, y;
const int maxn = 1e6 + 5;
long long b[2 * maxn];
long long poww(long long x, long long k) {
  long long t = 1;
  while (k) {
    if (k % 2) {
      t = (t * x) % mod;
    }
    x = (x * x) % mod;
    k /= 2;
  }
  return t;
}
long long c(long long x, long long y) {
  if (x < y) return 0;
  if (y == 0) return 1;
  long long t = 1;
  long long tmp = (b[x - y] * b[y]) % mod;
  t = (b[x] * poww(tmp, mod - 2)) % mod;
  return t;
}
void init() {
  b[1] = 1;
  b[0] = 1;
  for (int i = 1; i < 2 * maxn; i++) {
    b[i] = (b[i - 1] * i) % mod;
  }
}
int main() {
  std::ios::sync_with_stdio(false);
  cin.tie(0);
  cout.tie(0);
  init();
  cin >> n;
  while (n--) {
    cin >> x >> y;
    long long ans = poww(2, y - 1);
    for (int i = 2; i * i <= x; i++) {
      if (x % i == 0) {
        int k = 0;
        while (x % i == 0) k++, x /= i;
        ans = (ans * c(y + k - 1, y - 1)) % mod;
      }
    }
    if (x > 1) ans = (ans * y) % mod;
    cout << ans << endl;
  }
  return 0;
}
