牛客5 G Glad You Came

長度爲n的數組，m次操作， 每次操作把一個區間內小於v的數更新成v

m很大，但是據說mlogn能卡過去，也只能說一句隨機數據真的蠻快的

正解是rmq的變形？

用倍增的思想把一個區間分成nlogn個區間，這樣每一個區間都能分成2個區間。每次更新直接更新他所在的區間。

全部跟新完後按區間從大到小逆向更新一遍。

真是神奇的操作。

#include <bits/stdc++.h>

using namespace std;
#define go(i,a,b) for(int i=(a);i<=(b);i++)
#define us unsigned int
#define ll long long
#define N 100005
us x,y,z,w;
us rng(){
    x^=x<<11;
    x^=x>>4;
    x^=x<<5;
    x^=x>>14;
    w=x^y^z;
    x=y;y=z;z=w;
    return w;
}
int dp[N][23],mm[N],n,m;
int main()
{
    int T;
    mm[0]=-1;for(int i=1;i<N;i++)mm[i]=((i&(i-1))==0)?mm[i-1]+1:mm[i-1];
    cin>>T;
    while(T--){
        cin>>n>>m>>x>>y>>z;
        go(i,1,n)go(j,0,mm[n])dp[i][j]=0;
        while(m--){
            int l=rng()%n+1,r=rng()%n+1,v=rng()%(1<<30);

            if(l>r)swap(l,r);
            int k=mm[r-l+1];
            dp[l][k]=max(dp[l][k],v);
            dp[r-(1<<k)+1][k]=max(dp[r-(1<<k)+1][k],v);
        }
        for(int j=mm[n];j;j--)
            for(int i=1;i+(1<<j)-1<=n;i++){
                dp[i][j-1]=max(dp[i][j-1],dp[i][j]);
                dp[i+(1<<(j-1))][j-1]=max(dp[i+(1<<(j-1))][j-1],dp[i][j]);
            }
        ll ans=0;
        go(i,1,n)ans^=1ll*i*dp[i][0];
        cout<<ans<<endl;
    }
    return 0;
}