from collections import defaultdict
from math import isqrt,gcd
#==================================================
def is_square(x):
    y=isqrt(x)
    return y*y==x

def solve():
    N=int(input())
    A=list(map(int,input().split()))

    for i in range(N):
        if is_square(A[i]):
            A[i]=1

    for i in range(N):
        for j in range(i+1,N):
            g=gcd(A[i], A[j])
            A[i]//=g
            A[j]//=g

    return all([is_square(a) for a in A])

#==================================================
T=int(input())
for _ in range(T):
    print("Yes" if solve() else "No")