[Luogu P2507] [BZOJ 1237] [SCOI2008]配对

洛谷传送门

BZOJ传送门

题目描述

你有 $n$ 个整数$A_i$和 $n$ 个整数$B_i$。你需要把它们配对，即每个$A_i$恰好对应一个$B_{p[i]}$。要求所有配对的整数差的绝对值之和尽量小，但不允许两个相同的数配对。例如 A = { 5 , 6 , 8 } A=\{5,6,8\} A={5,6,8}， B = { 5 , 7 , 8 } B=\{5,7,8\} B={5,7,8}，则最优配对方案是$5\acute{o}8,6\acute{o}5,8\acute{o}7$，配对整数的差的绝对值分别为$2,2,1$，和为$5$。注意，$5\acute{o}5，6\acute{o}7，8\acute{o}8$是不允许的，因为相同的数不许配对。

输入输出格式

输入格式：

第一行为一个正整数$n$，接下来是 $n$ 行，每行两个整数$A_i$和$B_i$，保证所有

$A_i$各不相同，$B_i$也各不相同。

输出格式：

输出一个整数，即配对整数的差的绝对值之和的最小值。如果无法配对，输出$-1$。

输入输出样例

输入样例#1：

3
3 65
45 10
60 25

输出样例#1：

32

输入样例#2：

3
5 5
6 7
8 8

输出样例#2：

5

说明

30%的数据满足：$n\le 10^4$

100%的数据满足：$1\le n\le 10^5$，$A_i$和$B_i$均为$1$到$10^6$之间的整数。

解题分析

如果排序后两两对应不同， 显然对应匹配最好。

当我们发现排序后$a[i]=b[i]$的时候， 显然我们只会考虑和其相邻的两个数交换匹配的情况（因为交换的越多越不优）， 所以我们一遍$dp$考虑相邻三个数就好。

代码如下：

#include <cstdio>
#include <cstring>
#include <cstdlib>
#include <cctype>
#include <algorithm>
#include <cmath>
#define R register
#define IN inline
#define W while
#define gc getchar()
#define INF 100000000000ll
#define MX 100050
#define ll long long
template <class T>
IN void in(T &x)
{
    x = 0; R char c = gc;
    for (; !isdigit(c); c = gc);
    for (;  isdigit(c); c = gc)
    x = (x << 1) + (x << 3) + c - 48;
}
using namespace std;
IN ll judge(R int a, R int b) {return a == b ? INF : abs(a - b);}
int n;
ll dp[MX], a[MX], b[MX];
int main(void)
{
    in(n);
    for (R int i = 1; i <= n; ++i) in(a[i]), in(b[i]);
    sort(a + 1, a + 1 + n), sort(b + 1, b + 1 + n);
    if(n == 1 && a[1] == b[1]) return puts("-1"), 0;
    dp[1] = judge(a[1], b[1]);
    dp[2] = min(dp[1] + judge(a[2], b[2]), judge(a[1], b[2]) + judge(a[2], b[1]));
    for (R int i = 3; i <= n; ++i)
    {
        dp[i] = dp[i - 1] + judge(a[i], b[i]);
        dp[i] = min(dp[i], dp[i - 2] + judge(a[i - 1], b[i]) + judge(a[i], b[i - 1]));
        dp[i] = min(dp[i], dp[i - 3] + judge(a[i - 1], b[i - 2]) + judge(a[i], b[i - 1]) + judge(a[i - 2], b[i]));
        dp[i] = min(dp[i], dp[i - 3] + judge(a[i], b[i - 2]) + judge(a[i - 2], b[i - 1]) + judge(a[i - 1], b[i]));
        dp[i] = min(dp[i], dp[i - 3] + judge(a[i], b[i - 2]) + judge(a[i - 1], b[i - 1]) + judge(a[i - 2], b[i]));
    }
    printf("%lld", dp[n]);
}