import sys
input=lambda:sys.stdin.readline().rstrip()
def gcd(a,b):
    return gcd(a%b,b%a) if min(a,b,abs(a-b)) else max(a,b)
def isSquare(a):
	m,M=0,a+1
	while M-m>1:
		mid=(M+m)//2
		if mid**2==a:
			return 1
		elif mid**2<a:
			m=mid
		else:
			M=mid
	return 0
def solve():
	N=int(input())
	A=list(map(int,input().split()))
	for i in range(N):
		for j in range(i):
			temp=gcd(A[i],A[j])
			A[i]//=temp
			A[j]//=temp
	return ["No","Yes"][all([isSquare(i) for i in A])]
T=int(input())
for i in range(T):
	print(solve())