排序最小交换次数

1.只能相邻交换，冒泡排序O(n^2)-逆序对O(nlogn)
2.允许自由交换，无重复数值，选择排序O(n^2)-预排序记录目标位置O(nlogn)
3. 允许自由交换，有重复数值

//author:uni_xxxyd
//主要功能为求序列到达有序/到达目标序列的最少自由交换次数（相邻交换参考逆序对算法）
//数值可以有重复 
//复杂度nlogn 可过1e6 

//输入格式：
//n
//n个数 原序列
//*n个数 目标序列 

#include<cstdio>
#include<set>
#include<iostream>
#include<cstring>
#include<algorithm>
#define debug cout<<"****************"<<endl
using namespace std;
typedef long long ll;
inline void ju() {
	freopen("***.in","r",stdin);
	freopen("***.out","w",stdout);
}
inline ll rd() {
	ll x=0,f=1; char c=getchar();
	while((c<'0'||c>'9')&&(c!='-')) c=getchar();
	if(c=='-') f=-1,c=getchar();
	while(c>='0'&&c<='9') {
		x=(x<<3)+(x<<1)+c-'0';
		c=getchar();
	}
	return x*f;
}
const int N=1e6+10;
struct spair{
	int x,y;
} a[N];
bool cmp(spair l,spair r) {
	return l.y<r.y;
}
int b[N],l[N],r[N];
set<int> S[N];
int main() {
	int n=rd(),ans=0;
	for(int i=1;i<=n;++i) a[i].x=rd(); //原序列 
	for(int i=1;i<=n;++i) a[i].y=rd(); //目标序列，无目标则改为升/降序序列（可用其它nlogn算法完成），须保证和原序列仅有位置上的差异 
	//for(int i=1;i<=n;++i) a[i].y=i;
	sort(a+1,a+n+1,cmp);
	
	//离散化，b数组为去重后的a[].y 
	for(int i=1;i<=n;++i) b[i]=a[i].y;
	int m=unique(b+1,b+n+1)-b-1;
	
	//for(int i=1;i<=n;++i) printf("%d %d\n",a[i].x,a[i].y);
	
	//预处理各数值应位于的区间 
	for(int k=1;k<=m;++k) {
		l[b[k]]=n;
		r[b[k]]=1;
	}
	for(int i=1;i<=n;++i) {
		l[a[i].y]=min(l[a[i].y],i);
		r[a[i].y]=max(r[a[i].y],i);
		if(a[i].x!=a[i].y) S[a[i].x].insert(i);
	}
	
	//开始排序 
	for(int k=1;k<=m;++k) {
		//printf("%d %d\n",l[b[k]],r[b[k]]);
		
		//第一遍匹配，尝试寻找交换后双方都能到达正确位置的一组交换 
		for(int i=l[b[k]];i<=r[b[k]];++i) {
			if(a[i].x!=a[i].y) {
				auto it=S[a[i].y].lower_bound(l[a[i].x]);
				if(it!=S[a[i].y].end()&&*it<=r[a[i].x]) {
					int p=*it;
					S[a[i].y].erase(it);
					S[a[i].x].erase(i);
					S[a[i].x].insert(p);
					char tmp=a[i].x;
					a[i].x=a[p].x;
					a[p].x=tmp;
					++ans;
				}
			}
		}
		
		//第二遍匹配，交换后只有当前位置上的数正确，被换走的数等待之后的交换 
		for(int i=l[b[k]];i<=r[b[k]];++i) {
			if(a[i].x!=a[i].y) {
				auto it=S[a[i].y].begin();
				while(it!=S[a[i].y].end()&&a[*it].x==a[*it].y) S[a[i].y].erase(it++);
				int p=*it;
				S[a[i].y].erase(it);
				S[a[i].x].erase(i);
				S[a[i].x].insert(p);
				char tmp=a[i].x;
				a[i].x=a[p].x;
				a[p].x=tmp;
				++ans;
			}
		}
	}
	printf("%d\n",ans);
	return 0;
}