bzoj 4000: [TJOI2015]棋盘

矩乘，题意坑爹，那个第一行是中间那一行，所以记录一行的状态就可以了。f[S]表示最后一行状态为S的方案数，然后暴力看能否转移。

   
   

    
    

     
     #include<cstdio>
    
    
    
    

     
     #include<cstdlib>
    
    
    
    

     
     #include<cstring>
    
    
    
    

     
     #include<algorithm>
    
    
    
    

     
      
    
    
    
    

     
     #define ll long long
    
    
    
    

     
     #define ui unsigned int 
    
    
    
    

     
     #define inf 1e9
    
    
    
    

     
     #define eps 1e-8
    
    
    
    

     
     #define md
    
    
    
    

     
     #define N 70
    
    
    
    

     
     using namespace std;
    
    
    
    

     
     int m;
    
    
    
    

     
     bool c[20],d[20];
    
    
    
    

     
     int list[3][20];
    
    
    
    

     
     bool safe[N];
    
    
    
    

     
     struct Ju
    
    
    
    

     
     {
    
    
    
    

     
     	ui a[N][N];
    
    
    
    

     
     	int x,y;
    
    
    
    

     
     } a,ju,mid,ans;
    
    
    
    

     
     bool operator *= (Ju &a,Ju &b)
    
    
    
    

     
     {
    
    
    
    

     
     	memset(mid.a,0,sizeof(mid.a));
    
    
    
    

     
     	mid.x=a.x; mid.y=b.y;
    
    
    
    

     
     	for (int i=0;i<a.x;i++)
    
    
    
    

     
     	  for (int k=0;k<a.y;k++)
    
    
    
    

     
     	    if (a.a[i][k])
    
    
    
    

     
     	      for (int j=0;j<b.y;j++)
    
    
    
    

     
     	        mid.a[i][j]+=a.a[i][k]*b.a[k][j];
    
    
    
    

     
     	a=mid;
    
    
    
    

     
     }
    
    
    
    

     
      
    
    
    
    

     
     void kpow(Ju &a,int b)
    
    
    
    

     
     {
    
    
    
    

     
     	memset(ans.a,0,sizeof(ans.a));
    
    
    
    

     
     	ans.x=a.x; ans.y=a.y;
    
    
    
    

     
     	for (int i=0;i<a.x;i++) ans.a[i][i]=1;
    
    
    
    

     
     	while (b)
    
    
    
    

     
     	{
    
    
    
    

     
     		if (b&1) ans*=a;
    
    
    
    

     
     		a*=a; b>>=1;
    
    
    
    

     
     	}
    
    
    
    

     
     	a=ans;
    
    
    
    

     
     }
    
    
    
    

     
     	    
    
    
    
    

     
      
    
    
    
    

     
     bool ok(bool a[],bool b[],int c[])
    
    
    
    

     
     {
    
    
    
    

     
     	for (int i=0;i<m;i++)
    
    
    
    

     
     	{
    
    
    
    

     
     		if (!a[i]) continue;
    
    
    
    

     
     		for (int j=1;j<=c[0];j++)
    
    
    
    

     
     		{
    
    
    
    

     
     			int t=i+c[j];
    
    
    
    

     
     			if (t>=0&&t<m&&b[t]) return 0;
    
    
    
    

     
     		}
    
    
    
    

     
     	}
    
    
    
    

     
     	return 1;
    
    
    
    

     
     }
    
    
    
    

     
      
    
    
    
    

     
     bool get_safe(int S)
    
    
    
    

     
     {
    
    
    
    

     
     	for (int i=0;i<m;i++) c[i]=(S&(1<<i))>0;
    
    
    
    

     
     	return ok(c,c,list[1]);
    
    
    
    

     
     }
    
    
    
    

     
      
    
    
    
    

     
     bool check(int s1,int s2)
    
    
    
    

     
     {
    
    
    
    

     
     	if (!safe[s1]||!safe[s2]) return 0;
    
    
    
    

     
     	for (int i=0;i<m;i++) c[i]=(s1&(1<<i))>0,d[i]=(s2&(1<<i))>0;
    
    
    
    

     
     	return ok(c,d,list[2])&&ok(d,c,list[0]);
    
    
    
    

     
     }
    
    
    
    

     
      
    
    
    
    

     
     int main()
    
    
    
    

     
     {
    
    
    
    

     
     	int n,P,K;
    
    
    
    

     
     	scanf("%d%d",&n,&m);
    
    
    
    

     
     	scanf("%d%d",&P,&K);
    
    
    
    

     
     	for (int i=0;i<3;i++)
    
    
    
    

     
     	  for (int j=0;j<P;j++)
    
    
    
    

     
     	  {
    
    
    
    

     
     	  	int x;
    
    
    
    

     
     	  	scanf("%d",&x);
    
    
    
    

     
     	  	if (x&&!(i==1&&j==K)) list[i][++list[i][0]]=j-K;
    
    
    
    

     
     	  }
    
    
    
    

     
     	int S=(1<<m);
    
    
    
    

     
     	for (int i=0;i<S;i++) safe[i]=get_safe(i);
    
    
    
    

     
     	ju.x=ju.y=S;
    
    
    
    

     
     	for (int i=0;i<S;i++)
    
    
    
    

     
     	  for (int j=0;j<S;j++)
    
    
    
    

     
     	    ju.a[i][j]=check(i,j);
    
    
    
    

     
     	kpow(ju,n);
    
    
    
    

     
     	a.x=1; a.y=S; a.a[0][0]=1; a*=ju;
    
    
    
    

     
     	ui ans=0;
    
    
    
    

     
     	for (int i=0;i<S;i++) ans+=a.a[0][i];
    
    
    
    

     
     	printf("%u\n",ans);
    
    
    
    

     
     	return 0;
    
    
    
    

     
     }