#include<bits/stdc++.h>
using namespace std;

typedef long long ll;
typedef unsigned long long ull;

const int N=5e5+9,M=20009,mod=998244353;
const ll inf=1e16;
int n,m,a[N],b[N];
int K,c[N];
struct Info{
	ll ans;int cnt;
	Info(ll x=inf,int y=0){ans=x,cnt=y;}
};
Info operator+ (Info x,Info y)
{
	Info ans;
	ans.ans=min(x.ans,y.ans);
	ans.cnt=0;
	if(ans.ans==x.ans) ans.cnt=x.cnt;
	if(ans.ans==y.ans) (ans.cnt+=y.cnt)%=mod;
	return ans;
}
Info operator+ (Info x,int y){x.ans+=y;return x;}
Info dp[M][200];
vector<int> vec[M];
void work()
{
	scanf("%d%d",&n,&m);
	for(int i=1;i<=n;i++) scanf("%d",&a[i]);
	for(int i=1;i<=n;i++) scanf("%d",&b[i]);
	int mx=0,mxa=0;
	for(int i=1;i<=m;i++) mx=max(mx,a[i]+b[i]),mxa=max(mxa,a[i]);
	for(int i=1;i<=n;i++) vec[a[i]].push_back(i);
	scanf("%d",&K);
	for(int i=1;i<=K;i++) scanf("%d",&c[i]);
	for(int j=0;j<=mx+K+K;j++) for(int i=0;i<=K;i++)  dp[j][i]=Info();
	
	dp[0][0]=Info(0,1);
	Info ans=Info();
	for(int i=0;i<=mx+K+K;i++)
	for(int j=0;j<=K;j++)
	if(dp[i][j].ans!=inf)
	{
//		cerr<<i<<" "<<j<<" "<<dp[i][j].ans<<endl;
		if(i>=mx+j) {ans=ans+dp[i][j];continue;}
		for(int k=1;k<=K;k++)
		{
			int tmp=j;
			for(int l=i+1;l<=i+k;l++)
			for(int p:vec[l]) tmp=max(tmp,i+k+b[p]-mx);
			dp[i+k][tmp]=dp[i+k][tmp]+(dp[i][j]+c[k]);
		}
	}
	printf("%lld %d\n",ans.ans,ans.cnt);
	for(int i=1;i<=mxa;i++) vec[i].clear();
}
signed main()
{
	int T;
	scanf("%d",&T);
	while(T--) work();
}