【逆序对相关/数学】【P1966】【NOIP2013D1T2】 火柴排队
Description
涵涵有两盒火柴,每盒装有 $n$ 根火柴,每根火柴都有一个高度。 现在将每盒中的火柴各自排成一列, 同一列火柴的高度互不相同, 两列火柴之间的距离定义为:$ \sum (a_i-b_i)^2$
其中$ a_i$ 表示第一列火柴中第$ i $个火柴的高度,$b_i$表示第二列火柴中第 $i$ 个火柴的高度。
每列火柴中相邻两根火柴的位置都可以交换,请你通过交换使得两列火柴之间的距离最小。请问得到这个最小的距离,最少需要交换多少次?如果这个数字太大,请输出这个最小交换次数对 $99,999,997 $取模的结果。
## Input
共三行,第一行包含一个整数$ n$,表示每盒中火柴的数目。
第二行有$ n $个整数,每两个整数之间用一个空格隔开,表示第一列火柴的高度。
第三行有 $n$ 个整数,每两个整数之间用一个空格隔开,表示第二列火柴的高度。
## Output
一个整数,表示最少交换次数对 $99,999,997$ 取模的结果。
## Sample Input
```
4
2 3 1 4
3 2 1 4
```
## Sample Output
```
1
```
## Hint
对于 $100\%$的数据,$1 ≤ n ≤ 100,000,0 ≤$火柴高度$≤ maxlongint$
## Solution
和[NOIP2011D2T1聪明的质检员](https://www.cnblogs.com/yifusuyi/p/9531622.html)一样,既然题目给了一个表达式并且瞪眼看不出这个式子有啥可做的,那么就化简这个式子。
把式子平方展开,则$ \sum (a_i-b_i)2~=~\sum~a_i2+\sumb_i^2+2\times\suma_i*b_i$。
显然\(\sum~a_i^2~+~\sum~b_i^2\)已经固定了,所以我们要做的是最小化\(\sum~a_i*b_i\)
根据排序不等式,上式最小当且仅当\(\forall i\),a,b第\(i\)大的数在同一位置。
那么我们按照\(a_i\)的顺序给\(b_i\)排序。即新建一个序列\(C\),将a,b离散化后使得\(C\)满足\(C_{a_i}=b_i\)。它的数学意义是\(a\)中第i大的数对应位置上\(b\)的大小。
显然,达到合法的情况应该满足\(\forall i,C_i=i\)。
显然一起移动A,B和只移动B是等价的,那么本题就变为交换\(C\)中相邻元素用最少的步数满足\(C\)中元素按照升序排序。
咦?这不就是逆序对嘛?
其实我也是才知道这就是逆序对的意义
证明可以使用数学归纳法。
先将最大的元素移动到后面去。显然元素列他后面都是都是比他小的,所以把他移动到最后的步数就等于关于它的逆序对数。
他到最后去以后,下面的问题等价于将剩下n-1个数升序排列。
由于最大的数和后面每个数都进行了交换,所以除第n个数以外整个数列的相对的大小位置是不变的
这样就使用数学归纳法证明了上述问题的最小步数就是逆序对的个数。
Code
#include<cstdio>
#include<algorithm>
#define rg register
#define ci const int
#define cl const long long int
typedef long long int ll;
namespace IO {
char buf[50];
}
template<typename T>
inline void qr(T &x) {
char ch=getchar(),lst=' ';
while(ch>'9'||ch<'0') lst=ch,ch=getchar();
while(ch>='0'&&ch<='9') x=(x<<1)+(x<<3)+(ch^48),ch=getchar();
if (lst=='-') x=-x;
}
template<typename T>
inline void write(T x,const char aft,const bool pt) {
if(x<0) {putchar('-');x=-x;}
int top=0;
do {
IO::buf[++top]=x%10+'0';
x/=10;
} while(x);
while(top) putchar(IO::buf[top--]);
if(pt) putchar(aft);
}
template <typename T>
inline T mmax(const T a,const T b) {if(a>b) return a;return b;}
template <typename T>
inline T mmin(const T a,const T b) {if(a<b) return a;return b;}
template <typename T>
inline T mabs(const T a) {if(a<0) return -a;return a;}
template <typename T>
inline void mswap(T &a,T &b) {T temp=a;a=b;b=temp;}
const int maxn = 100010;
const int MOD = 99999997;
int n;
int a[maxn],b[maxn],c[maxn],d[maxn],frog[maxn],temp[maxn];
ll ans;
inline bool cmp1(const int &_a,const int &_b) {
return a[_a]<a[_b];
}
inline bool cmp2(const int &_a,const int &_b) {
return b[_a]<b[_b];
}
void ms(ci,ci);
int main() {
qr(n);
for(rg int i=1;i<=n;++i) qr(a[i]);
for(rg int i=1;i<=n;++i) qr(b[i]);
for(rg int i=1;i<=n;++i) c[i]=d[i]=i;
std::sort(c+1,c+1+n,cmp1);std::sort(d+1,d+1+n,cmp2);
for(rg int i=1;i<=n;++i) frog[c[i]]=d[i];
ms(1,n);
write(ans,'\n',true);
return 0;
}
void ms(ci l,ci r) {
if(l>=r) return;
int mid=(l+r)>>1;
ms(l,mid);ms(mid+1,r);
int i=l,j=mid+1,p=l;
while((i<=mid) && (j<=r)) {
if(frog[i]<frog[j]) temp[p++]=frog[i++];
else temp[p++]=frog[j++],ans=(ans+mid-i+1)%MOD;
}
while(i<=mid) temp[p++]=frog[i++];
while(j<=r) temp[p++]=frog[j++];
for(i=l;i<=r;++i) frog[i]=temp[i];
}
Summary
1、遇见棘手的表达式,先化简再说
2、将一个序列交换相邻元素用最小步数使得序列按照升序排列的答案就是序列中的逆序对个数
3、代码中给出了离散化的方法以及归并排序求逆序对的方法。
【逆序对相关/数学】【P1966】【NOIP2013D1T2】 火柴排队的更多相关文章
- 洛谷P1966 【火柴排队】
题解 P1966 [火柴排队] 说明: 在数学中有个公式: (a1-b1)^2+(a2-b2)^2<(a2-b1)^2+(a1-b2)^2 (你可以自己试着证一下) 两列火柴对应的两根火柴在各列 ...
- 【题解】洛谷P1966 [NOIP2013TG] 火柴排队(树状数组+逆序对)
次元传送门:洛谷P1966 思路 显然在两排中 每排第i小的分别对应就可取得最小值(对此不给予证明懒) 所以我们只在意两排的火柴是第几根 高度只需要用来进行排序(先把两个序列改成有序的方便离散化) 因 ...
- 【洛谷P1966】火柴排队
火柴排队 题目链接 ∑(ai−bi)^2=∑ai^2-2*∑ai*bi+∑bi^2 显然∑ai^2+∑bi^2是不变的,我们要让 2*∑ai*bi最大,才能使原式最小 然后我们一眼就可以看出来, ...
- [hdu4911]逆序对相关
思路:由于只能交换相邻的数,所以每次最多减小1个逆序对(且如果存在逆序对那么肯定可以减小1个)!于是乎..就是统计逆序对的裸题了.树状数组或归并都行. #pragma comment(linker, ...
- [NOIP2013提高&洛谷P1966]火柴排队 题解(树状数组求逆序对)
[NOIP2013提高&洛谷P1966]火柴排队 Description 涵涵有两盒火柴,每盒装有 n 根火柴,每根火柴都有一个高度. 现在将每盒中的火柴各自排成一列, 同一列火柴的高度互不相 ...
- luogu P1966 火柴排队 (逆序对)
luogu P1966 火柴排队 题目链接:https://www.luogu.org/problemnew/show/P1966 显然贪心的想,排名一样的数相减是最优的. 证明也很简单. 此处就不证 ...
- P1966 火柴排队——逆序对(归并,树状数组)
P1966 火柴排队 很好的逆序对板子题: 求的是(x1-x2)*(x1-x2)的最小值: x1*x1+x2*x2-2*x1*x2 让x1*x2最大即可: 可以证明将b,c数组排序后,一一对应的状态是 ...
- P1966 火柴排队(逆序对)
P1966 火柴排队 题目描述 涵涵有两盒火柴,每盒装有 n 根火柴,每根火柴都有一个高度. 现在将每盒中的火柴各自排成一列, 同一列火柴的高度互不相同, 两列火柴之间的距离定义为: ∑(ai-bi) ...
- 洛谷P1966 火柴排队 贪心+离散化+逆序对(待补充QAQ
正解: 贪心+离散化+逆序对 解题报告: 链接在这儿呢quq 这题其实主要难在想方法吧我觉得?学长提点了下说用贪心之后就大概明白了,感觉没有很难 但是离散化这里还是挺有趣的,因为并不是能很熟练地掌握离 ...
随机推荐
- Python里//与/的区别?
1.Python里面//的作用是除法取整,也就是直接取整数部分 例如:5//6=0; 56//3=18 2.而/的作用是直接进行常规的除法运算 例如:56/8=7 程序运算实例如下:
- TPO-16 C1 Reserve the room for a rehearsal
TPO-16 C1 Reserve the room for a rehearsal 第 1 段 1.Listen to a conversation between a Student and a ...
- (C#)设计模式之装饰模式
1.装饰模式 动态的给一个对象添加一些额外的职责,就添加功能来说,装饰模式比生成子类更加灵活.*装饰模式是为已有功能动态添加更多功能的一种方式.*装饰模式将原有类中的核心职责与装饰功能分离.简化了原有 ...
- leetcode-生成括号(回溯算法)
转载出处:https://blog.csdn.net/yanerhao/article/details/68561290 生成括号 给出 n 代表生成括号的对数,请你写出一个函数,使其能够生 ...
- 从零开始的Python学习Episode 6——字符串操作
字符串操作 一.输出重复字符串 print('smile'*6) #输出6个smile 二.通过引索输出部分字符串 print('smile'[1:]) print('smile'[1:3]) #输出 ...
- [HNOI2018]转盘
[HNOI2018]转盘 给你一个 \(n\) 元环, 你可以在 \(0\) 时刻从任意一个位置出发, 每一秒可以选择往后或者留在原地每个点有个参数 \(T_i\) , 当你走到 \(i\) 的时间 ...
- 剑指offer-二叉搜索树的后序遍历序列23
题目描述 输入一个整数数组,判断该数组是不是某二叉搜索树的后序遍历的结果.如果是则输出Yes,否则输出No.假设输入的数组的任意两个数字都互不相同. class Solution: def Verif ...
- 线性代数之——微分方程和 exp(At)
本节的核心是将常系数微分方程转化为线性代数问题. \[\frac{du}{dt}=\lambda u \quad 的解为 \quad u(t) = Ce^{\lambda t}\] 代入 \(t=0\ ...
- CSS动画@-webkit-keyframes
@-webkit-keyframes:以百分比来规定改变发生的时间,或者通过关键词 "from" 和 "to",等价于 0% 和 100%.0% 是动画的开始时 ...
- HDU 2487 Ugly Windows(暴力)(2008 Asia Regional Beijing)
Description Sheryl works for a software company in the country of Brada. Her job is to develop a Win ...