    #include <bits/stdc++.h>
#include<map>
#define all(x) x.begin(),x.end()
#define rall(x) x.rbegin(),x.rend()
#define pb(x) emplace_back(x)
#define forn(i,N) for(int i=0;i<N;i++)
#define rep(i,N) for(int i=0;i<=N;i++)

#include<vector>
#include<set>
#define ll long long int
#define S second

using namespace std;
#define IOS ios::sync_with_stdio(0); cin.tie(0); cout.tie(0);
using ld=long double;
using vi=vector<int>;
using vll=vector<ll>;
using vd=vector<double>;
using vvi=vector<vi>;
int ctoi(char c) {return c-'0';}
int slen(string s) {return s.length();}
int vsize(vi a) {return (int)a.size();}
int in() {int x;cin>>x;return x;}
#define fast ios_base::sync_with_stdio(false);cin.tie(NULL);
#define ld long double
// I failed countless times
// I am still failing,
// I don't know how to win but, 
// I know how to not give up ✌

typedef pair<ll,ll> pll;
typedef vector<pll> vpll;
typedef map<ll,ll> mpll;
typedef set<ll> stll;


void kmin(ll n, ll m){
    ll val = n / m; 
    n = n - (val*m);
    ll val1 = (m-n)*(val-1)*(val)/2;
    ll val2 = n*val*(val+1)/2;
    ll val3 = val1 + val2;
    cout<<val3<<" ";
}
void kmax(ll n, ll m){
    n = n -(m-1);
    ll val = n*(n-1)/2;
    cout<<val;
    cout<<endl;
}
const int M = 1000000  ;

 
ll prime(ll x){
    for(ll i=2;i*i<=x;i++){
        if(x%i==0)return i;
    }
    return x;
}
bool cmp(pair<ll, pair<ll,ll> > x, pair<ll,pair<ll,ll> >y){
    return x>y;
}

vector<pair<ll, ll> > v;
set<ll> s;
ll n;
bool check(ll x) {
	for (int i = 1; i < n; i++) {
		if (v[i].first % x != 0 && v[i].second % x != 0) {
			return 0;
		}
	}
	return 1;
}
void get_prime(ll x) {
	if (x % 2 == 0)
		s.insert(2);
	while (x % 2 == 0)x = x / 2;
	for (int i = 3; i <= sqrt(x); i += 2) {
		if (x % i == 0) {
			s.insert(i);
			while (x % i == 0)
				x = x / i;
		}
	}
	if (x > 2) {
		s.insert(x);
	}
}


vll ar[1000000];
stll cycle;
vll path;
bool vis[1000000];
 
void dfs(ll node,ll par){
	if(vis[node]==true){
		if(cycle.size()) return;
 
		ll i=path.size()-1;
 
		while(i>=0 && path[i]!=node){
			cycle.insert(path[i]);
			i--;
		}
		cycle.insert(node);
		return;
	}
 
	vis[node]=true;
	path.pb(node);
	for(auto x:ar[node]){
		if(x!=par){
			dfs(x,node);
		}
	}
	vis[node]=false;
	path.pop_back();
}
const int mod=1e9+7; 
ll dfs1(ll node,ll par){
	ll cnt=1;
 
	for(auto x:ar[node]){
		if(x!=par){
			cnt+=dfs1(x,node);
		}
	}
 
	return cnt;
}
ll power(ll a,ll b)
{
    ll res=1LL;
    a=a%mod;
    while(b>0)
    {
        if(b&1) res=(res*a)%mod;
        b=b/2;
        a=(a*a)%mod;
    }
    return res;
} 
 

 
 ll a[1000000];
int isValley(int i) {
    return (i > 0 && i < n - 1 && a[i] < a[i - 1] && a[i] < a[i + 1]);
}
 
int isHill(int i) {
    return (i > 0 && i < n - 1 && a[i] > a[i - 1] && a[i] > a[i + 1]);
}
 
void solve() {
    ll n;
    cin >> n;
    map<ll, ll> mp;
    for(int i = 0;i < 2 * n;++i) {
        ll x; 
        cin >> x;
        mp[x]++;
    }
    vector<ll> d;
    for(auto i : mp) {
        if(i.second % 2) {
            cout << "NO";
            return;
        }
        for(int j = 0;j < i.second / 2;++j) {
            d.push_back(i.first);
        }
    }
    sort(d.rbegin(),d.rend());
    ll y = 0;
    ll prev=0;
    int nc = n;
    for(int i = 0;i < nc;++i) {
        if((d[i] - 2LL * y) % (2 * n) != 0){
            cout << "NO";
            return;
        }
        long long x = (d[i] - 2LL * y) / (2 * n);
        if(x <= 0 || prev==x) {
            cout << "NO";
            return;
        }
        prev=x;
        y += x;
        n--;
    }
    cout << "YES";
}


int main()
{ fast;
    ll t=1;
   cin>>t;
    
    while(t--){
      
     
      solve();
      cout<<"\n";
     
    }

    return 0;
}



