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