分块，莫队算法浅析

分块 CDoj 1324

#include<bits/stdc++.h>
using namespace std;
typedef long long ll;
const int maxn=1e5+7;
int block,belong[maxn],num,l[maxn],r[maxn],n,q;
ll a[maxn],Max[maxn];
void build(){
    block=sqrt(n);
    num=n/block;if(n%block)num++;
    for(int i=1;i<=num;i++){
        l[i]=(i-1)*block+1;
        r[i]=i*block;
    }
    r[num]=n;
    for(int i=1;i<=n;i++){
        belong[i]=(i-1)/block+1;
    }
    for(int i=1;i<=num;i++){//首先O(n)一遍寻找每一块的最大值
        for(int j=l[i];j<=r[i];j++){
            Max[i]=max(Max[i],a[j]);
        }
    }
}
void update(int x,int y){
    a[x]+=y;
    Max[belong[x]]=max(Max[belong[x]],a[x]);
}
ll query(int x,int y){
    ll ans=0;
    if(belong[x]==belong[y]){
        for(int i=x;i<=y;i++){
            ans=max(a[i],ans);
        }
        return ans;
    }
    for(int i=x;i<=r[belong[x]];i++)
        ans=max(ans,a[i]);
    for(int i=belong[x]+1;i<belong[y];i++)
        ans=max(ans,Max[i]);
    for(int i=l[belong[y]];i<=y;i++)
        ans=max(ans,a[i]);
    return ans;
}
int main(){
    scanf("%d%d",&n,&q);
    build();int op,l,r;
    for(int i=1;i<=q;i++){
        scanf("%d%d%d",&op,&l,&r);
        if(op==1)update(l,r);
        else printf("%lld\n",query(l,r));
    }
}

莫队

codeforces E - XOR and Favorite Number

#include<bits/stdc++.h>
using namespace std;
const int maxn=2e6+5;
typedef long long ll;
struct node{
    int l,r,id;
}q[maxn];
int a[maxn],num[maxn],belong[maxn];
ll sum,ans[maxn];
int n,m,k;
bool cmp(node a,node b){
    if(belong[a.l]==belong[b.l])
        return a.r<b.r;
    return belong[a.l]<belong[b.l];
}
void add(int x){
    sum+=num[a[x]^k];
    num[a[x]]++;
}
void del(int x){
    num[a[x]]--;
    sum-=num[a[x]^k];
}
int main(){

    cin>>n>>m>>k;

    int block=sqrt(n);
    for(int i=1;i<=n;i++){
        cin>>a[i];
        a[i]=a[i]^a[i-1];
        belong[i]=i/block;
    }
    for(int i=1;i<=m;i++){
        cin>>q[i].l>>q[i].r;
        q[i].id=i;
    }
    num[0]=1;
    sort(q+1,q+m+1,cmp);
    int l=0,r=0;
    for(int i=1;i<=m;i++){
        while(r<q[i].r){
            r++;
            add(r);
        }
        while(l+1>q[i].l){
            l--;
            add(l);
        }
        while(l+1<q[i].l){
            del(l);
            l++;
        }
        while(r>q[i].r){
            del(r);
            r--;
        }
        ans[q[i].id]=sum;
    }
    for(int i=1;i<=m;i++){
        cout<<ans[i]<<endl;
    }
    return 0;
}