#include <bits/stdc++.h>
using namespace std;
#include <atcoder/all>
using namespace atcoder;
#define rep(i, n) for (ll i = 0; i < (ll)(n); i++)
#define rep1(i, n) for (int i = 1; i < (int)(n+1); i++)
#define repn(i,a,b) for (int i = (int)(a) ; i < (int)(b+1); i++)
#define repv(i, n) for (int i = (int)(n-1); i >= 0; i--)
#define all(x) (x).begin(),(x).end()
#define rall(x) (x).rbegin(),(x).rend()
#define sz(x) ((int)(x).size())
#define cauto const auto&
#define bit(n) (1LL<<(n))
#define uni(v) v.erase( unique(v.begin(), v.end()), v.end() );
using ll = long long;
using P = pair<ll,ll>;
using PP = pair<int, P>;
using Graph = vector<vector<int>>;
using mint = modint998244353;
#ifdef LOCAL
    #include <debug.hpp>
    #define debug(...) debug::multi_print(#__VA_ARGS__, __VA_ARGS__)
#else
    #define debug(...) (static_cast<void>(0))
#endif

int main(){
    cin.tie(0);
    ios::sync_with_stdio(false);
    int quary;
    cin >> quary;
    rep(_,quary){
        ll n;
        cin >> n;
        vector<ll> a(n);
        rep(i,n) cin >> a[i];
        rep(i,n){
            repn(j,i+1,n-1){
                ll gcd = __gcd(a[i],a[j]);
                if(gcd != 1){
                    a[i] /= gcd;
                    a[j] /= gcd;
                }
            }
        }
        bool ok = true;
        rep(i,n){
            ll s = static_cast<ll> (sqrt(a[i]) + 0.5);
            if(s*s != a[i]){
                ok = false;
                break;
            }
        }
        cout << (ok ? "Yes" : "No") << "\n";
    }
    return 0;
}