搜索——BZOJ3990/Luogu3322 [SDOI2015]排序

http://www.lydsy.com/JudgeOnline/problem.php?id=3990

https://www.luogu.org/problem/show?pid=3322

部分文字转自lc233，博主添加文字会另外注明

http://blog.csdn.net/largecub233/article/details/67633368

假设分块的长度是2^i

那我们把i从1到n枚举；

枚举到i时，会分成(2^n)/(2^i)块；

我们就让这几个块内的值递增，且每个相邻的数相差1；

如果我们不能维护成功，那么return；

怎么维护呢；

如果有n个块时不递增的；

那么n>2显然就无法维护了；

n=1我们可以试试把这个块分成两半，然后交换两半，看看可不可以；

（博主注：这里的可不可以表示连续递增）

n=2那我们把两个大块分成4个小块，两两组合，组合时要保证只能交换块一次；

（博主注：这里的两两交换判断意思是：

1.第一大块的前半部分（小块）和第二大块的前半部分（小块）交换后是否连续递增

2.第一大块的后半部分（小块）和第二大块的后半部分（小块）交换后是否连续递增

3.第一大块的前半部分（小块）和第二大块的后半部分（小块）交换后是否连续递增

4.第一大块的后半部分（小块）和第二大块的前半部分（小块）交换后是否连续递增）

博主注：最后贴一下我的代码，表示又被C++下标坑了。。。

#include<bits/stdc++.h>
using namespace std;
int n,a[10001];
long long ans=0,jie[101];
inline void jh(int x,int y,int l){for(int i=0;i<l;i++)swap(a[x+i],a[y+i]);}
inline void dfs(int x,int sum){
if(x==n){
ans+=jie[sum];
return;
}
int s[3]={0};
int jzq=1<<x+1;
for(int i=1,j=1<<n;i<=j;i+=jzq)if(a[i+(jzq>>1)-1]+1!=a[i+(jzq>>1)]){
if(s[0]==2)return;
s[++s[0]]=i;
}
if(!s[0]){dfs(x+1,sum);return;}
if(s[0]==1){
int i=s[1];
if(a[i+jzq-1]+1==a[i]){
jh(i,i+(jzq>>1),jzq>>1);
dfs(x+1,sum+1);
jh(i,i+(jzq>>1),jzq>>1);
}
return;
}
int i=s[1],j=s[2];
if(a[i+(jzq>>1)-1]+1==a[j+(jzq>>1)]&&a[j+(jzq>>1)-1]+1==a[i+(jzq>>1)]){
jh(i,j,jzq>>1);
dfs(x+1,sum+1);
jh(i,j,jzq>>1);
jh(i+(jzq>>1),j+(jzq>>1),jzq>>1);
dfs(x+1,sum+1);
jh(i+(jzq>>1),j+(jzq>>1),jzq>>1);
}
if(a[j+(jzq>>1)-1]+1==a[i]&&a[j+jzq-1]+1==a[i+(jzq>>1)]){
jh(i,j+(jzq>>1),jzq>>1);
dfs(x+1,sum+1);
jh(i,j+(jzq>>1),jzq>>1);
}
if(a[i+(jzq>>1)-1]+1==a[j]&&a[i+jzq-1]+1==a[j+(jzq>>1)]){
jh(j,i+(jzq>>1),jzq>>1);
dfs(x+1,sum+1);
jh(j,i+(jzq>>1),jzq>>1);
}
}
int main()
{
scanf("%d",&n);
for(int i=1,j=1<<n;i<=j;i++)scanf("%d",&a[i]);
jie[1]=1;for(int i=2;i<=n;i++)jie[i]=jie[i-1]*i;
dfs(0,0);
printf("%lld",ans);
return 0;
}