题目描述
给出一个长度为 m 的序列 A, 请你求出有多少种 1...n 的排列, 满足 A 是它的一个 LIS.
输入格式
第一行两个整数 n, m. 接下来一行 m 个整数, 表示 A.
输出格式
一行一个整数表示答案.
样例
Input:
5 3
1 3 4
Output:
11
数据范围与提示
对于前 30% 的数据, n ≤ 9;
对于前 60% 的数据, n ≤ 12;
对于 100% 的数据, 1 ≤ m ≤ n ≤ 15.
我们可以很自然的想到用f[S][s]表示目前选了S中的数,并且LIS的状态是s的方案数。这样的话转移就考虑把还没出现的数加入,再考虑其对s的影响就可以找到下一个状态了,但是这里还有两个需要注意的事情:
1.如何保证最后的LIS中有a[] ?
这个还是比较好解决的,我们只需要再两个地方判定就行了:1.如果加入的数是给出的a[]中的一个元素a[i],我们需要判断S中是否包含所有a[1],,,a[i-1];2.最后我们计算答案的时候,只能加s状态的LIS=m的f[2^n-1][s] ,因为这样才能保证a[]是LIS中的一个。
2.看起来状态是有 (2^n)^2 个的,会不会炸掉?
简单分析就可以知道,s肯定是S的一个子集,因为只有出现过了才可能在LIS的队列中出现。所以其实状态最多只有 3^n个。(但我真的写不出3进制的操作啊,只好用map强行套个log)
(然后我来回答一下为什么我考试的时候这个题爆零了23333,考试的时候我的程序就和现在的很接近了,但是我把题目看成LIS必须只能是给出的数列而不是有很多LIS共存,然后就怒WA10个点 (出题人来解释一下为什么这样能过样例233333))
#include<iostream> #include<cstdio> #include<cstdlib> #include<cstring> #include<cmath> #include<tr1/unordered_map> #include<algorithm> #define ll long long using namespace std; using namespace std::tr1; const int maxn=40005; unordered_map<int,int> mmp[maxn]; int ci[66],n,m,sum[maxn],a[20]; int to[maxn][20],pos[20],tmp[20]; inline void init(){ ci[0]=1; for(int i=1;i<=30;i++) ci[i]=ci[i-1]<<1; for(int i=1;i<ci[n];i++) sum[i]=sum[i^(i&-i)]+1; for(int i=1;i<ci[n];i++) for(int j=1;j<=n;j++){ to[i][j]=i; for(int k=j-1;k<n;k++) if(ci[k]&i){ to[i][j]^=ci[k]; break; } to[i][j]|=ci[j-1]; } } inline void solve(){ for(int i=1;i<=n;i++) if(pos[i]<=1) mmp[ci[i-1]][ci[i-1]]=1; for(int S=1;S<ci[n];S++) for(int s=S,TS,ts,now;s;s=(s-1)&S) if(mmp[S].count(s)){ now=mmp[S][s]; for(int k=1;k<=n;k++) if(!(ci[k-1]&S)){ if(pos[k]&&(S&tmp[pos[k]])!=tmp[pos[k]]) continue; TS=S|ci[k-1]; ts=to[s][k]; mmp[TS][ts]+=now; } } int ans=0; for(int i=1;i<ci[n];i++) if(sum[i]==m) ans+=mmp[ci[n]-1][i]; printf("%d ",ans); } int main(){ // freopen("arg.in","r",stdin); // freopen("arg.out","w",stdout); scanf("%d%d",&n,&m),init(); for(int i=1;i<=m;i++) scanf("%d",a+i),pos[a[i]]=i; tmp[1]=0; for(int i=2;i<=m;i++) tmp[i]=tmp[i-1]|ci[a[i-1]-1]; solve(); return 0; }