bzoj3990 [SDOI2015]排序（搜索）

题目链接

一开始看见题目蒙了
方案数？那不多了去了吗
又读了一遍题，发现：每种操作最多可以执行一次
而N<=12，也就是说操作数最多就是12次，数据范围很小

首先我们很容易发现一个操作序列是否合法与序列的顺序是无关的
因此我们只要计算出序列需要进行多少种不同的操作（x种），最后答案就是x！

我手写了几个序列，发现都不能通过n次操作变成递增序列（这就有点尴尬了）：

6 7 8 5 4 2 1 3

这n个操作可以说是互不影响，因为每一种都是针对一种特定的区间，不同长度的区间是互不影响的
就拿这个序列来说吧，如果我们考虑第一种操作，我发现至少要交换成：7 8 5 6 3 4 1 2
才能考虑进行别的种类的操作，不然通过其他操作是怎么也不能使其有序的，但是要达到这种状态至少需要5步
因此这个序列无解

4 2 3 1

这个序列显然有解（只要通过一次第一种操作）

以上规律翻译成一种规范一点的语言：

如果我们进行第i种操作，我们要保证所有长度为2^i的区间中最多有两个区间不是递增且连续的，
如果有两个以上这样的区间，那么这个序列就是无解了

• 如果只有一个区间不递增，那么我们就交换ta的前后部分
• 如果有两个区间（x，y），（u，w）不递增，我们尝试4种交换方式：x—>u，x—>w，y—>u，y—>w
  （当中只有一种合法）

从小到大dfs每种操作，对于一个操作暴力判断是否合法/需要

具体剪枝和判断见上文

---

我一开始的想法是从小到大循环i，暴力判断是否需要这种操作
如果存在的话，计数器cnt++，最后的答案cnt！
但是呢，但是WA了。。。

Q.

难道一个序列需要什么操作是不确定的吗？
换句话说，某个操作能用若干个其他操作代替吗？
或者说，会因为进行不同的合法操作而导致原来有解的序列无解吗？

A.

经过我疯狂的对拍后，构造出了一组数据：

2
3 2 1 4

ans=3
但是我YY的方法跑出了：1

方案：
1. change(1,3)
2. change((3,2),(1,4)) change(2,4)
3. change(2,4) change((3,2),(1,4))

显然我YY的做法计算的是最小步数，然而合法的方案可以不用最小步数

tip

回溯的时候不要忘了要change回来
这道题一句话题解就是手玩猜结论
省选的搜索题都不好做。。。

#include<cstdio>
#include<cstring>
#include<iostream>
#define ll long long

using namespace std;

const int N=5000;
ll jc[15],ans=0;
int n,a[N];

int xu(int l,int r)    //判断区间是否有序 
{
    for (int i=l;i<r;i++)
        if (a[i]!=a[i+1]-1) return 0;
    return 1;
}

void change(int l,int r,int len)
{
    for (int i=1;i<=len;i++)
        swap(a[i+l-1],a[i+r-1]);
}

void dfs(int t,int cnt)
{
    if (t>n)
    {
        ans+=jc[cnt];
        return;
    }
    int len=1<<t;
    int wh1=-1,wh2=-1;
    for (int i=1;i<=(1<<n);i+=len)
        if (!xu(i,i+len-1)) 
            if (wh1==-1) wh1=i;
            else if (wh2==-1) wh2=i;
            else return;             //不合法    

    if (wh1==-1) dfs(t+1,cnt);
    else if (wh2==-1) {              //只有一个不递增序列 
        change(wh1,wh1+len/2,len/2);
        dfs(t+1,cnt+1);
        change(wh1,wh1+len/2,len/2);
    }
    else {                           //尝试四种交换方法 
        change(wh1,wh2,len/2); 
        if (xu(wh1,wh1+len-1)&&xu(wh2,wh2+len-1)) dfs(t+1,cnt+1); 
        change(wh1,wh2,len/2);

        change(wh1+len/2,wh2,len/2); 
        if (xu(wh1,wh1+len-1)&&xu(wh2,wh2+len-1)) dfs(t+1,cnt+1); 
        change(wh1+len/2,wh2,len/2);

        change(wh1,wh2+len/2,len/2); 
        if (xu(wh1,wh1+len-1)&&xu(wh2,wh2+len-1)) dfs(t+1,cnt+1); 
        change(wh1,wh2+len/2,len/2);

        change(wh1+len/2,wh2+len/2,len/2); 
        if (xu(wh1,wh1+len-1)&&xu(wh2,wh2+len-1)) dfs(t+1,cnt+1); 
        change(wh1+len/2,wh2+len/2,len/2);
    }
}

int main()
{
    jc[1]=jc[0]=1;
    for (int i=2;i<=12;i++) jc[i]=jc[i-1]*(ll)i;

    scanf("%d",&n);
    for (int i=1;i<=(1<<n);i++)
        scanf("%d",&a[i]);

    dfs(1,0);
    printf("%lld",ans);
    return 0;
}