#include<bits/stdc++.h>
using namespace std;
const int b=RAND_MAX;
struct jgt{long long a,b;}a[3005];
bool cmp(jgt x,jgt y){return x.a-x.b<y.a-y.b;}
int ran(int l,int r){return rand()%(r-l+1)+1;}
__int128 ans,suma,sumb;
int c[3005],n;
void mnth(int T)
{
    random_device rd;
    mt19937 rnd(rd());
    uniform_int_distribution<int>s(1,n);
    for(int i=1;i<=T;i++)
    {
        int x=s(rnd);
        if(c[x]) c[x]=0,suma-=a[x].a,sumb+=a[x].b;
        else c[x]=1,suma+=a[x].a,sumb-=a[x].b;
        __int128 now=(__int128)suma*sumb;
        if(now>ans) ans=now;
        else
        {
            double tmp=0.001;
            if(rand()>1.0*tmp*b)
            {
                if(c[x])
                c[x]=0,suma-=a[x].a,sumb+=a[x].b;
                else c[x]=1,suma+=a[x].a,sumb-=a[x].b;
            }
        }
    }
}
void print(__int128 x)
{if(!x) return;print(x/10);putchar('0'+(int)x%10);}
int main()
{
    //freopen("1.in","r",stdin);
    //freopen("1.out","w",stdout);
    int tim,A,B;cin>>tim;
    cin>>n>>A>>B;int T=2e7/tim;
    while(tim--)
    {
        long long sb=0;int id=0;
        for(int i=1;i<=n;i++) scanf("%lld",&a[i].a);
        for(int i=1;i<=n;i++) scanf("%lld",&a[i].b);
        suma=0,sumb=0;ans=0;
        for(int i=1;i<=n;i++)
        {
            if(a[i].a>a[i].b)
            suma+=a[i].a,c[i]=1;
            else sumb+=a[i].b,c[i]=0;
        }
        ans=(__int128)suma*sumb;
        mnth(T);print(ans);puts("");
    }
    return 0;
}