xsy 排序

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

一个整数表示答案

---

Solution

注意 交换位置是相邻的几个数(eg,相邻的12)而不是真正的位置(eg.a[1]与a[5])

首先我们很容易发现一个操作序列是否合法与序列的顺序是无关的

因此我们只需要确定某个操作序列中每个操作选不选就行了 那么这类操作序列对答案的贡献就是选择的操作数的阶乘

我们从小到大DFS，对于第i次操作我们将序列分成$2^{(n-i)}$段，每段长度$2^i$

我们找到序列中不是连续递增的段，如果这样的段超过2个，显然就废了

如果没有这样的段，就不需要执行这个操作

如果有一个这样的段，判断将这个段的前半部分和后半部分交换后是否连续递增，如果是就交换然后继续DFS

如果有两个这样的段，判断四种交换情况然后DFS

如果有三个这样的段，没法解决，直接返回

Code

#include<bits/stdc++.h>
using namespace std;
typedef long long ll;
int a[10010],n,m,tot;
ll ans,fac[13];
bool check(int pos,int len){
    for(int i=1;i<len;i++) 
    if(a[pos+i]!=a[pos+i-1]+1){
        //cout<<pos<<" "<<i<<" "<<a[pos+i]<<" "<<a[pos+i-1]<<" "<<len<<endl; 
        return false;
    }
    return true;
}
void swp(int u,int v,int len){
    for(int i=0;i<len;i++) 
    swap(a[u+i],a[v+i]);
}
void dfs(int dep){
    //cout<<dep<<endl;
    if(dep>n){
        if(check(1,m)) 
        ans+=fac[tot];
        return;
    }
    int len=1<<dep,x=0,y=0;
    for(int i=1;i<=m;i+=len)
    if(check(i,len)==0){
        if(!x) x=i;
        else if(!y) y=i;
        else{
            //cout<<"???"<<" "<<x<<" "<<y<<" "<<i<<" "<<len<<" "<<check(x,2)<<endl;
            return;
        }
    }
    if(x==0&&y==0) dfs(dep+1);
    if(x&&y==0){
        tot++;
        swp(x,x+(len>>1),(len>>1));
        dfs(dep+1);
        swp(x,x+(len>>1),(len>>1));
        tot--;
    }
    if(x&&y){
        tot++;
        for(int i=0;i<=1;i++)
        for(int j=0;j<=1;j++){
            swp(x+i*(len>>1),y+j*(len>>1),(len>>1));
            if(check(x,len)&&check(y,len)) dfs(dep+1);
            swp(x+i*(len>>1),y+j*(len>>1),(len>>1));
        }
        tot--;
    }
}
int main(){
    fac[0]=1;
    for(int i=1;i<=12;i++) 
    fac[i]=fac[i-1]*i;
    scanf("%d",&n);
    m=(1<<n);
    for(int i=1;i<=m;i++) 
    scanf("%d",&a[i]);//,cout<<i<<" "<<a[i]<<endl;;
    dfs(1);
    printf("%lld",ans);
}