题目链接:http://acm.hdu.edu.cn/showproblem.php?pid=5145
#include<bits/stdc++.h>
using namespace std;#define debug puts("YES");
#define rep(x,y,z) for(int (x)=(y);(x)<(z);(x)++)
#define read(x,y) scanf("%d%d",&x,&y)
#define ll long long#define lrt int l,int r,int rt
#define lson l,mid,rt<<1
#define rson mid+1,r,rt<<1|1
#define root l,r,rt
const int maxn =3e4+5;
const int mod=1e9+7;
ll powmod(ll x,ll y){ll t; for(t=1;y;y>>=1,x=x*x%mod) if(y&1) t=t*x%mod; return t;}
ll gcd(ll x,ll y){return y?gcd(y,x%y):x;}
/*
题目大意;给定一个序列和几个查询，每次查询一段区间其有多少种不同的组合，
因为有重复的数所以要去个重复度。组合数学+莫队去处理。
我刚开始没用逆元疯狂WA，
后来想明白了，在取模的环境下很多性质都变了，
比如可能答案不是整除关系，这样就有影响了。
这样在原来的代码中加个组合递推逆元模板就OK了。*/
int d[maxn],blocks;
struct qy
{int l,r,id;bool operator<(const qy& y) const{if(l/blocks==y.l/blocks) return r<y.r;return l/blocks<y.l/blocks;}
}q[maxn];
ll ans[maxn],cnt[maxn];
int n,m;
///组合数逆元和阶乘
ll inv[maxn],fac[maxn];
void init()
{fac[0]=1;for(int i=1;i<maxn;i++) fac[i]=fac[i-1]*i%mod;inv[0]=inv[1]=1; for(int i=2;i<maxn;i++) inv[i]=(mod-mod/i)*inv[mod%i]%mod;
}int main()
{init();int t;scanf("%d",&t);while(t--){scanf("%d%d",&n,&m);for(int i=1;i<=n;i++){scanf("%d",&d[i]);}blocks=sqrt(n+0.5);for(int i=0;i<m;i++){scanf("%d%d",&q[i].l,&q[i].r);q[i].id=i;}sort(q,q+m);///对询问排序int l=1,r=0;ll ret=1;memset(cnt,0,sizeof(cnt));for(int i=0;i<m;i++){while(r<q[i].r){r++;cnt[d[r]]++;ret=ret*inv[cnt[d[r]]]%mod;}while(r>q[i].r){ret=ret*cnt[d[r]]%mod;cnt[d[r]]--;r--;}while(l<q[i].l){ret=ret*cnt[d[l]]%mod;cnt[d[l]]--, l++;}while(l>q[i].l){l--; cnt[d[l]]++;ret=ret*inv[cnt[d[l]]]%mod;}ans[q[i].id]=ret*fac[r-l+1]%mod;}for(int i=0;i<m;i++) printf("%lld\n",ans[i]);}return 0;
}