【BZOJ3990】【SDOI2015】排序

Description

 小A有一个1-2^N的排列A[1..2^N],他希望将A数组从小到大排序,小A可以执行的操作有N种,每种操作最多可以执行一次,对于所有的i(1<=i<=N),第i中操作为将序列从左到右划分为2^{N-i+1}段,每段恰好包括2^{i-1}个数,然后整体交换其中两段.小A想知道可以将数组A从小到大排序的不同的操作序列有多少个,小A认为两个操作序列不同,当且仅当操作个数不同,或者至少一个操作不同(种类不同或者操作位置不同).

下面是一个操作事例:

N=3,A[1..8]=[3,6,1,2,7,8,5,4].

第一次操作,执行第3种操作,交换A[1..4]和A[5..8],交换后的A[1..8]为[7,8,5,4,3,6,1,2].

第二次操作,执行第1种操作,交换A[3]和A[5],交换后的A[1..8]为[7,8,3,4,5,6,1,2].

第三次操作,执行第2中操作,交换A[1..2]和A[7..8],交换后的A[1..8]为[1,2,3,4,5,6,7,8].

Input

第一行,一个整数N

第二行,2^N个整数,A[1..2^N]

Output

一个整数表示答案

Sample Input

3
7 8 5 6 1 2 4 3

Sample Output

6

HINT

100%的数据, 1<=N<=12.

题解

非常好的一道搜索题。

首先可以发现如果最后选了j个操作，那么这j个操作的顺序是不影响最后结果的，所以直接加上个数的阶乘就行了。所以就判断这个操作是不是选就行了。我们找到序列中不是连续递增的段，如果这样的段超过2个就不行，因为以后分的段都比这次的段大。所以如果没有这样的段，就不需要执行这个操作如果有一个这样的段，判断将这个段的前半部分和后半部分交换后是否连续递增，如果是就交换然后继续DFS。最后把答案累加就行了。

#include<iostream>
#include<cstdio>
#include<cstring>
#include<algorithm>
#include<cmath>
#define ll long long
using namespace std;
const int N=10010;
ll ans,f[N];
int n,a[N];
inline bool pan(int x,int y){
	for(int i=x;i<y;i++)
	if(a[i]!=a[i+1]-1) return 0;
	return 1; 
}
inline void swap(int x,int y,int l){
	for(int i=0;i<l;i++) swap(a[x+i],a[y+i]);
}
void dfs(int i,int cnt){
	if(i==n){
		ans+=f[cnt];
		return ;
    }
    int b[10],tot=0;
    for(int j=1;j<=(1<<n);j+=(1<<(i+1))){
		if(!pan(j,j+(1<<(i+1))-1)){
			if(tot==4) return ;
			b[++tot]=j;b[++tot]=j+(1<<i);
		}
    }
    if(!tot) dfs(i+1,cnt);
    if(tot==2){
		if(a[b[2]]+(1<<i)==a[b[1]]){
			swap(b[1],b[2],1<<i);
			dfs(i+1,cnt+1);
			swap(b[1],b[2],1<<i);
		}
    }
    if(tot==4){
		if(a[b[3]]+(1<<i)==a[b[2]]&&a[b[1]]+(1<<i)==a[b[4]]){
			swap(b[1],b[3],1<<i);
			dfs(i+1,cnt+1);
			swap(b[1],b[3],1<<i);
		}
		if(a[b[1]]+(1<<i)==a[b[3]]&&a[b[2]]+(1<<i)==a[b[4]]){
			swap(b[2],b[3],1<<i);
			dfs(i+1,cnt+1);
			swap(b[2],b[3],1<<i);
		}
		if(a[b[4]]+(1<<i)==a[b[2]]&&a[b[3]]+(1<<i)==a[b[1]]){
			swap(b[1],b[4],1<<i);
			dfs(i+1,cnt+1);
			swap(b[1],b[4],1<<i);
		}
		if(a[b[1]]+(1<<i)==a[b[4]]&&a[b[3]]+(1<<i)==a[b[2]]){
			swap(b[2],b[4],1<<i);
			dfs(i+1,cnt+1);
			swap(b[2],b[4],1<<i);
		}
    }
}
int main(){
   scanf("%d",&n);
   for(int i=1;i<=(1<<n);i++) scanf("%d",&a[i]);
   f[0]=1;
   for(int i=1;i<=n;i++) f[i]=f[i-1]*i;
   
   
   dfs(0,0);
   printf("%d\n",ans);
   return 0;
}