#include <bits/stdc++.h>
using namespace std;
using ll = long long;
using P = pair<ll,ll>;
#define fix(x) fixed << setprecision(x)
#define asc(x) x, vector<x>, greater<x>
#define rep(i, n) for(ll i = 0; i < n; i++)
#define all(x) (x).begin(),(x).end()
template<class T>bool chmin(T&a, const T&b){if(a>b){a=b;return 1;}return 0;}
template<class T>bool chmax(T&a, const T&b){if(a<b){a=b;return 1;}return 0;}

int main(){
    cin.tie(nullptr);
    ios::sync_with_stdio(false);
    int t;
    cin >> t;
    while(t--){
        int n;
        cin >> n;
        vector<ll> a(n);
        rep(i,n) cin >> a[i];
        sort(all(a));
        bool ans = true;
        rep(i,n){
            for(int j=i+1;j<n;j++){
                ll g = __gcd(a[i], a[j]);
                a[i] /= g; a[j] /= g;
            }
            ll ok = 1, ng = 1e9+1;
            while(ng-ok>1){
                ll md = (ng+ok)/2;
                if(md*md<=a[i]) ok = md;
                else ng = md;
            }
            if(ok*ok!=a[i]){
                ans = false;
                break;
            }
        }
        cout << (ans?"Yes\n":"No\n");
    }
    return 0;
}