莫队算法

/*当我们知道s(n,m)的值之后便可以通过它得到四个递推公式：s（n,m-1）=s(n,m)-c(n,m);
s(n,m+1)=s(n,m)+c(n,m+1);2*s(n-1.m)=s(n,m)+c(n-1,m);s(n+1,m)=2*s(n,m)-c(n,m).
由此，如果我们知道S(n,m)，我们就可以知道在O(1)的时间复杂度内求出S(n+1,m)、S(n-1,m)、S(n,m+1)
和S(n,m-1)这四个的数值。那么，我们可以把n和m看作区间询问的左右端点，对于读入的T个询问，我们可以按照
莫队算法的区间移动的方式来解决这道题，之后就是莫队算法模板题了。*/
/*代码*/
#include<bits/stdc++.h>
using namespace std;
#define ll long long
const ll mod=1e9+7;
const int maxn=1e5+10;
ll fac[maxn],inv[maxn],ans[maxn];
struct node
{
    int n,m,id;
};
vector<node>G[maxn];
bool cmp(const node &a,const node &b)
{
    return a.n<b.n;
}
ll pow1(ll a,ll b)
{
    ll r=1;
    while(b)
    {
        if(b&1) r=r*a%mod;
        a=a*a%mod;
        b/=2;
    }
    return r;
}
void init()
{
    fac[0]=1;
    for(ll i=1; i<maxn; i++) fac[i]=fac[i-1]*i%mod;
    inv[maxn-1]=pow1(fac[maxn-1],mod-2);
    for(ll i=maxn-2; i>=0; i--) inv[i]=inv[i+1]*(i+1)%mod;
}
ll C(int n,int m)
{
    return fac[n]*inv[m]%mod*inv[n-m]%mod;
}
int main()
{
    init();
    int T;
    scanf("%d",&T);
    int unit=sqrt(maxn);
    for(int i=1; i<=T; i++)
    {
        node e;
        e.id=i;
        scanf("%d%d",&e.n,&e.m);
        G[e.m/unit+1].push_back(e);
    }
    for(int i=1; i<=unit+2; i++)
    {
        int len=G[i].size();
        if(!len)continue;
        sort(G[i].begin(),G[i].end(),cmp);
        int x=G[i][0].n,y=0;
        ll cnt=1;
        int ss=G[i].size();
        for(int j=0; j<ss; j++)
        {
            node e=G[i][j];
            while(x<e.n) cnt=(2*cnt-C(x++,y)+mod)%mod;
            while(y<e.m) cnt=(cnt+C(x,++y))%mod;
            while(y>e.m) cnt=(cnt-C(x,y--)+mod)%mod;
            ans[e.id]=cnt;
        }
    }
    for(int i=1; i<=T; i++)
        printf("%lld\n",ans[i]);
    return 0;
}