#include<bits/stdc++.h>
using namespace std;
typedef long long ll;
typedef long double db;
const int N=33333;const ll INF=1e18;
int T,n,k,f[N],cnt[N][4];vector<int> state;
struct node{db a,b,c,d;};
node operator *(const node&a,const db&b){
    return {a.a*b,a.b*b,a.c*b,a.d*b};
}
node operator /(const node&a,const db&b){
    return {a.a/b,a.b/b,a.c/b,a.d/b};
}
node operator +(const node&a,const node&b){
    return {a.a+b.a,a.b+b.b,a.c+b.c,a.d+b.d};
}
node operator -(const node&a,const node&b){
    return {a.a-b.a,a.b-b.b,a.c-b.c,a.d-b.d};
}
db v[10][10];node b[10],t[N][5];int tag[N];
void gauss(int st){
    tag[st]=1;
    for(int i=1;i<=4;i++){
        int p;for(p=i;p<=4;p++)if(abs(v[p][i])>1e-6)break;
        if(p==5){tag[st]=0;return;}
        swap(v[i][1],v[p][1]);swap(v[i][2],v[p][2]);
        swap(v[i][3],v[p][3]);swap(v[i][4],v[p][4]);swap(b[i],b[p]);
        for(int j=i+1;j<=4;j++){
            db d=v[j][i]/v[i][i];if(abs(d)<1e-6)continue;
            v[j][1]-=v[i][1]*d;v[j][2]-=v[i][2]*d;v[j][3]-=v[i][3]*d;
            v[j][4]-=v[i][4]*d;b[j]=b[j]-(b[i]*d);
        }
    }
    for(int i=4;i>=1;i--){
        t[st][i]=b[i];
        for(int j=4;j>i;j--)
            t[st][i]=t[st][i]-(t[st][j]*v[i][j]);
        t[st][i]=t[st][i]/v[i][i];
    }
}
void init(){
    n=15;k=4;
    for(int i=1;i<(1<<n);i++)
        for(int j=1;j<=n;j++)if((i>>j-1)&1)
            for(int c=0;c<k;c++)cnt[i][c]+=((j>>c)&1);
    for(int i=0;i<(1<<n);i++)f[i]=1;
    for(int i=1;i<(1<<n);i++)
        for(int j=i;j;j=(j-1)&i){
            int flag=1;
            for(int c=0;c<k;c++)
                if(cnt[j][c]*2!=cnt[i][c])flag=0;
            if(flag)f[i]=0;
        }
    for(int i=1;i<(1<<n);i++)
        for(int j=0;j<n;j++)if((i>>j)&1)f[i]&=f[i^(1<<j)];

    for(int i=1;i<(1<<n);i++)if(f[i]){
        int cnt=0;for(int j=1;j<=n;j++)if((i>>j-1)&1)++cnt;
        if(cnt!=4)continue;
        for(int j=1,m=0;j<=n;j++)if((i>>j-1)&1){
            ++m;for(int l=1;l<=k;l++)v[l][m]=((j>>l-1)&1);
        }
        b[1]={1,0,0,0};b[2]={0,1,0,0};b[3]={0,0,1,0};b[4]={0,0,0,1};
        gauss(i);
    }
    for(int i=1;i<(1<<n);i++)if(f[i])
        for(int j=i;j;j=(j-1)&i)if(f[j]&&j!=i)f[j]=0;
    for(int i=1;i<(1<<n);i++)if(f[i])state.push_back(i);
}
ll a[5],c[20],ans;int x[10],y[10],m;db w[10];
ll calc(){
    int st=(1<<x[1]-1)^(1<<x[2]-1)^(1<<x[3]-1)^(1<<x[4]-1);
    if(!tag[st])return -INF;ll sum=0;
    y[1]=x[1];y[2]=x[2];y[3]=x[3];y[4]=x[4];sort(y+1,y+5);
    for(int i=1;i<=4;i++){
        db r=t[st][i].a*w[1]+t[st][i].b*w[2]+
             t[st][i].c*w[3]+t[st][i].d*w[4];
        ll z=floor(r+0.5);sum+=z*c[y[i]];
        if(z<0||abs(r-z)>1e-6)return -INF;
    }
    return sum;
}
void solve(int st){
    m=0;for(int i=1;i<=n;i++)if((st>>i-1)&1)x[++m]=i;
    if(m==4){w[1]=a[1];w[2]=a[2];w[3]=a[3];w[4]=a[4];ans=max(ans,calc());}
    else{
        for(int i=1;i<=m;i++){
            swap(x[i],x[m]);
            for(int z=0;z<=3;z++){
                w[1]=a[1]-(x[m]&1)*z;
                w[2]=a[2]-((x[m]>>1)&1)*z;
                w[3]=a[3]-((x[m]>>2)&1)*z;
                w[4]=a[4]-((x[m]>>3)&1)*z;
                ans=max(ans,calc()+c[x[m]]*z);
            }
            swap(x[i],x[m]);
        }
    }
}
int main(){
    // freopen("schrodingersbomb.in","r",stdin);
    // freopen("schrodingersbomb.out","w",stdout);
    init();scanf("%d",&T);
    while(T--){
        scanf("%d",&k);n=(1<<k)-1;ans=0;
        for(int i=1;i<=n;i++)scanf("%lld",c+i);while(n<15)++n,c[n]=0;
        for(int i=1;i<=k;i++)scanf("%lld",a+i);while(k<4)++k,a[k]=0;
        for(int x:state)solve(x);printf("%lld\n",ans);
    }
    return 0;
}