BZOJ3591: 最长上升子序列

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

code：

#include
#include
#include
#include
#include
#include
#include
#include
#include
#include
#include
#include
#include
#include
#include
#include
#include
#include
#include
#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<1<1<1<int f[mask];
queue<int>q;

int ans;

int main()
{
    pw[0]=1; for(int i=1;i1]*3;
    for(int i=0;i<15;i++) to[1<scanf("%d%d",&n,&m);
    for(int i=1;i<=m;i++)
    {
        scanf("%d",&a[i]);
        for(int j=1;j1<1;
    }
    int al=1<0]=-1; to3[0]=-1;
    for(int i=0;ifor(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<1;
            if(t[j]) ev|=1<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<1;
            if(t[j]) now|=1<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;
}