hdu 5145 NPY and girls 莫队

题意：

给你1-n属于的班级

给你一个[l,r]区间

问你如果要访问这个区间内所有的女生

有多少种走不同教室的方法

思路：

和小z差不多 只不过这个维护的是阶乘

找出来公式之后莫队直接离线处理

莫队更多的是离线排序优化的思想

把所有查询排序处理 然后逐个处理 可以应用到很多方面

#include<bits/stdc++.h>
#define cl(a,b) memset(a,b,sizeof(a))
#define debug cerr<<#a<<"=="<<a<<endl;
using namespace std;
typedef long long ll;
typedef pair<int,int> pii;

const int maxn=3e4+10;
const int mod=1e9+7;

int n,m;
ll inv[maxn];
ll a[maxn];
ll res[maxn];
ll num[maxn];

struct query
{
    int l,r,id,sqr;
    bool operator < (const query &a) const
    {
        if(a.sqr==sqr) return r<a.r;
        return l<a.l;
    }
} q[maxn];

ll quick_mod(ll a,ll x)
{
    ll ans=1;
    while(x)
    {
        if(x&1)
        {
            ans=(ans*a)%mod;
        }
        a=(a*a)%mod;
        x>>=1;
    }
    return ans;
}

ll getinv(ll x)
{
    return quick_mod(x,mod-2);
}

void init()
{
    for(int i=1; i<=maxn; i++)
    {
        inv[i]=getinv(i);
    }
}

void solve()
{
    int l=1,r=1;
    ll ans=1;
    cl(num,0);
    num[a[1]]++;
    for(int i=1; i<=m; i++)
    {
        while(r<q[i].r)
        {
            r++;
            num[a[r]]++;
            ans=ans*(r-l+1)%mod*inv[num[a[r]]]%mod;
        }
        while(l>q[i].l)
        {
            l--;
            num[a[l]]++;
            ans=ans*(r-l+1)%mod*inv[num[a[l]]]%mod;
        }
        while(r>q[i].r)
        {
            ans=ans*num[a[r]]%mod*inv[r-l+1]%mod;
            num[a[r]]--;
            r--;
        }
        while(l<q[i].l)
        {
            ans=ans*num[a[l]]%mod*inv[r-l+1]%mod;
            num[a[l]]--;
            l++;
        }
        res[q[i].id]=ans;
    }
}

int main()
{
    int T;
    scanf("%d",&T);
    init();
    while(T--)
    {
        scanf("%d%d",&n,&m);
        int bk=(int)sqrt(1.0*n);
        for(int i=1; i<=n; i++)
        {
            scanf("%lld",&a[i]);
        }
        int l,r;
        for(int i=1; i<=m; i++)
        {
            scanf("%d%d",&l,&r);
            q[i]= {l,r,i,l/bk};
        }
        sort(q+1,q+m+1);
        solve();
        for(int i=1;i<=m;i++)
        {
            printf("%lld
",res[i]);
        }
    }
    return 0;
}/*

2
4 2
1 2 1 3
1 3
1 4
1 1
1
1 1

*/