zoj-4053(2018ICPC青岛网络赛K题)启发式分裂

题意:每次删除区间内的一个数导致这个区间不连续,让你求出连续区间的最大逆序对,因为这题是强制在线的,所以没办法莫队,而且卡常卡的比较斤也没办法分块,所以我们可以使用一个启发式分裂,每次分裂计算两段区间的逆序对然后插入一个multiset容器里面最后查询接着查询最大值,我每次选取比较小的一块区间计算我这个小区间有多少个逆序对,最后原本区间逆序对减去p对左区间和右区间的贡献加上比较小的区间和比较大的区间里面选两个数的逆序对数.剩下的就是单独左间的贡献加上单独右区间的贡献,那么你再计算比较小的区间的贡献,最后再减一减就是贡献,

每次选取比较小的区间用主席树求逆序对时间复杂度大概是O(nlog^2n),最后用一个v数组记录一下区间以这个点开头的逆序对是多少我们就可以做到O(1)查询某段区间的逆序对,删除的话是nlog^2

#include<iostream>
#include<cstring>
#include<algorithm>
#include<queue>
#include<vector>
#include<cstdio>
#include<cmath>
#include<set>
#include<map>
#include<complex>
#include<cstdlib>
#include<ctime>
#include<stack>
#include<bitset>
using namespace std;
#define mes(a,b) memset(a,b,sizeof(a))
#define rep(i,a,b) for(int i = a; i <= b; i++)
#define dec(i,a,b) for(int i = a; i >= b; i--)
#define pb push_back
#define mk make_pair
#define fi first
#define se second
#define ls rt<<1
#define rs rt<<1|1
#define lson ls,L,mid
#define rson rs,mid+1,R
#define lowbit(x) x&(-x)
typedef double db;
typedef long long int ll;
typedef pair<int,int> pii;
typedef complex<double> cd;
typedef unsigned long long ull;
const ll inf = 0x3f3f3f3f;
const int mx = 1e5+5;
const int mod = 1e9+7;
const int x_move[] = {1,-1,0,0,1,1,-1,-1};
const int y_move[] = {0,0,1,-1,1,-1,1,-1};
int n,m;
struct node{
	int l,r;
	int sum;
}T[mx<<6];
int cnt;
int id[mx];
int a[mx];
ll v[mx<<2];
ll ans;
void built(int &rt,int l,int r){
	rt = cnt++;
	T[rt].sum = 0;
	if(l==r)
		return ;
	int mid = (l+r)/2;
	built(T[rt].l,l,mid);
	built(T[rt].r,mid+1,r);
}
void update(int prt,int &rt,int l,int r,int p){
	rt = cnt++;
	T[rt] = T[prt];
	T[rt].sum++;
	if(l==r)
		return;
	int mid = (l+r)/2;
	if(p>mid) update(T[prt].r,T[rt].r,mid+1,r,p);
	else update(T[prt].l,T[rt].l,l,mid,p);
}
int query(int prt,int rt,int l,int r,int p){
	if(p>=r)
		return 0;
	if(p<l)
		return T[rt].sum-T[prt].sum;
	int mid = (l+r)/2;
	return query(T[prt].r,T[rt].r,mid+1,r,p)+query(T[prt].l,T[rt].l,l,mid,p);
}
int query1(int prt,int rt,int l,int r,int p){
	if(p<=l)
		return 0;
	if(p>r)
		return T[rt].sum-T[prt].sum;
	int mid = (l+r)/2;
	return query1(T[prt].r,T[rt].r,mid+1,r,p)+query1(T[prt].l,T[rt].l,l,mid,p);
}
set<int>st;
multiset<ll>mst;
void del(int l,int r,int L,int R,int p){
	ll sum = 0;
	if(l!=r){
		if(R<l) sum += query1(id[l-1],id[r-1],1,n,a[p]);
		else sum += query(id[l-1],id[r-1],1,n,a[p]);
	}
	if(L!=R){
		if(R<l)
			sum += query(id[L-1],id[R-1],1,n,a[p]);
		else
			sum += query1(id[L-1],id[R-1],1,n,a[p]);
	}
	mst.erase(mst.find(ans));
	for(int i = l; i < r; i++)
		if(R<l)
			sum += query(id[L-1],id[R-1],1,n,a[i]);
		else
			sum += query1(id[L-1],id[R-1],1,n,a[i]);
	ans -= sum;
	sum = 0;
	for(int i = l; i < r; i++)
		sum += query(id[l-1],id[i],1,n,a[i]);	
	ans -= sum;
	if(l!=r)
		v[l] = sum,mst.insert(sum);
	if(L!=R)
		v[L] = ans,mst.insert(ans);
}
int main(){
	//freopen("1.in","r",stdin);
	//freopen("1.out","w",stdout);
	int t,q,ca = 1;
	scanf("%d",&t);
	while(t--){
		scanf("%d",&n);
		for(int i = 1; i <= n; i++)
			scanf("%d",&a[i]);
		cnt = 1;
		built(id[0],1,n);
		ans = 0;
		for(int i = 1; i <= n; i++){
			update(id[i-1],id[i],1,n,a[i]);
			ans += query(id[0],id[i],1,n,a[i]);
		}
		v[1] = ans;
		st.clear();
		st.insert(0);
		st.insert(n+1);
		mst.clear();
		mst.insert(ans);
		for(int i = 1; i <= n; i++){
			ll p;
			scanf("%lld",&p);
			printf("%lld%c",*(--mst.end()),i==n?'\n':' ');
			p = p^(*(--mst.end()));
			auto it2 = st.lower_bound(p);
			if(p<1||p>n||*it2==p)
				continue;
			auto it1 = it2;
			it1--;
			ans = v[*it1+1];
			if(p-*it1 > *it2-p)
				del(p+1,*it2,*it1+1,p,p);
			else
				del(*it1+1,p,p+1,*it2,p);
			st.insert(p);
		}
	}
	return 0;
}