BZOJ 3591 最长上升子序列

题目链接

题意：
给定 $k$ 个数，求所有长度为 $n$ 且最长上升子序列的排列为给定的 $k$ 个数的方案数。
$n<=15$

首先想到可以用状压解决排列的顺序问题，而仅通过 $0,1$ 无法表达出最长上升子序列的长度，而再开一维进行表示也无法进行有效转移。那么考虑在求是 $LIS$ 时我们会对一个数组不断进行覆盖，每个数字会逐个填入或被覆盖，那么如果我们用 $0,1,2$ 三个状态分别表示未用、已用、已被覆盖，即可通过此得到最长上升子序列的长度及其他序列信息，并可以轻松转移。而考虑到复杂度 $O(n{3}^n)$ 已经很大，在对每个三进制状态进行枚举时，按从小到大的顺序进行枚举，另一个计算最长上升子序列长度的指针满足单调性，因此不会增加额外复杂度。

#include<bits/stdc++.h>
using namespace std;

const int inf=0x3f3f3f3f;
int dp[15000000];
int is[20],pre[20],ppow[20],a[20],lis[20];

int main(){
    ios::sync_with_stdio(false);
    cin.tie(0);

    int n,K;
    cin>>n>>K;
    for(int i=0;i<K;i++){
        cin>>a[i];
        if(i>0) pre[a[i]]=a[i-1];
    }
    ppow[1]=is[0]=1;
    for(int i=2;i<=n+1;i++) ppow[i]=ppow[i-1]*3;

    int res=0;
    dp[0]=1;
    for(int i=0;i<(ppow[n+1]);i++){
        if(!dp[i]) continue;
        int cnt=1,tot=0;
        for(int j=1,temp=i;j<=n;j++,temp/=3){
            is[j]=0;
            if(temp%3) is[j]=1,tot++;
            if(temp%3==1) lis[cnt++]=j;
        }
        if(tot==n){
            res+=dp[i];
            continue;
        }
        lis[cnt]=inf;
        int k=1;
        for(int j=1;j<=n;j++){
            if(is[j] || !is[pre[j]]) continue;
            while(lis[k]<j) k++;
            if(k>K) continue;
            if(k==cnt) dp[i+ppow[j]]+=dp[i];
            else dp[i+ppow[j]+ppow[lis[k]]]+=dp[i];
        }
    }
    cout<<res<<endl;
}