ARC093F - Dark Horse

题目链接
题目相当于有长度为$2^0,2^1,...,2^{n-1}$这$n$个区间，要在上面填上$2,...,n$，使得每个区间的最小值都不为给定的$m$个数的方案数（最后再乘个$n$）。
显然最小值不为某些数是不好计算的，考虑容斥。又注意到$m$很小，枚举哪些人一定是区间最小的，按照编号从大到小考虑那些人。记$f[i][j]$为考虑了前$i$个人，$n$个区间分别是否已经有一定是最小值的人，只考虑这些区间的排列情况的方案数，转移时枚举这个人是否一定最小，以及如果是，这个人在哪个区间，乘上组合数即可。
代码：

#include<bits/stdc++.h>
using namespace std;
typedef long long ll;
const int N=17,M=(1<<N);
const ll P=1e9+7;
ll fc[M],iv[M];
ll A(int n,int m){
	if(m<0 || m>n) return 0;
	return fc[n]*iv[n-m]%P;
}
ll fpw(ll a,ll x){
	ll s=1;
	for(;x;x>>=1,(a*=a)%=P) if(x&1) (s*=a)%=P;
	return s;
}
int n,m,a[N],p[N];
ll f[N][M],ans;
int main()
{
	//freopen(".in","r",stdin);
	//freopen(".out","w",stdout);
	fc[0]=1; for(int i=1;i<M;i++) fc[i]=fc[i-1]*i%P;
	iv[M-1]=fpw(fc[M-1],P-2);
	for(int i=M-2;~i;i--) iv[i]=iv[i+1]*(i+1)%P;
	cin>>n>>m; p[0]=1; for(int i=1;i<=n;i++) p[i]=(p[i-1]<<1);
	for(int i=1;i<=m;i++) scanf("%d",&a[i]);
	reverse(a+1,a+m+1); f[0][0]=1;
	for(int i=0;i<m;i++){
		for(int j=0;j<p[n];j++) if(f[i][j]){
			//printf("i=%d j=%d f=%lld\n",i,j,f[i][j]);
			ll x=f[i][j]; int s=0;
			for(int k=0;k<n;k++) if(j&p[k]) s+=p[k];
			(f[i+1][j]+=x)%=P;
			for(int k=0;k<n;k++) if(!(j&p[k])){
				//cout<<k<<" "<<x*A(p[n]-a[i+1]-s,p[k]-1)%P*p[k]<<endl;
				(f[i+1][j+p[k]]+=x*A(p[n]-a[i+1]-s,p[k]-1)%P*p[k])%=P;
			}
		}
	}
	for(int i=0;i<p[n];i++){
		int s=0,t=1;
		for(int j=0;j<n;j++){
			if(!(i&p[j])) s+=p[j];
			else t=-t;
		}
		(ans+=P+fc[s]*f[m][i]%P*t)%=P;
	}
	cout<<ans*p[n]%P<<endl;
	return 0;
}