[SDOI2015]排序

题目

传送门 to usOJ

传送门 to luogu（不想搜了）

思路

大佬的博客讲得很有道理，但是我这种蒟蒻看不懂，自闭了

因为第 $x$ 种操作可以完全把第 $x-1$ 种操作给“套住”，所以，对于合法的操作方案，交换两种操作的先后，仍然存在一个合法操作序列。比如，你先将 $[x,x+2^k)$ 跟别人换了，后来又将 $[x,x+2^r)$ 挪走了（满足 $r>k$ ）；你也可以先将 $[x,x+2^r)$ 挪走，然后在 原来的$[x,x+2^r)$中交换前$2^k$个元素。

再用一个形象的譬喻。你可以先到厕所（大区间移动），再和别人换一条裤子（小区间交换）；也可以先把裤子交换，然后再走到厕所，就是可能被朋友发现。

起到关键作用的是，大区间移动包含了小区间交换——你到厕所的时候，你的裤子其实是跟着你一起到了厕所的。并没有拟人呢

有了这个思考之后，我们不妨从第 $1$ 种操作考虑到第 $n$ 种。在第 $i$ 种操作进行完之后，假如还有一个区间 $[x,x+2^i)$ 非 连续递增正整数，并且这个区间是成功地被划分出来的区间（即，$2^i|x$ ），那么靠以后的操作，能力挽狂澜吗？

显然是不能的，因为你再也无法解决 $[x,x+2^i)$ 的内部顺序了，以后的区间交换只会越来越长，并且没有分割线会落在这个区间内。

所以，我们考虑第 $i$ 种操作的时候，一定要把每一个划分出来的小区间搞成连续递增。最后，乘一个 $m!$ ，其中 $m$ 是操作的个数——诚如我前面所说，交换顺序仍然存在一个合法方案。

代码就很简单了，代码还得好好想想咋写。交换的时候得枚举一下呢。

代码

#include <cstdio>
#include <iostream>
#include <algorithm>
using namespace std;
typedef long long int_;
inline int readint(){
	int x; scanf("%d",&x); return x;
}
inline void writeint(int x){
	printf("%d",x);
}

# define FOR(i,n) for(int i=0; i<(n); ++i)

const int MaxN = 12;

int a[1<<MaxN], n;
long long ans = 0;
long long jc[MaxN+1] = {1};

void change(int x,int y,int l){
	FOR(i,l) swap(a[x+i],a[y+i]);
} // 本质是swap 
void dfs(int l,int x){
// 操作的长度为 l ，已经操作了 x 次
	if(l == 1<<n) // 考虑了所有操作
		return void(ans += jc[x]);
	int L = l<<1, f = 0, p[4];
	for(int i=0; i<(1<<n); i+=L)
    	if(a[i]%L or a[i+l] != l+a[i]){
			p[++ f] = i; // 记录一下
			if(f > 2) return ;
		} // 超过两个当场去世
	if(f == 0) dfs(L,x); // 不操作
	else if(f == 1){
		change(p[1],p[1]+l,l);
		dfs(L,x+1); // 交换一次
		change(p[1],p[1]+l,l);
	}
	else FOR(i,2) FOR(j,2){ // 枚举如何交换
		change(p[1]+i*l,p[2]+j*l,l);
		if(a[p[1]]%L == 0) // 四个并列条件 
		if(a[p[2]]%L == 0)
		if(a[p[1]]+l == a[p[1]+l])
		if(a[p[2]]+l == a[p[2]+l])
			dfs(L,x+1);
		change(p[1]+i*l,p[2]+j*l,l);
	}
}

int main(){
	n = readint();
	FOR(i,n) jc[i+1] = jc[i]*(i+1);
	FOR(i,1<<n) a[i] = readint()-1;
	dfs(1,0); printf("%lld\n",ans);
	return 0;
}