#include <bits/stdc++.h>
#define ll long long
const int N = 1e5, M = 1e5;

int n, m, S, T, c[N];

const ll INFLL = 1e18;
struct flow {
  int cnt = 1, hd[N], nxt[M << 1], to[M << 1]; ll limit[M << 1], cst[M << 1];
  void add(int u, int v, ll w, ll c) {
    nxt[++cnt] = hd[u], hd[u] = cnt, to[cnt] = v, limit[cnt] = w, cst[cnt] = c;
    nxt[++cnt] = hd[v], hd[v] = cnt, to[cnt] = u, limit[cnt] = 0, cst[cnt] = -c;
  }
  int fr[N], in[N]; ll fl[N], dis[N];
  std::pair<ll, ll> mcmf(int n, int s, int t) {
    ll flow = 0, cost = 0;
    while(1) {
      std::queue<int> q;
      std::fill(dis, dis + n + 1, INFLL);
      q.push(s), fl[s] = INFLL, dis[s] = 0;
      while(!q.empty()) {
        int t = q.front();
        q.pop(), in[t] = 0;
        for(int i = hd[t]; i; i = nxt[i]) {
          int it = to[i]; ll d = dis[t] + cst[i];
          if(limit[i] && d < dis[it]) {
            dis[it] = d, fl[it] = std::min(fl[t], limit[i]), fr[it] = i;
            if(!in[it]) in[it] = 1, q.push(it);
          }
        }
      }
      if(dis[t] >= INFLL) return {flow, cost};
      flow += fl[t], cost += dis[t] * fl[t];
      for(int u = t; u != s; u = to[fr[u] ^ 1]) limit[fr[u]] -= fl[t], limit[fr[u] ^ 1] += fl[t];
    }
  }
  void init(int n) {
    cnt = 1; std::fill(hd, hd + n + 1, 0);
  }
};
struct bounded_flow {
  int e, u[M], v[M], lo[M], hi[M], cst[M];
  void add(int _u, int _v, ll l, ll r, ll c) {
    if(c < 0) {
      u[++e] = _u, v[e] = _v, lo[e] = r, hi[e] = r, cst[e] = c;
      u[++e] = _v, v[e] = _u, lo[e] = 0, hi[e] = r - l, cst[e] = -c;
    }
    else u[++e] = _u, v[e] = _v, lo[e] = l, hi[e] = r, cst[e] = c;
  }
  flow g;
  std::pair<ll, ll> mcmf(int n, int s, int t) {
    int ss = 0, tt = n + 1;
    static ll w[N];
    std::fill(w + 1, w + n + 1, 0);
    int flow = 0, cost = 0, tot = 0;
    g.init(n + 1);
    for(int i = 1; i <= e; i++) {
      w[u[i]] -= lo[i], w[v[i]] += lo[i];
      cost += lo[i] * cst[i];
      g.add(u[i], v[i], hi[i] - lo[i], cst[i]);
    }
    for(int i = 1; i <= n; i++)
      if(w[i] > 0) g.add(ss, i, w[i], 0), tot += w[i];
      else if(w[i] < 0) g.add(i, tt, -w[i], 0);
    g.add(t, s, INFLL, 0);
    std::pair<ll, ll> res = g.mcmf(n + 1, ss, tt);
    if (res.first != tot) return {-1, -1};
    cost += res.second;
    flow += g.limit[g.hd[s]];
    g.hd[s] = g.nxt[g.hd[s]], g.hd[t] = g.nxt[g.hd[t]];
    res = g.mcmf(n + 1, s, t);
    return {flow + res.first, cost + res.second};
  }
  std::pair<ll, ll> mcf(int n, int s, int t) {
    int ss = 0, tt = n + 1;
    static ll w[N];
    std::fill(w + 1, w + n + 1, 0);
    int flow = 0, cost = 0, tot = 0;
    g.init(n + 1);
    for(int i = 1; i <= e; i++) {
      w[u[i]] -= lo[i], w[v[i]] += lo[i];
      cost += lo[i] * cst[i];
      g.add(u[i], v[i], hi[i] - lo[i], cst[i]);
    }
    for(int i = 1; i <= n; i++)
      if(w[i] > 0) g.add(ss, i, w[i], 0), tot += w[i];
      else if(w[i] < 0) g.add(i, tt, -w[i], 0);
    g.add(t, s, INFLL, 0);
    std::pair<ll, ll> res = g.mcmf(n + 1, ss, tt);
    if (res.first != tot) return {-1, -1};
    return {res.first, res.second + cost};
  }
  void init() { e = 0; }
} f;

void Solve() {
    scanf("%d %d", &n, &m);
    n += n;
    S = n + n + m + 1;
    T = S + 1;
    f.init();
    for (int i = 1; i <= m; i++) {
        int x;
        scanf("%d", &x);
        f.add(S, i, x, n / 2, 0);
    }
    for (int i = 1; i <= n; i++) {
        scanf("%d", &c[i]);
    }
    for (int i = 1; i <= n; i++) {
        int w;
        scanf("%d", &w);
        f.add(c[i], i + m, 0, 1, -w);
        f.add(i + m, i + n + m, 0, 1, 0);
    }
    for (int i = 1; i < n; i++) {
        f.add(i + n + m, i + 1 + n + m, (i + 1) / 2, i, 0);
    }
    f.add(n + n + m, T, n / 2, n / 2, 0);

    auto res = f.mcf(n + n + m + 2, S, T);
    if (res.first == -1 && res.second == -1) {
        puts("-1");
    } else {
        printf("%lld\n", -res.second);
    }
}

int main() {
    int tt;
    scanf("%d", &tt);
    while (tt--) {
        Solve();
    }

    return 0;
}