CR751_C 分治DP

本题的数列A次序不变，而B次序任意，显然可以令B有序排列，当B升序时，显然B内的inversion对数最少，于是我们天然的能够想到令B升序很可能是最优的选择。

而实际上，它确实是，比如对于序列12，13，1，5 有当13与1交换时inversion数可以降低1，而对任意的序列a1,a2,...,ai,...,aj,...,an有ai>aj则将ai与aj互换位置之后，inversion的变化量主要发生在i<=t<=j这一区间中，本不与ai构成inversion对的，在新的序列中仍然不与aj构成inversion对，与aj不构成inversion对的，在新序列中也仍不与ai构成inversion对，于是显然，这样的交换是使序列的inversion对不增的，于是显然当序列B升序的时候，出现最优情况。

于是我们容易证明这样的序列具有决策单调性，如图：

 于是，我们可以采取分治DP的策略计算决策变量。

#include<bits/stdc++.h>
using namespace std;
typedef long long ll;
struct add_bint{
	#define MAXN 2000010
	ll D1[MAXN],D2[MAXN];
	ll Num;
	void _update(int i,ll k)
	{
		int x=i;
		while(i<=Num)
		{
			D1[i]+=k;
			D2[i]+=k*(x-1);
			i+=i&-i;
		}
	}
	ll _getsum(int i)
	{
		ll res=0;
		int x=i;
		while(i)
		{
			res+=x*D1[i]-D2[i];
			i-=i&-i;
		}
		return res;
	}
	void build(int l,int r){
		Num=r+1;
		for(int i=0;i<=r+1;++i){
		    D1[i]=D2[i]=0;
		}
	}
	inline void add(int l,int r,ll k){
		_update(l,k);
		_update(r+1,-k);
	}
	inline ll sum(int l,int r){
		return _getsum(r)-_getsum(l-1);
	}
}bul;
#define maxn 2000010
ll n,m;
ll A[maxn],B[maxn];
ll OPT[maxn];
pair<ll,ll> ans[maxn];
typedef pair<ll,ll> pii;
bool sortbysecond(const pii&x,const pii&y){
	if(x.second==y.second){
		return x.first<y.first;
	}
	return x.second<y.second;
}
bool sortbyfirst(const pii&x,const pii&y){
	if(x.first==y.first){
		return x.second<y.second;
	}
	return x.first<y.first;
}
int id;
void descret(pii*al,int num){
	sort(al,al+num,sortbysecond);
	id=0;
	ll val=-1;
	for(int i=0;i<num;++i){
		if(al[i].second==val){
			al[i].second=id;
		}
		else{
			val=al[i].second;
			al[i].second=++id;
		}
	}
	sort(al,al+num,sortbyfirst);
}
void DC_DP(int l,int r,int optl,int optr){
	if(l>r){
		return;
	}
	int mm=l+r>>1;
	int opt=-1;
	int deno=1000000000;
	int tem=0;
	for(int i=optl;i<optr;++i){
		if(A[i]<B[mm])++tem;
	}
	for(int i=optl;i<=optr;++i){
		if(tem<deno){
			deno=tem;
			opt=i;
		}
		if(i==optr) break;
		if(A[i]<B[mm])--tem;
		if(A[i]>B[mm])++tem;
	}
	OPT[mm]=opt;
	DC_DP(l,mm-1,optl,opt);
	DC_DP(mm+1,r,opt,optr);
}
void work(){
	scanf("%lld%lld",&n,&m);
	for(int i=1;i<=n;++i){
		scanf("%lld",&A[i]);
	}
	for(int i=1;i<=m;++i){
		scanf("%lld",&B[i]);
	}
	sort(B+1,B+m+1);
	DC_DP(1,m,1,n+1);
	int denote=0;
	for(int i=1;i<=n;++i){
		ans[denote++]={2*i+1,A[i]};
	}
	for(int i=1;i<=m;++i){
		ans[denote++]={2*OPT[i],B[i]};
	}
	descret(ans,denote);
	ll result=0;
	bul.build(1,id+2);
	for(int i=0;i<denote;++i){
		result+=bul.sum(ans[i].second+1,id+1);
		bul.add(ans[i].second,ans[i].second,1);
	}
	printf("%lld\n",result);
}
 
int main(){
	ll T;
	scanf("%lld",&T);
	while(T--){
		work();
	}
}