【逆序对相关/数学】【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 这题其实主要难在想方法吧我觉得?学长提点了下说用贪心之后就大概明白了,感觉没有很难 但是离散化这里还是挺有趣的,因为并不是能很熟练地掌握离 ...
随机推荐
- MySQL☞in语句
in语句: 1)列名 in(数值1,数值2,数值3…):求出满足该列的多个列值 格式: select 列名1,列名2 from 表名 where 列名 in (数值1,数值2,数值3...) 如下图 ...
- 创建一个socket服务实时统计在线人数
主要是两个文件,一个是后端文件,一个是前端文件: 后端文件:有人登录了,就通知所有的正在访问的页面,把总人数+1:反之-1: 前端文件:有人登录了,通知后端,页面关闭了,通知后端,同时接收后端派发来的 ...
- MySQL三方面优化
第一方面:30种mysql优化sql语句查询的方法1.对查询进行优化,应尽量避免全表扫描,首先应考虑在 where 及 order by 涉及的列上建立索引. 2.应尽量避免在 where 子句中使用 ...
- LeetCode 98——验证二叉搜索树
1. 题目 2. 解答 2.1. 方法一 我们初始化根节点的范围为长整形数据的最小最大值 \([LONG\_MIN,LONG\_MAX]\),则其左子节点的取值范围为 \([LONG\_MIN,根节点 ...
- 基于AdaBoost算法——世纪晟结合Haar-like特征训练人脸检测识别
AdaBoost 算法是一种快速人脸检测算法,它将根据弱学习的反馈,适应性地调整假设的错误率,使在效率不降低的情况下,检测正确率得到了很大的提高. 系统在技术上的三个贡献: 1.用简单的Haa ...
- STM32F4 编程手册学习1_编程模型
STM32F4 programming manual_1 1. 处理器模式与特权等级 处理器模式分为以下两种: 线程模式: 用来执行应用软件: 处理器从reset出来时,进入线程模式: CONTROL ...
- 软工实践Alpha冲刺(2/10)
队名:我头发呢队 组长博客 作业博客 杰(组长) 过去两天完成了哪些任务 查看了讯飞开放平台的sdk 查阅了Google Material Design 2的官方文档 接下来的计划 继续打磨UI界面: ...
- [solution]xdebug正确配置,但不显示错误信息
一开始以为是配置问题,其实不是,折腾了好久,貌似中文网页很少有人提到这事,更别提解决之道! 最好还是用英文关键词google之:得如下网页 https://bugs.launchpad.net/ubu ...
- DDB与DIB
DB与DIB的区别是什么?觉得书上介绍的有点抽象.不容易理解.他们两者之间的区别的“物理意义” [“现实意义”]——姑且这么叫吧,呵呵!被这个问题困扰了很久,所以今天决定好好查资料总结一下,把它彻底搞 ...
- ORA-00933 SQL命令未正确结束 INSERT INTO ... SELECT
最近在修改数据库存储过程时,出现了一个ORA-00933错误, 执行的是 INSERT INTO...SELECT 语句,具体语句如下: INSERT INTO BASP_DX.QLR@GT(BDCD ...