【逆序对相关/数学】【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 这题其实主要难在想方法吧我觉得?学长提点了下说用贪心之后就大概明白了,感觉没有很难 但是离散化这里还是挺有趣的,因为并不是能很熟练地掌握离 ...
随机推荐
- 对网页进行截图(selenium)
import os def insert_img(driver,file_name): #获取当前路径,并转换为字符串 base_dir=str(os.path.dirname(__file__)) ...
- Error -26377: No match found for the requested parameter
Error -26377: No match found for the requested parameter
- PyMySQL连接MySQL数据库
首先, 添加PyMySQL模块: 代码: import pymysql db = pymysql.connect(host="localhost", user="root ...
- 适配iPhoneX、iPhoneXs、iPhoneXs Max、iPhoneXr 屏幕尺寸及安全区域
此篇文章是对上一篇文章(http://www.ifiero.com/index.php/archives/611)的进一步补充,主要说明如何适配Apple的最新三款手机iPhoneXs.iPhoneX ...
- flask源码走读
Flask-Origin 源码版本 一直想好好理一下flask的实现,这个项目有Flask 0.1版本源码并加了注解,挺清晰明了的,我在其基础上完成了对Werkzeug的理解部分,大家如果想深入学习的 ...
- 关于java中的“error: bad operand types for binary operator ”
今天做这个题目的时候(142. O(1) Check Power of 2),遇到一个错误“ bad operand types for binary operator '&' ”. 先贴一下 ...
- python numpy数据相减
numpy数据相减,a和b两者shape要一样,然后是对应的位置相减.要不然,a的shape可以是(1,m),注意m要等于b的列数. import numpy as np a = [ [0, 1, 2 ...
- github项目切换远程https到ssh通道
github 切换远程https到ssh通道 github 每个仓库有两类地址:https和ssh通道. https通道获取代码方便,提交过程中每次都需要输入用户名和密码. ssh通道需要提前配置号s ...
- 七:Web Application Proxy
yarn自带了web接口,默认是和RM一起的(8088端口).但是为了减少从web接口受到的攻击,可以把Web接口单独放在别的机器上. 设置下web代理就行了 Configurations Confi ...
- [C++] Class (part 1)
The fundamental ideas behind classes are data abstraction and encapsulation. Data abstraction is a p ...