//24 -> 2^3 * 3 -> 2 * 3のように置き換えてOK。 //各素因数の個数 mod 2の累積和を考えると、これをキーにして数えたいが、2 3 5 7 11 …のようなケースで膨らむ。 //集合をハッシュで管理できないか? -> ゾブリストハッシュ! //今回のように、2回同じ要素が来ると元に戻る性質のときに使えるハッシュ。 //unsigned long longで乱数生成すれば、衝突の心配も少ない。 #include #include #include #define rep(i, n) for(i = 0; i < n; i++) using namespace std; const int MAX = 1000000; int n; int a[200000]; bool isPrime[MAX + 1]; vector ps[MAX + 1]; unsigned long long Rnd[MAX + 1]; mt19937_64 mt(2521); void init() { for (int i = 2; i <= MAX; i++) isPrime[i] = true; for (int i = 2; i <= MAX; i++) { if (isPrime[i]) { for (int j = i * 2; j <= MAX; j += i) { isPrime[j] = false; int tmp = j; int cnt = 0; while (tmp % i == 0) { tmp /= i; cnt++; } if (cnt % 2 == 1) { ps[j].push_back(i); } } ps[i].push_back(i); } } for (int i = 0; i <= MAX; i++) { Rnd[i] = mt(); } } map cnts; signed main() { int i, j; init(); scanf("%d", &n); rep(i, n) scanf("%d", a + i); unsigned long long zhash = 0; cnts[zhash]++; rep(i, n) { rep(j, ps[a[i]].size()) { int p = ps[a[i]][j]; zhash ^= Rnd[p]; } cnts[zhash]++; } long long ans = 0; for (map::iterator it = cnts.begin(); it != cnts.end(); it++) { int cnt = it->second; ans += (long long)cnt * (cnt - 1) / 2; } printf("%lld\n", ans); return 0; }