jzoj5600 Arg

Description

给出一个长度为 m 的序列 A, 请你求出有多少种 1…n 的排列, 满足 A 是它的一个 LIS.

Input

第一行两个整数 n,m.
接下来一行 m 个整数, 表示 A.

Output

一行一个整数表示答案.

Sample Input

5 3
1 3 4

Sample Output

11

Data Constraint

对于前 30% 的数据, n ≤ 9;
对于前 60% 的数据, n ≤ 12;
对于 100% 的数据, 1 ≤ m ≤ n ≤ 15.

Solution

数据范围很状压
我们知道很久很久以前有一个大牛发明了特别的LCS算法
维护一个f数组，f[i]表示长度为i的LCS的末尾的最小值
那么我们每次更新一下就可以了
现在我们考虑状压
定义一个三进制状态
0表示某个数x不在当前集合中，1表示x在集合中并且在f数组中，2表示在集合中但不在f数组中
那么我们可以由此转移
由于状态比较多枚举的时候可以加一点小优化
具体可以看代码里面的注释

Code

#include <cstdio>
#include <cstring>
#include <algorithm>
#define fo(i,a,b) for(i=a;i<=b;i++)
using namespace std;
typedef long long ll;
const ll inf=14348910;
ll n,m,i,j,fg,cnt,num,ans,x,bit[40],a[40],vis[40],f[inf],k;
int main(){
    freopen("1.in","r",stdin);
    freopen("1.out","w",stdout);
    scanf("%lld%lld",&n,&m); bit[0]=1;
    fo(i,1,m) scanf("%lld",&a[i]),vis[a[i]]=i;
    fo(i,1,n) bit[i]=bit[i-1]*3;
    f[0]=1;
    fo(i,0,bit[n]-1) 
      if (f[i])
        fo(j,1,n) {
          if ((i/bit[j-1])%3>0) continue;
          if (vis[j]&&vis[j]!=1) 
            if ((i/bit[a[vis[j]-1]-1])%3==0) continue; //剪枝，如果LCS中前面的没填就填后面的                     显然不行
          fg=0;
          fo(k,j+1,n)
              if ((i/bit[k-1])%3==1) {
                f[i+bit[k-1]+bit[j-1]]+=f[i],fg=1;
                //如果k在f中,但是j比k小，那么j可以替换k，可以转移
                break;}
         if (!fg) 
           f[i+bit[j-1]]+=f[i];
           //如果没有可以替换的，直接接在后面
      }
    fo(i,bit[n-1],bit[n]-1) {
        x=i,cnt=0,num=0;
        while (x) {
            if (x%3==0) break; //如果还有未填的可以直接break，这样剪枝的效果很好
            if (x%3==1) cnt++; else if (x%3==2)num++;
            x/=3;
        }
        if (cnt==m&&cnt+num==n) ans+=f[i];
    }
    printf("%lld\n",ans);
}