#include<bits/stdc++.h>
using namespace std;
typedef long long ll;
const int N=3e4+10,K=110,mod=998244353;
struct node{
	ll val;
	int cnt;
	node add(int x){return node{val+x,cnt};}
}f[N][K];
void operator +=(node &a,node b)
{
	if(a.cnt==0||a.val>b.val) return a=b,void();
	else if(b.cnt==0||a.val<b.val) return;
	a.cnt+=b.cnt;
	if(a.cnt>mod) a.cnt-=mod;
}
int T,n,m,k,t[N],h[N][K],c[K],a[N*100],b[N*100];
int main()
{
	scanf("%d",&T);
	for(int i=1;i<N;i++) t[i]=-1e9;
	while(T--)
	{
		scanf("%d%d",&n,&m),m=0;
		for(int i=1;i<=n;i++) scanf("%d",&a[i]);
		for(int i=1;i<=n;i++) scanf("%d",&b[i]);
		for(int i=1;i<=n;i++)
		{
			t[b[i]]=max(t[b[i]],a[i]);
			m=max(m,a[i]+b[i]);
		}
		scanf("%d",&k);
		for(int i=1;i<=k;i++) scanf("%d",&c[i]);
		for(int i=0;i<=m+2*k;i++)
			for(int j=0;j<=k;j++) f[i][j]=node{0,0};
		for(int i=0;i<=m+k;i++)
		{
			for(int j=1;j<=k;j++)
				h[i][j]=max(t[i+j]+i,h[i][j-1]);
			for(int j=1;j<=k;j++)
				h[i][j]=max(0,h[i][j]+j-m);
		}
		f[0][0]=node{0,1};
		for(int i=0;i<=m+k;i++)
			for(int x=0;x<=k;x++)
				for(int y=1;y<=k;y++)
					f[i+y][max(x,h[i][y])]+=f[i][x].add(c[y]);
		node ans=node{0,0};
		for(int i=0;i<=2*k;i++)
			for(int j=0;j<=min(i,k);j++) ans+=f[i+m][j];
		printf("%lld %d\n",ans.val,ans.cnt%mod);
		for(int i=1;i<=m;i++) t[i]=-1e9;
	}
}