#include<bits/stdc++.h>
#define int long long
#define MOD 1000000007
using namespace std;
int T,n,l,r,s[5003],ans,lft,rgt;
int getmod(int X,int Y){return ((X%Y)+Y)%Y;}
void add(int &x,int y){x=(x+y)%MOD;return;}
int getsum(int X){
	X%=MOD;
	return (X*(X+1)/2)%MOD;
}
int getsum2(int X){
	X%=MOD;
	return (X*(X+1)%MOD)*((2*X+1)*166666668ll%MOD)%MOD;
}
int getsum3(int X){
	X%=MOD;
	return (X*(X+1)/2%MOD)*(X*(X+1)/2%MOD)%MOD;
}
int calc(int X,int L,int R,int A,int B){
	int ret=0,v0=0,v1=0,v2=0,v3=0;
	X+=L;
	A-=L;
	B-=L;
	X%=MOD;A%=MOD;B%=MOD;
	v0=-(X*(A*B%MOD)%MOD);
	v1=n*((B*X+A*X-A*B)%MOD)%MOD;
	v2=(n*n%MOD)*(A+B-X)%MOD;
	v3=-(n*n%MOD)*n%MOD;
	add(ret,v0*(((R-L+n)/n)%MOD));
	add(ret,v1*getsum((R-L)/n));
	add(ret,v2*getsum2((R-L)/n));
	add(ret,v3*getsum3((R-L)/n));
	return ret;
}
void sol(){
	cin>>n>>l>>r;
	for(int i=0;i<n;i++){
		cin>>s[i];
		s[i]-=i;
	}
	ans=0;
	for(int i=0;i<n;i++){
		for(int j=0;j<n;j++){
			if(i==j)continue;
			if(getmod(s[i]-s[j],n)!=getmod(j-i,n))continue;
			if(s[i]<s[j])continue; 
			lft=(l/n)*n+i;
			rgt=(r/n)*n+j;
			while(lft<l)lft+=n;
			while(rgt>r)rgt-=n;
			rgt-=(s[i]-s[j]);
			if(lft<=rgt)ans=(ans+2*calc(s[i],lft,rgt,(l-1),(r+1)-(s[i]-s[j])))%MOD;
		}
	}
	for(int i=0;i<n;i++){
		lft=(l/n)*n+i;
		rgt=(r/n)*n+i;
		while(lft<l)lft+=n;
		while(rgt>r)rgt-=n;
		if(lft<=rgt)ans=(ans+calc(s[i],lft,rgt,l-1,r+1))%MOD;
	}
	ans%=MOD;
	ans+=MOD;
	ans%=MOD;
	cout<<ans<<'\n';
	return;
}
signed main(){
	ios::sync_with_stdio(false);
	cin>>T;
	while(T--)sol();
	return 0;
}