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

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

#include<iostream>
#include<cstdio>
#include<cstring>
#include<algorithm>
#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]<max(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;s<mi[n];s++)
    if(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;s<mi[n];s++)
    if(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;
}
评论
添加红包

请填写红包祝福语或标题

红包个数最小为10个

红包金额最低5元

当前余额3.43前往充值 >
需支付:10.00
成就一亿技术人!
领取后你会自动成为博主和红包主的粉丝 规则
hope_wisdom
发出的红包
实付
使用余额支付
点击重新获取
扫码支付
钱包余额 0

抵扣说明:

1.余额是钱包充值的虚拟货币,按照1:1的比例进行支付金额的抵扣。
2.余额无法直接购买下载,可以购买VIP、付费专栏及课程。

余额充值