【代码超详解】UVA 11212 Editing a book 编辑图书(迭代加深搜索 IDDFS ,150 ms)

一、题目描述

已知一些段落,用 1 ~ n 编号,2 ≤ n ≤ 9 。给出开始时各段落的顺序,要求通过单个剪贴板对段落剪切粘贴,使得顺序为最小排列。单次剪切可以剪切一个段落或多个连续段落。因为只有一个剪贴板,所以在一次粘贴操作前不能连续剪切两次。
输入不超过 20 个测试样例,每个样例有两行:第一行是一个整数 n ,接下来 1 行有 n 个数,它们均落在离散闭区间 [1, n] 内。用空格分隔两个相邻的数,代表最开始时的段落排序。输出一行,写出样例编号和最少需要做的剪切/粘贴的操作数。一次剪切和粘贴记 1 次操作。n = 0 时输入结束。具体格式见样例。

Sample Input

6
2 4 1 5 3 6
5
3 4 5 1 2
0

Sample Output

Case 1: 2
Case 2: 1

二、算法分析说明与代码编写指导

本算法和 AC 代码是根据博文 https://blog.csdn.net/weizhuwyzc000/article/details/47292669 改写的。
由于数据量不大(n ≤ 9),因此可以考虑直接暴力求解。我们可以简单粗暴地模拟整个尝试各种剪切、粘贴步骤的过程。
迭代加深搜索是一种 DFS 的衍生算法,可以看做是 DFS 和 BFS 的混合体之一。与 DFS 不同,使用 IDDFS 搜索时,会限制搜索的深度,当找不到解时才继续深入搜索。
以本体为例,如果直接采用 DFS ,为了找到操作步数最少的剪贴方法,可能需要在差不多搜索完全部的可行步骤时才能统计出最优解。既然要找步骤最少的解,那么这个解通常不会在很深层。也就是说我们没有必要深入深层进行全图搜索。
如果直接采用 BFS ,我们需要将所有的步骤都保存下来,然后再进行全图搜索。虽然搜索时可以较快地找到步数最少的解,但是空间需求很大,而且存图读图也需要不少的时间。
采用迭代加深搜索只需在 DFS 的基础上作一定量的改动,结合下面的算法详解可以看出其空间复杂度不比深搜差太多,但是比广搜小很多。当然,该方法也会重复搜索许多路径,所以一般在运行时间限制充裕时使用。
我们用 a 存储输入的段落序列,t 代表样例编号,dmax 代表最大搜索深度。
我们在启动搜索时需要依次传入 dmax = 0, 1, 2, …, n - 1 来实现 IDDFS 。
至于搜索函数的写法,我们先将搜索的终止条件写在前方,这样在每次进入深一层的搜索时都会检查是否需要终止搜索。

当达到搜索深度限制时,判断序列是否变为最小排列1、2、3、……、n即可。
我们用 h 表示序列中“不正确”的数字个数。具体解释如下:


到这里 h 和 sorted 函数就编写完毕。可见,如果某次输入的不正确数字个数较多,判定即使用光全部搜索层数都无法纠正,那么这次尝试在该深度下第一次由 main 函数调用搜索时就会被这个 if 语句拦住,直接结束搜索,于是层数限制增加 1 层,重新搜索。
当尝试过程中发现即使达到搜索深度限制时仍然无法找出解,就返回 false ,函数返回到浅一层,在浅一个深度更换最后一次步骤,继续搜索。

我们用 a 表示输入的段落序列。我们用一个离散闭区间 [i, j] 来表示需要剪切并粘贴的连续段落,未剪切的段落用字符数组 b 保存,a0 用于在尝试接下来的剪切粘贴(记为动作 A)前先将改动前的字符串保存下来,当深层尝试失败后,就要还原到动作 A 之前作新的尝试。
l1 代表未剪切的段落的长度,我们将其复制到 b 中。然后我们模拟将剪切下来的段落在剩下的字符中的各个位置(从左到右)都尝试粘贴。在一个位置尝试粘贴完毕后,就进行深一层的尝试(在未达到深度限制时)。如果在某一层能够成功将段落编号按最小排列排好,就返回 true ,并结束全部的递归函数。由于深度限制从 0 慢慢放开到 n-1,所以一旦找到第一个解,肯定是步骤最少的解。记住:搜索深度指的是已经进行的剪切粘贴的操作次数
搜索函数的传入值 d 代表已经进行的搜索层数,开始调用时,需要传 0 ,因为还未搜索过。
搜索最大深度是 n-1 的原因是最多只要(n-1)步就能将段落序列整理成1、2、……、n:
第一步,将 1 提到全部文段之前;第二步、将 2 提到 1 之后;第三步、将 3 提到 2 之后;……第(n-1)步:将 n 提到 n-1 之后。
已经找到结果后,直接跳出 main 函数里调用搜索函数的循环,输出结果。

三、AC 代码(150 ms)

#include<cstdio>
#include<algorithm>
#pragma warning(disable:4996)
using namespace std;
unsigned n, n0, a[9], t = 0, dmax;
inline unsigned h() {
	unsigned h0 = 0; for (unsigned i = 0; i < n0; ++i)if (a[i] + 1 != a[i + 1])++h0;
	if (a[n0] != n)++h0;
	return h0;
}
inline bool sorted() {
	for (unsigned i = 0; i < n0; ++i)if (a[i + 1] - a[i] != 1)return false;
	return true;
}
inline bool Enumerate(const unsigned& d, const unsigned& dmax) {
	if (d == dmax)return sorted();
	if (3 * (dmax - d) < h())return false;
	unsigned b[9], a0[9]; copy(a, a + 9, a0);
	for (unsigned i = 0; i < n; ++i) {
		for (unsigned j = i; j < n; ++j) {
			unsigned l1 = 0;
			for (unsigned k = 0; k < i; ++k) { b[l1] = a[k]; ++l1; }
			for (unsigned k = j + 1; k < n; ++k) { b[l1] = a[k]; ++l1; }
			for (unsigned k = 0; k < l1; ++k) {
				unsigned l2 = 0;
				for (unsigned p = 0; p < k; ++p) { a[l2] = b[p]; ++l2; }
				for (unsigned p = i; p <= j; ++p) { a[l2] = a0[p]; ++l2; }
				for (unsigned p = k; p < l1; ++p) { a[l2] = b[p]; ++l2; }
				if (Enumerate(d + 1, dmax) == true)return true;
				copy(a0, a0 + 9, a);
			}
		}
	}
	return false;
}
int main() {
	for (;;) {
		scanf("%u", &n); if (n == 0)break;
		n0 = n - 1; for (unsigned i = 0; i < n; ++i)scanf("%u", &a[i]);
		for (dmax = 0; dmax < n0; ++dmax) { if (Enumerate(0, dmax) == true)break; }
		printf("Case %u: %u\n", ++t, dmax);
	}
	return 0;
}
  • 0
    点赞
  • 1
    收藏
    觉得还不错? 一键收藏
  • 0
    评论

“相关推荐”对你有帮助么?

  • 非常没帮助
  • 没帮助
  • 一般
  • 有帮助
  • 非常有帮助
提交
评论
添加红包

请填写红包祝福语或标题

红包个数最小为10个

红包金额最低5元

当前余额3.43前往充值 >
需支付:10.00
成就一亿技术人!
领取后你会自动成为博主和红包主的粉丝 规则
hope_wisdom
发出的红包
实付
使用余额支付
点击重新获取
扫码支付
钱包余额 0

抵扣说明:

1.余额是钱包充值的虚拟货币,按照1:1的比例进行支付金额的抵扣。
2.余额无法直接购买下载,可以购买VIP、付费专栏及课程。

余额充值