[NOIP2013]火柴排队

题目描述 涵涵有两盒火柴，每盒装有 n 根火柴，每根火柴都有一个高度。现在将每盒中的火柴各自排成一列，同一列火柴的高度互不相同，两列火柴之间的距离定义为：∑(ai-bi)^2, i=1~n，其中ai表示第一列火柴中第i个火柴的高度，bi表示第二列火柴中第i个火柴的高度。

每列火柴中相邻两根火柴的位置都可以交换，请你通过交换使得两列火柴之间的距离最小。请问得到这个最小的距离，最少需要交换多少次？如果这个数字太大，请输出这个最小交换次数对 99,999,997 取模的结果。

输入格式 共三行，第一行包含一个整数 n，表示每盒中火柴的数目。

第二行有 n 个整数，每两个整数之间用一个空格隔开，表示第一列火柴的高度。

第三行有 n 个整数，每两个整数之间用一个空格隔开，表示第二列火柴的高度。

输出格式 输出共一行，包含一个整数，表示最少交换次数对 99,999,997 取模的结果。

样例输入 4 2 3 1 4 3 2 1 4 样例输出 1 样例输入 4 1 3 4 2 1 7 2 4 样例输出 2 注释 对于 10%的数据， 1 ≤ n ≤ 10；

对于 30%的数据，1 ≤ n ≤ 100；

对于 60%的数据，1 ≤ n ≤ 1,000；

对于 100%的数据，1 ≤ n ≤ 100,000，0 ≤火柴高度≤ 2^31-1

[分析]

根据题意，两列火柴的总距离就等于$$\sum{a^2} + \sum{b^2} - 2 * \sum{a_i * b_i}$$
对于一个给定的序列，前两项均为常数，于是我们的目标就转为最大化所有对应的$a_i * b_i$之和。
这里需要一个结论：\[ \\ \quad \forall a_1 \leq a_2, b_1 \leq b_2\,\in\, \mathbb{R} ，
\quad a_1*b_1 + a_2*b_2 \geq a_1*b_2 + a_2 * b_1 \]
人生苦短，这个简单的结论我们就不证了~ 2333333（好吧其实直接把不等式两边减一下再因式分解就好……）
那么我们只需要把在对应序列中名次相同的两项匹配起来就好了……
于是我们就固定A不动(因为对A或B操作在对“两个同名次的项在两序列中的相对位置"的影响上是等价的)，把B中每个元素都与A中对应名次的元素对齐。不难发现，$$\forall i \le j \in [1, n]，若order_i \ge order_j，则i与j至少有一次交换.$$ 那么最少的交换次数就等于order序列中的逆序对个数>_< 又是经典模型，随便套个算法就搞定了= =

1 #include <cstdio>
2 #include <cctype>
3 #include <algorithm>
4 #include <cstring>
5 #include <queue>
6 #include <cmath>
7
8 #if defined DEBUG
9 FILE *in = fopen("temp", "r");
#define out stdout
#else
FILE *in = fopen("MatchNOIP2013.in", "r");
FILE *out = fopen("MatchNOIP2013.out", "w");
#endif

inline void getint(int &x){
char c = fgetc(in);
while(!isdigit(c))c = fgetc(in);
x = c - '0';
while(isdigit(c = fgetc(in)))x = x * 10 - '0' + c;
}
typedef long long LL;
using namespace std;
inline int lowbit(int x){return x & -x;}
/*=====================================*/
const int maxn = 100000 + 5, mod = 99999997;
int n, A[maxn], B[maxn], t1[maxn], t2[maxn], ans = 0;

inline bool tmpa(int a, int b){return A[a] < A[b];}
inline bool tmpb(int a, int b){return B[a] < B[b];}

//namespace Fenwick{
int Arr[maxn] = {0};
inline void insert(int x){
++Arr[x];
int i = x + lowbit(x);
while(i <= n){
++Arr[i];
i += lowbit(i);
}
}
inline int getS(int x){
int ans = Arr[x], i = x ^ lowbit(x);
while(i){
ans += Arr[i];
i ^= lowbit(i);
}
return ans;
}
//}//Fenwick
//using Fenwick::getS;
//using Fenwick::insert;

inline void work(){
getint(n);
int i, ord[maxn];
for(i = 1;i <= n;++i)
getint(A[i]), t1[i] = i;
for(i = 1;i <= n;++i)
getint(B[i]), t2[i] = i;
sort(t1 + 1, t1 + n + 1, tmpa);
sort(t2 + 1, t2 + n + 1, tmpb);
for(i = 1;i <= n;++i)
ord[t1[i]] = t2[i];
for(i = n;i;--i){
ans = (ans + getS(ord[i])) % mod;
insert(ord[i]);
}
fprintf(out, "%d\n", (ans + mod) % mod);
}

int main(){

work();

return 0;
}Fenwick Tree