[NOIP 2005]-- 篝火晚会
额~~,对这组题感兴趣的具体的解题报告可以戳戳这里:http://wenku.baidu.com/view/878beb64783e0912a2162aa7.html?qq-pf-to=pcqq.c2c
【问题描述】
佳佳刚进高中,在军训的时候,由于佳佳吃苦耐劳,很快得到了教官的赏识,成为了“小教官”。在军训结束的那天晚上,佳佳被命令组织同学们进行篝火晚会。一共有n个同学,编号从1到n。一开始,同学们按照1,2,……,n的顺序坐成一圈,而实际上每个人都有两个最希望相邻的同学。如何下命令调整同学的次序,形成新的一个圈,使之符合同学们的意愿,成为摆在佳佳面前的一大难题。
佳佳可向同学们下达命令,每一个命令的形式如下:
(b1, b2,... bm -1, bm)
这里m的值是由佳佳决定的,每次命令m的值都可以不同。这个命令的作用是移动编号是b1,b2,…… bm –1,bm的这m个同学的位置。要求b1换到b2的位置上,b2换到b3的位置上,……,要求bm换到b1的位置上。
执行每个命令都需要一些代价。我们假定如果一个命令要移动m个人的位置,那么这个命令的代价就是m。我们需要佳佳用最少的总代价实现同学们的意愿,你能帮助佳佳吗?
【输入文件】
输入文件fire.in的第一行是一个整数n(3 <= n <= 50000),表示一共有n个同学。其后n行每行包括两个不同的正整数,以一个空格隔开,分别表示编号是1的同学最希望相邻的两个同学的编号,编号是2的同学最希望相邻的两个同学的编号,……,编号是n的同学最希望相邻的两个同学的编号。
【输出文件】
输出文件fire.out包括一行,这一行只包含一个整数,为最小的总代价。如果无论怎么调整都不能符合每个同学的愿望,则输出-1。
【样例输入】
4
3 4
4 3
1 2
1 2
【样例输出】
2
【数据规模】
对于30%的数据,n <= 1000;
对于全部的数据,n <= 50000。
解题思路:
由原来的环可以得到一个目标序列枚举可能的2*n个不同的序列,比如1 2 3 4 5 6 7 8
1 8 2 7 5 3 6 4 10 9、8 2 7 5 3 6 4 10 9 1、2 7 5 3 6 4 10 9 1 8、
7 5 3 6 4 10 9 1 8 2、5 3 6 4 10 9 1 8 2 7、……及9 10 4 6 3 5 7 2 8 1
1 9 10 4 6 3 5 7 2 8、8 1 9 10 4 6 3 5 7 2、2 8 9 10 4 6 3 5 7、……。
将每个目标序列都与初始序列作比较,计算出目标序列与初始序列有多少个元素的位置已相同
(位置已相同的就不必移动,其它的都得移动,最少总代价就是得移动的只移动一次),
所以,问题就是求目标序列与初始序列最多有多少个(记为max)元素的位置已相同,则(n-max)就是所要求的最少总代价。
按上述算法编写程序,可以正确求出各组测试数据的结果,但当n=50000时,程序运行的时间约为40s~60s,果断超时~~~
那么如何快速求出目标环与初始环间最多有多少个元素的位置可能相同,就是ac本题的关键~~
关于求目标序列和原序列的最多等价元素个数的办法如下:
只要求出目标环各元素值与初始环对应位置的元素值之差的结果相同的最多是多少个就可以了
为使差值为正数,可用用 (a[i]-b[i]+n)%n来求差值(b[i]表示初始序列,实际就是12345……n)
同时注意因为是一个环,那么与方向无关,所以顺时针,逆时针我们都用上述方法求解一次
代码如下:
![](https://images.cnblogs.com/OutliningIndicators/ContractedBlock.gif)
1 #include <stdio.h> 2 #include <string.h> 3 int map[50005][2], ptr[50005], vis[50005], mpt[50005], n, flag; 4 5 void mergy(){ 6 int i, j, tmp1, tmp2; 7 mpt[0] = vis[mpt[0]] = 1; 8 for (i = 1; i < n; i++){ 9 tmp1 = map[mpt[i - 1] - 1][0], tmp2 = map[mpt[i - 1] - 1][1]; 10 if (!vis[tmp1]) mpt[i] = tmp1, vis[tmp1] = 1; 11 else if (!vis[tmp2]) mpt[i] = tmp2, vis[tmp2] = 1; 12 else{ 13 flag = 1; 14 return; 15 } 16 } 17 } 18 19 int solve(){ 20 int i, tmp, maxn = 0; 21 memset(ptr, 0, sizeof(ptr)); 22 for (i = 0; i < n; i++){ 23 tmp = (n + mpt[i] - i) % n; 24 ptr[tmp]++; 25 } 26 for (i = 0; i<n; i++) 27 maxn = maxn>ptr[i] ? maxn : ptr[i]; 28 memset(ptr, 0, sizeof(ptr)); 29 for (i = 0; i < n; i++){ 30 tmp = (n + mpt[n - 1 - i] - i) % n; 31 ptr[tmp]++; 32 } 33 for (i = 0; i<n; i++) 34 maxn = maxn>ptr[i] ? maxn : ptr[i]; 35 return n - maxn; 36 } 37 38 int main(){ 39 scanf("%d", &n); 40 for (int i = 0; i < n; i++) scanf("%d%d", &map[i][0], &map[i][1]); 41 mergy(); 42 printf("%d\n", flag ? -1 : solve()); 43 return 0; 44 }