// #include // #include // #include // #include // #include // using namespace atcoder; #include using namespace std; constexpr long long INF_LL = 2000000000000000000LL; constexpr int INF = 1000000000; constexpr long long MOD = 998244353; #define all(x) x.begin(), x.end() #define REP(i, a, b) for(int i = a; i < b; i++) #define rep(i, n) REP(i, 0, n) typedef long long ll; typedef pair P; typedef vector vi; typedef vector vvi; typedef vector

vp; typedef vector vl; int dx[4] = {0, -1, 0, 1}; int dy[4] = {1, 0, -1, 0}; int sign[2] = {1, -1}; template bool chmax(T &a, T b) { if(a < b) { a = b; return 1; } return 0; } template bool chmin(T &a, T b) { if(a > b) { a = b; return 1; } return 0; } ll modpow(ll a, ll b, ll m) { if(b == 0) return 1; ll t = modpow(a, b / 2, m); if(b & 1) { return (t * t % m) * a % m; } else { return t * t % m; } } template T gcd(T m, T n) { if(n == 0) return m; return gcd(n, m % n); } struct edge { int to; ll cost; edge(int t, ll c) { to = t, cost = c; } }; typedef vector> graph; // using mint = modint998244353; // using mint = modint1000000007; // constexpr int MAX_COM = 1000001; // mint fac[MAX_COM], ifac[MAX_COM]; // void initfac() { // fac[0] = ifac[0] = 1; // REP(i, 1, MAX_COM) fac[i] = i * fac[i - 1]; // REP(i, 1, MAX_COM) ifac[i] = 1 / fac[i]; // } // mint nCr(int n, int r){ // if(r < 0 || n < r) return 0; // return fac[n] * ifac[n - r] * ifac[r]; // } // typedef ll S; // S op(S x, S y){ return max(x, y); } // S e(){ return 0; } // typedef int F; // S mapping(F f, S x){ // return f + x; // } // F composition(F f, F g){ // return f + g; // } // F id(){ return 0; } void solve(){ int n; cin >> n; vl a(n); rep(i, n) cin >> a[i]; ll s = 0, m = 0; rep(i, n){ s += a[i]; chmax(m, a[i]); } if(s % 3 == 0 && s / 3 >= m){ cout << "Yes" << endl; }else{ cout << "No" << endl; } } int main(){ cin.tie(0); std::ios_base::sync_with_stdio(false); std::cout << std::fixed << std::setprecision(16); // initfac(); int t; t = 1; cin >> t; while (t--) { solve(); } }