迭代加深法, 可以解决一些看起来更适合用BFS和回溯法的题,但是BFS和回溯法浪费空间,所以迭代加深法更能解决问题。
分析:
由于编辑的特殊性,只要没达到升序这个条件就可以永远不停的搜下去,最特别的还是那些的错误的路径,也就是说假如一开始就选择了错误的剪切方式,之后这条错误的路会一直往后延伸。。。
所以采用IDA*解法,IDA*相比与普通的dfs或bfs,多了一个估值函数,比如我从这一层到第n层,每层最多能干多少活,假如剩下的活到第n层怎么都做不完,就不要往下搜了。像这道题,每层最多有三个数字的后缀改变,也就是说最多能使三个数字的后继正确。
证明:设当前有四个数字 a b c d ,如果把b移到c后面,则改变了a、b、c三个数的后继
所以每次搜索前,当h>(maxd-d)*3 时剪枝,这里maxd-d即剩余搜索次数;
搜索的核心就是截取 [i, j] 插入剩余序列的第k个数字前
关于这IDA*的减枝策略,每遍历一个深度的时候,进行判断:
当前局面的估价函数值 + 当前深度 > 预定义最大搜索深度
的时候进行减枝。
就能这道题而言,假如我们定义一个数字是不是位置正确:这个数 x 是否 等于 这个数 后面的 数 y - 1,也就是
x =? y - 1,如果等于,说明这个数位置正确,如果不等于,说明这个数位置错误,位于最后一个位置的数的时候,判断他是不是等于n,比如:4,5,6,1,2,3 这个序列存在2个不正确位置数,分别是6(后面是1)和3(3不等于6)。
下面进行这道题的减枝的分析(也叫做启发函数),当你改变一个区间的位置,你会改变3个数的位置的正确性
比如 1,2,3,4,5,6.序列,你把2,3移动到6后面,那么1的后面变成了5, 而 6的后面编程了2,而3的后面变成 空了,所以每次移动一个区间,最多可以改变3个数的正确性,也就是说,对于这道题如果遍历到了一个深度, (还能遍历的深度 - 当前深度) *3 < 不正确数字的个数,那么就没有必要继续遍历了,因为往后你就是全把这些数字该对了也无法达到理想状态,即使再遍历也是做无用功。
知道这个之后时间复杂度的问题就得到解决了,下面我们只需要每次枚举该步的所有移动就可以了。
移动的话,实际就是2个相邻区间的交换,比如A B C D(字母代表区间),将A移动到C后面,也可以看成A 和(BC)互换,所以实质就是相邻区间的互换。
#include
#include
using namespace std;
const int maxn = 9;
int n, a[maxn];
bool is_sorted() {
for(int i = 0; i < n-1; i++)
if(a[i] >= a[i+1]) return false;
return true;
}// 判断是否到达目标状态
int h() {
int cnt = 0;
for(int i = 0; i < n-1; i++)
if(a[i]+1 != a[i+1]) cnt++;
if(a[n-1] != n) cnt++;
return cnt;
} //统计需要变更的总数
bool dfs(int d, int maxd) {
if(d*3 + h() > maxd*3) return false; //估价函数
if(is_sorted()) return true; //递归的边界条件
int b[maxn], olda[maxn];
memcpy(olda, a, sizeof(a)); //存放原图
for(int i = 0; i < n; i++)
for(int j = i; j < n; j++) {
int cnt = 0;
for(int k = 0; k < n; k++)
if(k < i || k > j) b[cnt++] = a[k];
for(int k = 0; k <= cnt; k++) {
int cnt2 = 0;
for(int p = 0; p < k; p++) a[cnt2++] = b[p]; //插入已经拍好序的
for(int p = i; p <= j; p++) a[cnt2++] = olda[p]; //插入
for(int p = k; p < cnt; p++) a[cnt2++] = b[p]; //插入终止状态
if(dfs(d+1, maxd)) return true;
memcpy(a, olda, sizeof(a));
}
}
return false;
}
int solve() {
if(is_sorted()) return 0;
int max_ans = 8;
for(int maxd = 1; maxd < max_ans; maxd++)
if(dfs(0, maxd)) return maxd;
return max_ans;
}
int main() {
int kase = 0;
while(scanf("%d", &n) == 1 && n) {
for(int i = 0; i < n; i++) scanf("%d", &a[i]);
printf("Case %d: %d\n", ++kase, solve());
}
return 0;
}