[2018雅礼3-25]arg dp套dp

考虑求LIS的算法，我们有一个数组didi表示上升子序列长度为ii时最后一项的最小值，然后每加入一个值xx，就替换掉dd中第一个≥x≥x的值。于是我们可以用一个三进制状态表示这个dp（00表示还未选，11表示选了且在dd中，22表示选了但不在dd中），转移即可。
复杂度O(3n∗nlogn)O(3n∗nlog⁡n)，实测把loewr_bound改成暴力扫的O(3n∗n2)O(3n∗n2)更快。。。
代码：

#include
#include
#include
#include
#define N 17
#define G(x,y) (x/mi[y-1]%3)
int n,m,a[N],R,f[14349000],id[14349000],mi[N],c[N],ans,nxt[N];
bool vis[N];
using namespace std;
int main()
{
    scanf("%d%d",&n,&m);
    for(int i=1;i<=m;i++)
    {
        scanf("%d",&a[i]);
        if(vis[a[i]]||(a[i]1,a[i-1]))||(a[i]>n)) {puts("0");return 0;}
        vis[a[i]]=1;
    }
    mi[0]=1;
    for(int i=1;i<=n;i++)
        mi[i]=mi[i-1]*3;    
    f[0]=1;
    for(int s=0;sif(f[s])
    {
        for(int i=1;i<=n;i++)
            c[i]=G(s,i);
        for(int i=1;i<=n;i++)
        {
            if(c[i]>0) continue;
            if(vis[i]&&(a[id[s]+1]!=i)) continue;
            int pos=-1,t=s;
            for(int j=i+1;j<=n;j++)
                if(c[j]==1) {pos=j;break;}
            t+=mi[i-1];if(pos!=-1) t+=mi[pos-1];
            f[t]+=f[s];
            id[t]=id[s]+(a[id[s]+1]==i); 
        }
    }
    for(int s=0;sif(f[s])
    {
        int cnt1=0,cnt2=0;
        for(int i=1;i<=n;i++)
        {
            if(G(s,i)>0) cnt2++;
            if(G(s,i)==1) cnt1++;
        }
        if(cnt1==m&&cnt2==n&&id[s]==m) ans+=f[s];
    }
    printf("%d",ans);
    return 0;
}