HYSBZ - 2038经典莫队算法题

无修改的莫队

emmm莫队的几条性质，必须是离线的，复杂度是n*sqrt(n)

就是通过预处理查询区间，然后从(l,r)转移到(ll,rr)，这样的复杂度是曼哈顿距离，即abs(l-ll)+abs(r-rr)

通过分块可以保证复杂度在n*sqrt(n)内

详细请参考：https://www.cnblogs.com/Paul-Guderian/p/6933799.html

#include<bits/stdc++.h>
#define fi first
#define se second
#define ll long long
#define mp make_pair
#define pb push_back
#define mod 1000000007
#define pii pair<int,int>
#define ls l,m,rt<<1
#define rs m+1,r,rt<<1|1
 
using namespace std;
 
const int g=10.0,eps=1e-;
const int N=+,maxn=+,inf=0x3f3f3f3f;
 
struct mo{
    int l,r,id;
    ll a,b;
}q[N];
int belong[N],co[N];
ll sum[N],ans;
bool cmp(mo a,mo b)
{
    return belong[a.l]==belong[b.l]?a.r<b.r:a.l<b.l;
}
bool cmp1(mo a,mo b)
{
    return a.id<b.id;
}
void re(int i,int add)
{
    ans-=sum[co[i]]*sum[co[i]];
    sum[co[i]]+=add;
    ans+=sum[co[i]]*sum[co[i]];
}
int main()
{
    /*ios::sync_with_stdio(false);
    cin.tie(0);*/
    int n,m;
    scanf("%d%d",&n,&m);
    int block=sqrt(n);
    for(int i=;i<=n;i++)
        scanf("%d",&co[i]),belong[i]=i/block+;
    for(int i=;i<=m;i++)
        scanf("%d%d",&q[i].l,&q[i].r),q[i].id=i;
    sort(q+,q++m,cmp);
    int l=,r=;
    for(int i=;i<=m;i++)
    {
        while(l<q[i].l)re(l,-),l++;
        while(l>q[i].l)re(l-,),l--;
        while(r<q[i].r)re(r+,),r++;
        while(r>q[i].r)re(r,-),r--;
        q[i].a=ans-(q[i].r-q[i].l+);
        q[i].b=(ll)(q[i].r-q[i].l+)*(q[i].r-q[i].l);
        ll x=__gcd(q[i].a,q[i].b);
        q[i].a/=x;q[i].b/=x;
    }
    sort(q+,q++m,cmp1);
    for(int i=;i<=m;i++)
        printf("%lld/%lld\n",q[i].a,q[i].b);
    return ;
}
/*******************
 
*******************/