預處理後主席樹維護

首先得出最後的答案為 \(\sum_{i=l}^{r}{min(right[i],r)-i+1}\)

step1

那麽首要問題就是如何求出\(right[i]\)

考慮當i--j-1是上升時使區間i--j是上升的

即sum[i-1]^sum[j-1]<=sum[i-1]^sum[j]

觀察到兩邊有差異的是sum[j-1]和sum[j] 也就意味著sum[j-1]和sum[j]的不同會對i的取值有限制

假設k為二進制下sum[j-1]與sum[j]最高的不同位
如果sum[j]此位為1對i的限制是sum[i-1]的此位不能為1
如果sum[j]此位為0對i的限制是sum[i-1]的此位不能為0

通過枚舉每一位的限制即可得\(ri[i]\)的最大合理值

step2

接下來就是利用主席數維護答案了

\(\sum_{i=l}^{r}{min(right[i],r)-i+1}\)

我們可以對於所有的\(ri[i]\)建設主席數 維護兩個值

1.所有\(ri[i]\)在i--j的範圍內總和sum

2.所有\(ri[i]\)在i--j的範圍內有幾個cnt

最後的答案及為l--r內\(ri[i]\)的值在l--r內的sum+l--r內\(ri[i]\)的值大於r的cnt\(\times\)r-l--r所有數字和+(r-l+1)

#include<bits/stdc++.h>
#define ll long long
using namespace std;
const ll maxn=10000005;
ll n,t,Q,x,y,l,r,tot,ans;
ll root[maxn],ri[maxn],lf[maxn],cnt[maxn],sum[maxn],a[maxn],f[maxn];
ll p[32][2];
ll read()
{
    ll ch=0,x=0;while(ch=getchar(),ch<‘0‘||ch>‘9‘);
    while(x=x*10+ch-48,ch=getchar(),ch>=‘0‘&&ch<=‘9‘);
    return x;
}
ll build(ll l,ll r)
{
    ll rt=++tot;
    if(l<r)
    {
        ll mid=(l+r)>>1;
        lf[rt]=build(l,mid);
        ri[rt]=build(mid+1,r);
    }
    return rt;
}
ll updata(ll pre,ll l,ll r,ll t)
{
    ll rt=++tot;lf[rt]=lf[pre];ri[rt]=ri[pre];sum[rt]=sum[pre]+t;cnt[rt]=cnt[pre]+1;
    if(l<r)
    {
        ll mid=(l+r)>>1;
        if(t<=mid)lf[rt]=updata(lf[pre],l,mid,t);
        else ri[rt]=updata(ri[pre],mid+1,r,t);
    }
    return rt;
}
ll getsum(ll x,ll y,ll l,ll r,ll L,ll R)
{
    if(L<=l&&r<=R)return sum[y]-sum[x];
    ll mid=(l+r)>>1,Tans=0;
    if(L<=mid)Tans+=getsum(lf[x],lf[y],l,mid,L,R);
    if(R>mid)Tans+=getsum(ri[x],ri[y],mid+1,r,L,R);
    return Tans;
}
ll getcnt(ll x,ll y,ll l,ll r,ll L,ll R)
{
    if(L<=l&&r<=R)return cnt[y]-cnt[x];
    ll mid=(l+r)>>1,Tans=0;
    if(L<=mid)Tans+=getcnt(lf[x],lf[y],l,mid,L,R);
    if(R>mid)Tans+=getcnt(ri[x],ri[y],mid+1,r,L,R);
    return Tans;
}
ll Sum(ll r,ll l)
{
    return r*(r-1)/2-l*(l-1)/2;
}
int main()
{
    n=read();t=read();
    for(ll i=1;i<=n;i++)a[i]=read(),a[i]^=a[i-1];
    memset(p,63,sizeof(p));Q=read();
    for(ll i=n;i>=1;i--)
    {
        f[i]=n;
        for(ll j=30;j>=0;j--)f[i]=min(f[i],p[j][(a[i-1]>>j)&1]-1);
        for(ll j=30;j>=0;j--)if(((a[i]>>j)&1)^((a[i-1]>>j)&1)){
            p[j][(a[i]>>j)&1]=min(p[j][(a[i]>>j)&1],i);break;
        }
    }
    root[0]=build(1,n);
    for(ll i=1;i<=n;i++)root[i]=updata(root[i-1],1,n,f[i]);
    for(ll i=1;i<=Q;i++)
    {
        x=read();y=read();
        x=(x+ans*t)%n+1;y=(y+ans*t)%n+1;l=min(x,y);r=max(x,y);
        //printf("%d %d %d\n",getsum(root[l-1],root[r],1,n,l,r),getcnt(root[l-1],root[r],1,n,r+1,n),Sum(r,l));
        ans=getsum(root[l-1],root[r],1,n,l,r)+r*getcnt(root[l-1],root[r],1,n,r+1,n)-Sum(r,l-1);
        printf("%lld\n",ans);
    }
    return 0;
}
 

[codechef]SnackDown 2017 Online Elimination Round Prefix XOR