【洛谷 P1966 】火柴排队【逆序对】【数学（映射）】

---

解题思路

好了回归题目，既然要求找到 $min(\Sigma [(ai-bi{)}^2])$，那么我们不妨变形一下：
$\Sigma [(ai-bi{)}^2]=\Sigma (a{i}^2-2\ast ai\ast bi+b{i}^2)$
由此可见，若想使$min(\Sigma [(ai-bi{)}^2])$最小，而和式中$a{i}^2+b{i}^2$是个定值，那么，就只能想办法让 $\Sigma 2\ast aibi$ 这一项最大

经过排序+离散化之后，假设得到了这样两个序列：
$A:(1,3,4,2)和B:(1,4,2,3)$。
我们要让A的第1大对应B的第1大，A的第2大对应B的第2大，A的第n大对应B的第n大.

我们以$a[i]$为关键字对$b[i]$排序，令$x[a[i]]=b[i]$，在a和b中构造一种映射关系。

若序列a与序列b相等，那么此时$x[a[i]]$应该等于a[i]的，也就是$x[i]=i$。那么也就是说如果我们想让序列a与序列b相等，那么我们需要让 $x$升序排列。

问题就变为：将原本乱的X序列升序排列的最少交换次数=求X序列的逆序对个数。

---

代码

#include<cstdio>
#include<iostream>
#include<queue>
#include<algorithm>
#include<cstring>
#include<cmath>
#define ll long long
using namespace std;

const int mod=1e8-3;
int n,ans,w[1000010],s[1000010];
struct c{
	int x,id;
}a[100010],b[100010];

bool cmp(c l,c r)
{
	return l.x<r.x;
}

void work(int l,int r)
{
	int mid=(l+r)/2;
	int i=l,j=mid+1,tot=l;
	while(i<=mid&&j<=r)
	{
		if(w[i]<=w[j])
			s[tot++]=w[i++];
		else 
		{
			s[tot++]=w[j++];
			ans=(ans+mid-i+1)%mod;
		}
	}
	while(i<=mid)
		s[tot++]=w[i++];
	while(j<=r)
		s[tot++]=w[j++];
	for(int i=l;i<=r;i++)
		w[i]=s[i];
}

void px(int l,int r)
{
	if(l==r)return;
	int mid=(l+r)/2;
	px(l,mid);
	px(mid+1,r);
	work(l,r);
}

int main(){
	scanf("%d",&n);
	for(int i=1;i<=n;i++)
	{
		scanf("%d",&a[i].x);
		a[i].id=i;
	}
	for(int i=1;i<=n;i++)
	{
		scanf("%d",&b[i].x);
		b[i].id=i;
	}
	sort(a+1,a+n+1,cmp);
	sort(b+1,b+n+1,cmp);
	for(int i=1;i<=n;i++)
		w[b[i].id]=a[i].id;
	px(1,n);
	printf("%d",ans);
}