NPY and girls

NPY and girls

题目链接：http://acm.hdu.edu.cn/showproblem.php?pid=5145

莫队算法

注意到没有修改区间的操作，使用莫队算法：将整个区间分成若干个块，将询问区间按块优先升序排序，同块内按区间右界升序排序，添加一个元素，满足条件的值sum就变为sum=(sum*times[a[r]]/(r-l+1));减少一个值同理。从第一个区间一直推到最后一个区间即可。（之前TLE一直以为是区间大小的问题，现在发现是这句话while(x<0)x+=M，不过因为这题M比较大，改改区间居然被我水过去了...）

代码如下：

#include<cstdio>
#include<cmath>
#include<algorithm>
#include<cstring>
#define N 30005
#define M 1000000007
#define LL long long
#define B 250/**((int)sqrt((double)n)) sqrt(n) will be TLE**/
#define mst(x) memset(x,0,sizeof(x))
using namespace std;
LL times[N];
LL T,n,m,x,y;
LL a[N];
struct nod{
    LL l,r,no,value;
}q[N];
LL exGCD(LL a,LL b){
    if(b==0){
        x=1;
        y=0;
        return a;
    }
    LL r=exGCD(b,a%b);
    LL tmp=x;
    x=y;
    y=tmp-(a/b)*y;
    return r;
}
bool cmp_by_range(nod x,nod y){
    if(x.l/B==y.l/B){
        return x.r<y.r;
    }else return (x.l/B<y.l/B);
}
bool cmp_by_no(nod x,nod y){
    return x.no<y.no;
}
void init(){
    mst(times);
    scanf("%I64d%I64d",&n,&m);
    for(LL i=0;i<n;++i)
        scanf("%I64d",&a[i]);
    for(LL i=0;i<m;++i){
        LL l,r;
        scanf("%I64d%I64d",&l,&r);
        l--,r--;
        q[i].l=l;
        q[i].r=r;
        q[i].no=i;
        q[i].value=0;
    }
    sort(q,q+m,cmp_by_range);
}
void solve(){
    LL l=q[0].l,r=q[0].r,sum=1;
    times[a[l]]++;
    for(LL i=l+1;i<=r;++i){
        times[a[i]]++;
        exGCD(times[a[i]],M);
        while(x<0)x+=M;
        sum=((sum*(i-l+1))%M*x)%M;
    }
    q[0].value=sum;
    for(LL i=1;i<m;++i){
        //l=q[i-1].l,r=q[i-1].r;
        while(r<q[i].r){
            r++;
            times[a[r]]++;
            exGCD(times[a[r]],M);
            while(x<0)x+=M;
            sum=((sum*(r-l+1))%M*x)%M;
        }
        while(r>q[i].r){
            exGCD((r-l+1),M);
            while(x<0)x+=M;
            sum=((sum*times[a[r]])%M*x)%M;
            times[a[r]]--;
            r--;
        }
        while(l<q[i].l){
            exGCD((r-l+1),M);
            while(x<0)x+=M;
            sum=((sum*times[a[l]])%M*x)%M;
            times[a[l]]--;
            l++;
        }
        while(l>q[i].l){
            l--;
            times[a[l]]++;
            exGCD(times[a[l]],M);
            while(x<0)x+=M;
            sum=((sum*(r-l+1)%M)*x)%M;
        }
        q[i].value=sum;
    }
}
int main(void){
    scanf("%I64d",&T);
    while(T--){
        init();
        solve();
        sort(q,q+m,cmp_by_no);
        for(LL i=0;i<m;++i)
            printf("%I64d
",q[i].value);
    }
}