hdu5737

首先思考一个朴素的做法

将b[]维护成一个线段树套有序表，每次修改a[]用线段树+lazy tag

并在线段树的子区间上在有序表中二分更新这段区间中a[i]>=b[i]的值，复杂度O(nlog^2)

有没有更优的做法？

考虑在一次修改操作中，查询的数是相同的，并且b[]这个有序表始终不变

因此在生成b[]的线段树套有序表过程中（其实就是归并排序的记录）

我们维护每个数在左右子区间中名次。

这样修改的时候我们只要对b[]排好序的整体做一次二分，找到b[m]<=x<b[m+1]

把修改成x当作修改成b[m]，这样我们就能O(1)地更新每个子区间的计数，问题得解

#include<bits/stdc++.h>
#define mp make_pair
using namespace std;
const int mo=1000000007;
vector< pair<int,int> > tr[100010*4];
int laz[100010*4],s[100010*4],a[100010],b[100010],c[100010];
int A,B,n,m,t;
int C = ~(1<<31), M = (1<<16)-1;
int rnd(int last,int &a,int &b)
{
    a = (36969 + (last >> 3)) * (a & M) + (a >> 16);
    b = (18000 + (last >> 3)) * (b & M) + (b >> 16);
    return (C & ((a << 16) + b)) % 1000000000;
}

void build(int i,int l,int r)
{
    laz[i]=-2;
    tr[i].clear();
    if (l==r)
    {
        s[i]=(a[l]>=b[l]);
    }
    else {
        int m=(l+r)>>1;
        build(i*2,l,m);
        build(i*2+1,m+1,r);
        s[i]=s[i*2]+s[i*2+1];
        int x=l,y=m+1,h=l;
        while (x<=m&&y<=r)
        {
            if (b[x]<b[y])
            {
                tr[i].push_back(mp(x-l,y-m-2));
                c[h++]=b[x++];
            }
            else {
                tr[i].push_back(mp(x-l-1,y-m-1));
                c[h++]=b[y++];
            }
        }
        while (x<=m)
        {
            tr[i].push_back(mp(x-l,r-m-1));
            c[h++]=b[x++];
        }
        while (y<=r)
        {
            tr[i].push_back(mp(m-l,y-m-1));
            c[h++]=b[y++];
        }
        for (int k=l; k<=r; k++) b[k]=c[k];
    }
}

void push(int i)
{
    int x=laz[i];
    laz[i*2]=(x==-1)?-1:tr[i][x].first;
    laz[i*2+1]=(x==-1)?-1:tr[i][x].second;
    s[i*2]=laz[i*2]+1;
    s[i*2+1]=laz[i*2+1]+1;
    laz[i]=-2;
}

void add(int i,int l,int r,int x,int y,int w)
{
    if (x<=l&&y>=r)
    {
        laz[i]=w;
        s[i]=w+1;
    }
    else {
        int m=(l+r)>>1;
        if (laz[i]>-2) push(i);
        if (x<=m) add(i*2,l,m,x,y,w==-1?w:tr[i][w].first);
        if (y>m) add(i*2+1,m+1,r,x,y,w==-1?w:tr[i][w].second);
        s[i]=s[i*2]+s[i*2+1];
    }
}

int ask(int i,int l,int r,int x,int y)
{
    if (x<=l&&y>=r) return s[i];
    else {
        int m=(l+r)>>1,ans=0;
        if (laz[i]>-2) push(i);
        if (x<=m) ans+=ask(i*2,l,m,x,y);
        if (y>m) ans+=ask(i*2+1,m+1,r,x,y);
        return ans;
    }
}

int main()
{
    int cas;
    scanf("%d",&cas);
    while (cas--)
    {
        scanf("%d%d%d%d",&n,&m,&A,&B);
        for (int i=1; i<=n; i++) scanf("%d",&a[i]);
        for (int i=1; i<=n; i++) scanf("%d",&b[i]);
        b[n+1]=2147483647;
        build(1,1,n);
        int ans=0;
        int last=0;
        for (int i=1; i<=m; i++)
        {
            int l=rnd(last,A,B)%n+1,r=rnd(last,A,B)%n+1,x=rnd(last,A,B)+1;
            if (l>r) swap(l,r);
            if ((l+r+x)%2==0)
            {
                last=ask(1,1,n,l,r);
                ans=(ans+1ll*i*last%mo)%mo;
            }
            else {
                int w=upper_bound(b+1,b+1+n,x)-b-2;
                add(1,1,n,l,r,w);
            }
        }
        printf("%d
",ans);
    }
}