hdu5145——莫队+多重集全排列

传送门

对于区间查询且不带修改，加上多重集合的公式，我们很容易想到莫队，边界的转移也很容易写出来。(多重集全排列$n!/(a_i!+...+a_n!)$)

#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;
}