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;
}