【bzoj 3990】排序(DFS)

传送门biu~
显然一个操作序列是否合法与操作序列的顺序是无关的,所以我们只需要知道某种操作选不选即可。可以从小操作向大操作枚举,这样一个合法的操作序列对答案的贡献为选择的操作数的阶乘。
对于第i种操作我们将序列分成 2 n i ,每段长度 2 i 。找到序列中不递增的段,如果这样的段超过两个,显然无法实现;如果不存在这样的段,则不用操作;如果这样的段有一个,判断将这个段前后交换能否有序;如果有两个,则判断四种情况是否成立。

#include<bits/stdc++.h>
using namespace std;
int n,ans,jc[15],a[1<<13];
inline bool check(int x,int k){
    for(int i=1;i<(1<<k);++i)
        if(a[x+i]!=a[x+i-1]+1)  return false;
    return true;
}
inline void swap(int x,int y,int k){
    for(int i=0;i<(1<<k);++i)
        swap(a[x+i],a[y+i]);
}
void dfs(int k,int now){
    if(k==n+1){
        ans+=jc[now];
        return;
    }
    int c1=0,c2=0;
    for(int i=1;i<=(1<<n);i+=(1<<k))
        if(!check(i,k)){
            if(!c1) c1=i;
            else if(!c2) c2=i;
            else return;
        }
    if(!c1&&!c2)    dfs(k+1,now);
    else if(c1&&!c2){
        swap(c1,c1+(1<<k-1),k-1);
        if(check(c1,k)) dfs(k+1,now+1);
        swap(c1,c1+(1<<k-1),k-1);
    }
    else{
        for(int x=0;x<=1;++x){
            for(int y=0;y<=1;++y){
                swap(c1+x*(1<<k-1),c2+y*(1<<k-1),k-1);
                if(check(c1,k)&&check(c2,k))    dfs(k+1,now+1);
                swap(c1+x*(1<<k-1),c2+y*(1<<k-1),k-1);
            }
        }
    }
}
int main(){
    scanf("%d",&n);
    jc[0]=1;for(int i=1;i<=n;++i)    jc[i]=jc[i-1]*i;
    for(int i=1;i<=(1<<n);++i) scanf("%d",&a[i]);
    dfs(1,0);
    printf("%d\n",ans);
    return 0;
}

猜你喜欢

转载自blog.csdn.net/zp1ng/article/details/80170036