【ACWing】180. 排书

题目地址：

https://www.acwing.com/problem/content/182/

给定$n$本书，编号为$1\sim n$。在初始状态下，书是任意排列的。在每一次操作中，可以抽取其中连续的一段，再把这段插入到其他某个位置。我们的目标状态是把书按照$1\sim n$的顺序依次排列。求最少需要多少次操作。

输入格式：
第一行包含整数$T$，表示共有$T$组测试数据。每组数据包含两行，第一行为整数$n$，表示书的数量。第二行为$n$个整数，表示$1\sim n$的一种任意排列。同行数之间用空格隔开。

输出格式：
每组数据输出一个最少操作次数。如果最少操作次数大于或等于$5$次，则输出5 or more。每个结果占一行。

数据范围：
$1\le n\le 15$

思路是迭代加深。可以模仿A*算法的办法，设一个启发函数$h$，如果当前状态是$s$，那么$h(s)$就是从$s$这个状态走到终点状态需要的步数的某个下界。迭代加深的英文叫做Iterative DFS，所以带A*算法优化的迭代加深就叫做IDA*。我们考虑一下$h$的选取。对于一个$1\sim n$的序列，如果某个位置$a_i+1\ne a_{i+1}$，那么我们称这个位置产生了一个错配。显然升序排列的错配数是$0$。而做一次操作，错配数最多减$3$（例如将$a_{i\sim j}$这一段挪到$a_k$后面，那么$a_{i-1},a_j,a_k$这三个数的错配有可能恰好全部消失，当然也有可能本来没有错配，这么一换产生了错配。但无论怎样，每次减去$3$个错配的理想操作，是操作数的一个下界），设$s$的错配数是$m(s)$，所以可以定义$h(s)=\lceil \frac{m(s)}{3}\rceil =\lfloor \frac{m(s)+2}{3}\rfloor$。DFS的过程即是每次都枚举所有的操作可能的情况。当当前状态$s$的深度加上$h(s)$已经大于了最大深度了，就可以直接返回false了，因为$s$这棵子树上一定没有深度小于等于最大深度的解。

注意这里还有一个可以优化的地方。在枚举区间该挪到哪里的时候，可以只枚举其向后挪的可能性而不必枚举向前挪，因为这可以由挪到的位置到区间左端点所成区间向后挪等价得到。

代码如下：

#include <iostream>
#include <cstring>
using namespace std;

const int N = 18;
int n, a[N];
// w[i]记录的是DFS到第i层的时候，预存的”现场“，以便回溯回来的时候”恢复现场“
int w[5][N];

// 启发函数
int f() {
    int cnt = 0;
    for (int i = 0; i < n - 1; i++)
        if (a[i + 1] != a[i] + 1) cnt++;

    return (cnt + 2) / 3;
}

bool dfs(int depth, int max_depth) {
    if (depth + f() > max_depth) return false;
	
	// 易知如果错配数是0，对应的就是终止状态
    if (f() == 0) return true;

	// 枚举操作的区间长度
    for (int len = 1; len <= n; len++)
    	// 枚举左端点
        for (int l = 0; l + len - 1 < n; l++) {
        	// 算出右端点
            int r = l + len - 1;
            // 枚举挪到哪个位置，只枚举向后挪
            for (int k = r + 1; k < n; k++) {
            	// 记录现场
                memcpy(w[depth], a, sizeof a);
                // 开始挪，算一下挪完得到的状态
                int y = l;
                for (int x = r + 1; x <= k; x++)
                    a[y++] = w[depth][x];
                for (int x = l; x <= r; x++)
                    a[y++] = w[depth][x];
				
				// 看一下本次操作能否得到深度不超过max_depth的解，如果能则返回true
                if (dfs(depth + 1, max_depth)) 
                    return true;
				
				// 恢复现场
                memcpy(a, w[depth], sizeof a);
            }
        }

    return false;
}

int main() {
    int T;
    cin >> T;
    while (T--) {
        cin >> n;
        for (int i = 0; i < n; i++) cin >> a[i];

        int depth = 0;
        while (depth < 5 && !dfs(0, depth)) depth++;
        
        if (depth >= 5) cout << "5 or more" << endl;
        else cout << depth << endl;
    }

    return 0;
}

时间复杂度指数级，空间$O(n)$。