HDU 5145 NPY and girls (莫队算法+组合逆元)

题目链接: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;
}