#include <bits/stdc++.h>
using namespace std;
const int maxn = 1e5 + 7;
int prim[maxn], vis[maxn], cnt;
void sprim() {
  for (int i = 2; i < maxn; i++) {
    if (!vis[i]) prim[++cnt] = i;
    for (int j = 1; j <= cnt && i * prim[j] < maxn; j++) {
      vis[i * prim[j]] = 1;
      if (i % prim[j] == 0) break;
    }
  }
}
struct Node {
  vector<pair<int, int>> vec;
  bool operator<(const Node &T) const {
    int i = 0, j = 0;
    while (i < vec.size() && j < T.vec.size()) {
      if (vec[i] != T.vec[j]) return vec[i] < T.vec[j];
      i++, j++;
    }
    return (i == vec.size()) && (j != T.vec.size());
  }
  bool operator==(const Node &T) const {
    if (vec.size() != T.vec.size()) return 0;
    for (int i = 0; i < vec.size(); i++) {
      if (vec[i] != T.vec[i]) return 0;
    }
    return 1;
  }
};
map<Node, int> Last;
int Preprim[maxn], n, x, k;
long long ans;
int main() {
  sprim();
  scanf("%d%d", &n, &k);
  for (int i = 1; i <= n; i++) {
    scanf("%d", &x);
    Node Fid, Add;
    Fid.vec.clear();
    Add.vec.clear();
    for (int j = 1; x != 1; j++) {
      if (x % prim[j] == 0) {
        int cnt = 0;
        while (x % prim[j] == 0) {
          x /= prim[j];
          cnt++;
        }
        if (cnt % k == 0) continue;
        Add.vec.push_back({prim[j], cnt % k});
        Fid.vec.push_back({prim[j], (k - (cnt % k)) % k});
      }
    }
    ans += Last[Fid];
    Last[Add]++;
  }
  printf("%lld\n", ans);
  return 0;
}
