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=;

 vector< pair<int,int> > tr[*];

 int laz[*],s[*],a[],b[],c[];

 int A,B,n,m,t;

 int C = ~(<<), M = (<<)-;

 int rnd(int last,int &a,int &b)

 {

     a = ( + (last >> )) * (a & M) + (a >> );

     b = ( + (last >> )) * (b & M) + (b >> );

     return (C & ((a << ) + b)) % ;

 }

 void build(int i,int l,int r)

 {

     laz[i]=-;

     tr[i].clear();

     if (l==r)

     {

         s[i]=(a[l]>=b[l]);

     }

     else {

         int m=(l+r)>>;

         build(i*,l,m);

         build(i*+,m+,r);

         s[i]=s[i*]+s[i*+];

         int x=l,y=m+,h=l;

         while (x<=m&&y<=r)

         {

             if (b[x]<b[y])

             {

                 tr[i].push_back(mp(x-l,y-m-));

                 c[h++]=b[x++];

             }

             else {

                 tr[i].push_back(mp(x-l-,y-m-));

                 c[h++]=b[y++];

             }

         }

         while (x<=m)

         {

             tr[i].push_back(mp(x-l,r-m-));

             c[h++]=b[x++];

         }

         while (y<=r)

         {

             tr[i].push_back(mp(m-l,y-m-));

             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*]=(x==-)?-:tr[i][x].first;

     laz[i*+]=(x==-)?-:tr[i][x].second;

     s[i*]=laz[i*]+;

     s[i*+]=laz[i*+]+;

     laz[i]=-;

 }

 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+;

     }

     else {

         int m=(l+r)>>;

         if (laz[i]>-) push(i);

         if (x<=m) add(i*,l,m,x,y,w==-?w:tr[i][w].first);

         if (y>m) add(i*+,m+,r,x,y,w==-?w:tr[i][w].second);

         s[i]=s[i*]+s[i*+];

     }

 }

 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)>>,ans=;

         if (laz[i]>-) push(i);

         if (x<=m) ans+=ask(i*,l,m,x,y);

         if (y>m) ans+=ask(i*+,m+,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=; i<=n; i++) scanf("%d",&a[i]);

         for (int i=; i<=n; i++) scanf("%d",&b[i]);

         b[n+]=;

         build(,,n);

         int ans=;

         int last=;

         for (int i=; i<=m; i++)

         {

             int l=rnd(last,A,B)%n+,r=rnd(last,A,B)%n+,x=rnd(last,A,B)+;

             if (l>r) swap(l,r);

             if ((l+r+x)%==)

             {

                 last=ask(,,n,l,r);

                 ans=(ans+1ll*i*last%mo)%mo;

             }

             else {

                 int w=upper_bound(b+,b++n,x)-b-;

                 add(,,n,l,r,w);

             }

         }

         printf("%d\n",ans);

     }

 }