20200721 T1 s1mple【容斥，划分】

题目描述

$n\le 17,q\le 10^5$

题目分析

先不考虑0的限制，相当于强制一些链连通，其它任意，这样很多排列状态就变成了本质相同的划分方案。

Code：

#include<bits/stdc++.h>
#define maxn 17
#define LL long long
using namespace std;
int n,N,Q,lg[1<<maxn],bit[1<<maxn];
LL f[1<<maxn][maxn],F[maxn+1][1<<maxn],P[maxn+1][maxn+1][1<<maxn],ans[1<<maxn],A[1<<maxn];
bool B[maxn][maxn];
char str[maxn+5];
void Mul(LL *a,LL *b,LL *c){for(int i=0;i<=N;i++) a[i]=b[i]*c[i];}
void FWT_OR(LL *a,int len){
	for(int i=2,l=1;i<=len;l=i,i<<=1) for(int j=0;j<len;j+=i) for(int k=j;k<j+l;k++) a[k+l]+=a[k];
}
void IFWT_AND(LL *a,int len){
	for(int i=2,l=1;i<=len;l=i,i<<=1) for(int j=0;j<len;j+=i) for(int k=j;k<j+l;k++) a[k]-=a[k+l];
}
int cnt[maxn+1];
int encode(int s){
	memset(cnt,0,sizeof cnt); int l=0;
	for(int i=0;i<n-1;i++) if(l++,!(s>>i&1)) cnt[l]++,l=0;
	cnt[l+1]++,s=0;
	for(int i=1;i<=n;i++) for(int j=1;j<=cnt[i];j++) s=s<<i|((1<<i-1)-1);
	return s;
}
int main()
{
	freopen("s1mple.in","r",stdin);
	freopen("s1mple.out","w",stdout);
	scanf("%d",&n),N=(1<<n)-1;
	for(int i=0;i<n;i++) {scanf("%s",str); for(int j=0;j<n;j++) B[i][j]=str[j]=='1';}
	for(int i=(bit[1]=1,2);i<1<<n;i++) lg[i]=lg[i>>1]+1,bit[i]=bit[i>>1]+(i&1);
	for(int i=0;i<n;i++) f[1<<i][i]=1;
	for(int s=1;s<N;s++)
		for(int t=s,i;t;t-=1<<i) if(f[s][i=lg[t]])
			for(int r=N^s,j;r;r-=1<<j) if(B[i][j=lg[r]])
				f[s|1<<j][j]+=f[s][i];
	for(int s=1;s<=N;s++){
		LL x=0,c=0; for(int i=0;i<n;i++) if(s>>i&1) c++,x+=f[s][i];
		F[c][s]=x;
	}
	for(int i=1;i<=n;i++) FWT_OR(F[i],N+1);
	for(int i=1;i<=n;i++){
		memcpy(P[i][1],F[i],sizeof F[i]);
		for(int j=2;i*j<=n;j++) Mul(P[i][j],P[i][j-1],F[i]);
	}
	memset(ans,-1,sizeof ans);
	for(int s=0;s<=N>>1;s++){
		int t=encode(s);
		if(~ans[t]) {ans[s]=ans[t];continue;}
		for(int i=0;i<=N;i++) A[i]=1;
		for(int i=1;i<=n;i++) if(cnt[i]) Mul(A,A,P[i][cnt[i]]);
		LL ret=0;
		for(int t=N;;t=(t-1)&N) {ret+=(n-bit[t]&1?-1:1)*A[t]; if(!t) break;}
		ans[t]=ans[s]=ret;
	}
	IFWT_AND(ans,(N>>1)+1);
	scanf("%d",&Q);
	while(Q--){
		scanf("%s",str); int s=0;
		for(int i=0;i<n-1;i++) s|=(str[i]-'0')<<i;
		printf("%lld\n",ans[s]);
	}
	cerr<<clock()<<endl;
}