题解【NOIP2013_火柴排序】

题解_NOIP2013_火柴排序

题意：

有两个序列a,b，每个序列的元素不同。任意一个序列的任意两个相邻的元素可以互换，要求对两个序列进行操作，使得∑(ai−bi)2\sum(a_i-b_i)^2∑(ai−bi)2最小，求最小操作次数。

题解：

对于这两个序列而言，要求我们最小化的是∑(ai−bi)2\sum(a_i-b_i)^2∑(ai−bi)2
展开这个式子，得到∑(ai2−2aibi+bi2)\sum(a_i^2-2a_ib_i+b_i^2)∑(ai2−2aibi+bi2) 也就是∑(ai2+bi2)−∑(2aibi)\sum(a_i^2+b_i^2)-\sum(2a_ib_i)∑(ai2+bi2)−∑(2aibi) 因为ai,bia_i,b_iai,bi都是定值，所以我们只需要最大化∑(2aibi)\sum(2a_ib_i)∑(2aibi)即可。

对于两个序列而言，全部升序/降序时总比乱序时的∑(2aibi)\sum(2a_ib_i)∑(2aibi)大。这一点很容易证明可以考虑转化为面积或者通过数学归纳法进行证明。

另外很显然的一点就是，对于序列a,b，对他们的对应元素进行同样的操作，也就是“一起进行某种操作”对答案没有影响。

那么，问题就转化为，对于序列a,b，求对b最少的操作次数使得对b操作后产生的b’序列和a“一起”进行某些操作后使得a变为升序后，a,b’的元素对应的大小位置相等（大小位置就是对于某个序列排完序以后某个元素的位置）。

那么，我们先对a，b序列进行一步奇妙的离散化。把序列中的元素映射为对应元素的下标，设a映射为c，b映射为d。然后，对吧c以a为基准排序，再把d以b为基准排序。得到新的序列c’,d’。此时这两个序列保存的是把a和b分别变成升序时，原来元素的新位置。那么，根据前面所说的，a和b“一起”进行某些操作时，对答案没有影响，我们已经成功的把问题转化为：只能交换序列中的相邻元素，使得d’=c’的最小交换次数。但是我们仍然不能很容易的进行求解。

那么，我们把c排序，再把d’以c’为基准排序。这时候得到的新序列c’’,d’‘相当于把c’,d’ “一起”进行了某些操作使得c’升序时的d’。到了这一步我们就会求只能交换序列中的相邻元素，使得d’=c’的最小交换次数了。

根据线性代数老师所言：

只交换相邻的两数让他们变成升序：交换次数为逆序对数
证明：考虑一个排列，先将其最大的数移至序列尾，次数为最大数与它后面的数组成的逆序对数。此时最后一位已定，再考虑前n-1个数的排列。最后整个序列的操作次数即为逆序对数。

代码：

#include <iostream>
#include <cstdio>
#include <algorithm>
using namespace std;
const int MAXN = 1e5+7,P=1e8-3;
template <typename _TP>
inline _TP read(_TP &X){char ch=0;int w;X=0;while(!isdigit(ch)){w=ch=='-';ch=getchar();}while(isdigit(ch)){X=(X<<1)+(X<<3)+(ch^48);ch=getchar();}X=w?-X:X;return X;
}
long long a[MAXN],b[MAXN],t[MAXN],pos[MAXN],w[MAXN];
inline bool cmpA(const int x,const int y){return a[x]<a[y];
}
inline bool cmpB(const int x,const int y){return b[x]<b[y];
}
long long cnt=0;
void merge_sort(const int l,const int r){if(r-l>1){int m=l+(r-l)/2;int p=l,q=m,i=l;merge_sort(l,m);merge_sort(m,r);while(p<m || q<r){if(q>=r || (p<m && w[p]<=w[q])){t[i++]=w[p++];}else {t[i++]=w[q++];cnt+=m-p;}}for(int i=l;i<r;i++){w[i]=t[i];}}
}
int c[MAXN],d[MAXN];
int main(){long long n;read(n);for(int i=1;i<=n;i++){read(a[i]);c[i]=i;}for(int i=1;i<=n;i++){read(b[i]);d[i]=i;}sort(c+1,c+n+1,cmpA);sort(d+1,d+n+1,cmpB);for(int i=1;i<=n;i++){w[c[i]]=d[i];}merge_sort(1,n+1);cout<<cnt%(P);return 0;
}
/*
关于序列操作：1.只交换相邻的两数让他们变成升序：交换次数为逆序对数原理需要咨询线性代数老师:考虑一个排列，先将其最大的数移至序列尾，次数为最大数与它后面的数组成的逆序对数此时最后一位已定，再考虑前n-1个数的排列2.交换任意的两数是使得序列升序交换次数为n-循环个数定义一个循环为3在5的位置，5在8的位置，8在3的位置
*/

题解【NOIP2013_火柴排序】— 问题与转化相关推荐

【题解】 [HEOI2016]排序题解（二分答案，线段树）
题目描述在2016年,佳媛姐姐喜欢上了数字序列.因而他经常研究关于序列的一些奇奇怪怪的问题,现在他在研究一个难题,需要你来帮助他.这个难题是这样子的:给出一个1到n的全排列,现在对这个全排列序列进行 ...
洛谷 P3244 / loj 2115 [HNOI2015] 落忆枫音题解【拓扑排序】【组合】【逆元】
组合计数的一道好题.什么非主流题目题目背景 (背景冗长请到题目页面查看) 题目描述不妨假设枫叶上有 \(n\) 个穴位,穴位的编号为 \(1\sim n\).有若干条有向的脉络连接着这些穴位. ...
ios 数组中的字典排序_题解 | 26.删除排序数组中的重复项
假设我们有一个传感器,会不停地向 HQ 传输已经归类好数据,数据的格式全部是数字(种类编号),且已经按照从小到大的顺序排列,大概是 [1,1,4,5,7,9] 这样的,但是我们的需求在于,希望知道有多 ...
【题解】（排序） —— POJ 0803:DNA排序
总时间限制: 1000ms 内存限制: 65536kB 描述现在有一些长度相等的DNA串(只由ACGT四个字母组成),请将它们按照逆序对的数量多少排序. 逆序对指的是字符串A中的两个字符A[i].A ...
【题解】（排序） —— POJ 0810:距离排序
总时间限制: 1000ms 内存限制: 65536kB 描述给出三维空间中的n个点(不超过10个),求出n个点两两之间的距离,并按距离由大到小依次输出两个点的坐标及它们之间的距离. 输入输入包括两 ...
AcWing 1191. 家谱树题解（拓扑排序）
AcWing 1191. 家谱树拓扑排序:建立一个队列,遍历所有点,将度数为1的点全部加入队列中,之后队内元素依次出队,对出队的点遍历所有临边,将这些临点的入度–,如果临点的入度为0,则加入队列中, ...
P3387 【模板】缩点洛谷 java题解连通图+拓扑排序
传送门: P3387 [模板]缩点 - 洛谷 | 计算机科学教育新生态 (luogu.com.cn)https://www.luogu.com.cn/problem/P3387 问题分析首先为了节省 ...
【题解】（排序） —— POJ 0811:牛的选举
总时间限制: 1000ms 内存限制: 65536kB 描述现在有N(1<=N<=50000)头牛在选举它们的总统,选举包括两轮:第一轮投票选举出票数最多的K(1<=K<=N ...
QLU_ACM 2021 专题训练（一）题解 [暴力、排序、贪心、二分]
欢迎大家来做QLU_ACM寒假组织的专题训练
题解：单词排序（水）
题目描述小红学会了很多英文单词,妈妈为了帮小红加强记忆,拿出纸.笔,把 N 个单词写在纸上的一行里,小红看了几秒钟后,将这张纸扣在桌子上.妈妈问小红:"你能否将这 N 个单词按照字典排列的 ...

题解【NOIP2013_火柴排序】— 问题与转化

题解_NOIP2013_火柴排序

题意：

题解：

代码：

题解【NOIP2013_火柴排序】— 问题与转化相关推荐

最新文章

热门文章