BZOJ3591: 最长上升子序列

考虑我们做最长上升子的过程，维护一个单调栈，每次加入一个元素，替换掉最前面一个比他大的数
我们可以dp这个加入元素的过程，用一个3进制的状态表示每个数在不在队列里且在不在栈里，可以用$O(n^22^n)$ 预处理对于栈中的每个状态，加入一个新的数之后的状态
然后做一个$O(n3^n)$ 的dp，注意要按当前序列内的元素数量从小到大dp

code：

#include<set>
#include<map>
#include<deque>
#include<queue>
#include<stack>
#include<cmath>
#include<ctime>
#include<bitset>
#include<string>
#include<vector>
#include<string>
#include<vector>
#include<cstdio>
#include<cstdlib>
#include<cstring>
#include<climits>
#include<complex>
#include<iostream>
#include<algorithm>
#define ll long long
#define lowbit(x) (x)&(-(x))
using namespace std;

const int maxn = 16;
const int mask = 14348907;
int pw[maxn];

int n,m;
int a[maxn],pre[maxn];
int nex[maxn],t[maxn],tp;
int to[1<<maxn],to3[1<<maxn],sum[1<<maxn],nxt[1<<maxn][maxn];
int f[mask];
queue<int>q;

int ans;

int main()
{
    pw[0]=1; for(int i=1;i<maxn;i++) pw[i]=pw[i-1]*3;
    for(int i=0;i<15;i++) to[1<<i]=i;

    scanf("%d%d",&n,&m);
    for(int i=1;i<=m;i++)
    {
        scanf("%d",&a[i]);
        for(int j=1;j<i;j++) pre[a[i]]|=1<<a[j]-1;
    }
    int al=1<<n; sum[0]=-1; to3[0]=-1;
    for(int i=0;i<al;i++)
    {
        for(int j=1;j<=n;j++) t[j]=i>>j-1&1;
        int ev=0;
        for(int j=1;j<=n;j++) 
        {
            if(!t[j])
                nxt[i][j]=(i^lowbit(i^ev))|1<<j-1;
            if(t[j]) ev|=1<<j-1;
        }

        sum[i]=sum[i^lowbit(i)]+1;
        to3[i]=to3[i^lowbit(i)]+pw[to[lowbit(i)]];
    }

    f[0]=1; q.push(0);
    while(!q.empty())
    {
        const int i=q.front(); q.pop();

        int &temp=f[i];
        int ii=0,now=0;
        for(int j=1;j<=n;j++)
        {
            t[j]=(i/pw[j-1])%3;
            if(t[j]==2) ii|=1<<j-1;
            if(t[j]) now|=1<<j-1;
        }
        bool ok=true;
        for(int j=1;j<=n;j++) if(!t[j])
        {
            ok=false;
            int nex=nxt[ii][j];
            if((pre[j]&now)==pre[j]&&sum[nex]<=m)
            {
                int ci=to3[now]+to3[nex]+pw[j-1];
                if(!f[ci]) q.push(ci);
                f[ci]+=temp;
            }
        }
        //printf("%d ",temp);
        //for(int i=1;i<=n;i++) printf("%d ",t[i]); putchar('\n');
        if(ok) ans+=temp;
    }
    printf("%d\n",ans);

    return 0;
}