传送门
对于区间查询且不带修改,加上多重集合的公式,我们很容易想到莫队,边界的转移也很容易写出来。(多重集全排列 n ! / ( a i ! + . . . + a n ! ) n!/(a_i!+...+a_n!) n!/(ai!+...+an!))
#include<bits/stdc++.h>
using namespace std;
typedef long long ll;
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;}
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;
}