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

const int N=155,H=1e9+7;
int n,l[N],r[N],tot=1,f[N][3],a[N<<2],len,ans,fac[N],ifac[N];

int adc(int a,int b){return a+b>=H?a+b-H:a+b;}
int dec(int a,int b){return a<b?a-b+H:a-b;}
int mul(int a,int b){return 1ll*a*b%H;}
void add(int &a,int b){a=adc(a,b);}
void del(int &a,int b){a=dec(a,b);}

int qpow(int a,int b=H-2){
  int res=1;
  while(b){if(b&1) res=mul(res,a);a=mul(a,a),b>>=1;}
  return res;
}

void init(){
  fac[0]=1;
  for(int i=1;i<N;i++) fac[i]=mul(fac[i-1],i);
  ifac[N-1]=qpow(fac[N-1]);
  for(int i=N-1;i>=1;i--) ifac[i-1]=mul(ifac[i],i);
}

int Len(int l,int r){return max(0,r-l+1);}

int F(int x){
  f[0][0]=1;
  for(int i=1;i<=n;i++){
    f[i][0]=mul(Len(l[i],min(r[i],x-1)),adc(f[i-1][0],f[i-1][2]));
    f[i][1]=mul(Len(max(l[i],x),r[i]),f[i-1][0]);
    f[i][2]=mul(Len(l[i],min(r[i],x-1)),f[i-1][1]);
  }
  int res=tot;
  for(auto i:{0,1,2}) del(res,f[n][i]);
  return res;
}

int sgn(int x){return x&1?H-1:1;}

int Lag(vector<int> &X,vector<int> &Y,int pos){
  int res=0,m=(int)X.size(),s=1;
  for(int i=0;i<m;i++) s=mul(s,pos-X[i]);
  for(int i=0;i<m;i++) add(res,mul(mul(s,qpow(pos-X[i])),mul(mul(sgn(i),Y[i]),mul(ifac[i],fac[m-i-1]))));
  return res;
}

void SOLVE(){
  cin>>n,tot=1,len=ans=0;
  for(int i=1;i<=n;i++) cin>>l[i],a[++len]=l[i]-1;
  for(int i=1;i<=n;i++) cin>>r[i],tot=mul(tot,Len(l[i],r[i])),a[++len]=r[i];
  sort(a+1,a+len+1),len=unique(a+1,a+len+1)-a-1;
  
  for(int i=1;i<len;i++){
    int L=a[i]+1,R=a[i+1];
    if(Len(L,R)<=n+2) for(int x=L;x<=R;x++) add(ans,F(x));
    else{
      vector<int> X(0),Y(0);
      for(int x=L,s=0;x<=L+n+1;x++) add(s,F(x)),X.pb(x),Y.pb(s);
      add(ans,Lag(X,Y,R));
    }
  }
  cout<<adc(ans,mul(a[1],tot))<<'\n';
}

int main(){
  ios::sync_with_stdio(false);cin.tie(0);cout.tie(0);
  int _;cin>>_,init();
  while(_--) SOLVE();
  return 0;
}