#include<bits/stdc++.h>
using namespace std;
typedef long long ll;
typedef long double db;
typedef pair<int,int> pr;
const int mod=998244353;
const db eps=1e-9;
ll b[20],ans;
int n,c[20],N,a[20],to[20],tc[20];
bool bad;
db f[9][9],tf[9][9];
void Getinv(){
	for(int i=0;i<n;i++)f[i][i+n]=1;
	for(int i=0;i<n;i++){
		int p=i;
		while(p<n&&fabs(f[p][i])<eps)p++;
		if(p>=n)return bad=1,void();
		for(int j=i;j<n+n;j++)swap(f[p][j],f[i][j]);
		for(int j=i+1;j<n;j++){
			db v=f[j][i]/f[i][i];
			for(int k=i;k<n+n;k++)f[j][k]-=v*f[i][k];
		}
	}
	for(int i=n-1;i>=0;i--){
		db v=1/f[i][i];
		for(int j=i;j<n+n;j++)f[i][j]*=v;
		for(int j=i-1;j>=0;j--){
			db x=f[j][i];
			for(int k=i;k<n+n;k++){
				f[j][k]-=x*f[i][k];
			}
		}
	}
}
void dfs(int x){
	if(x==N){
		ll sum=0;
		for(int i=0;i<N;i++)sum+=1ll*a[i]*b[i];
		ans=max(ans,sum);
		return ;
	}
	for(int i=0;;i++){
		a[x]+=i;
		bool ok=1;
		for(int j=0;j<n;j++){
			if((x+1)&(1<<j))tc[j]-=i;
			if(tc[j]<0)ok=0;
		}
		if(!ok){
			for(int j=0;j<n;j++)if((x+1)&(1<<j))tc[j]+=i;
			a[x]-=i;
			return ;
		}
		dfs(x+1);
		for(int j=0;j<n;j++)if((x+1)&(1<<j))tc[j]+=i;
		a[x]-=i;
	}
}
void Solve(){
	cin>>n,N=(1<<n)-1;
	for(int i=0;i<N;i++)cin>>b[i];
	for(int i=0;i<n;i++)cin>>c[i];
	ans=0;
	for(int i=0;i<(1<<N);i++){
		if(__builtin_popcount(i)!=N-n)continue;
		memset(f,0,sizeof(f)),bad=0;
		int id[20]={0},tot=0;
		for(int j=0;j<N;j++){
			if(i&(1<<j))continue;
			to[tot]=j,id[j]=(tot++);
			for(int k=0;k<n;k++){
				f[k][id[j]]=((j+1)>>k)%2;
				tf[k][id[j]]=f[k][id[j]];
			}
		}
		Getinv();
		if(bad)continue;
		memset(a,0,sizeof(a));
		bool nd=0;
		for(int i=0;i<n;i++){
			db v=0;
			for(int j=0;j<n;j++)v+=c[j]*f[i][j+n];
			if(fabsl(v-roundl(v))<eps)a[to[i]]=roundl(v);
			else a[to[i]]=floorl(v),nd=1;
			if(a[to[i]]<0)bad=1;
		}
		if(bad)continue;
		//if(nd)for(int i=0;i<N;i++)a[i]=max(a[i]-1,0);
		for(int i=0;i<n;i++){
			tc[i]=c[i];
			for(int j=0;j<N;j++)if((j+1)&(1<<i))tc[i]-=a[j];
		}
		dfs(0);
	}
	cout<<ans<<'\n';
}
int main(){
	int t;
	cin>>t;
	while(t--)Solve();
}