这题和之前的一个 状态压缩二,是很类似的,这道题的n比较大,开不下数组,无法使用之前的类似方法,但是由于k比较小,可以使用状态压缩,结合这道题目的前2个版本,考虑使用矩阵快速幂进行计算,关键是计算出状态转移的矩阵M,使用dfs,跑出状态转移矩阵
用二进制表示状态,
上一行的状态用y表示
当前这一行的状态用x表示
枚举当前这一行的每一个空的位置,放置骨牌,有3种情况,题目的提示讲解的很清楚
边界就是当前的这一行最后一个位置也放上了数,状态x就转移为状态y了(这个y是相对于当前行的下一行为x的时候,当前行就是y了)
#include<bits/stdc++.h>
using namespace std;
#define cl(a,b) memset(a,b,sizeof(a))
#define LL long long
#define pb push_back
#define gcd __gcd
#define For(i,j,k) for(int i=(j);i<k;i++)
#define lowbit(i) (i&(-i))
#define _(x) printf("%d\n",x)
const int maxn = 1e5+10;
const int inf = 1 << 28;
int M[1<<7][1<<7];
int k,n,m;
void dfs(int x,int y,int col){
if(col==k){
M[x][y]++;return ;
}
dfs(x<<1,y<<1|1,col+1);
dfs(x<<1|1,y<<1,col+1);
if(col+2<=k){
dfs((x<<2)+3,(y<<2)+3,col+2);//二进制11==3
}
}
int c[1<<7][1<<7];
void mul(int a[1<<7][1<<7],int b[1<<7][1<<7]){
cl(c,0);
for(int i=0;i<m;i++){
for(int j=0;j<m;j++){
for(int k=0;k<m;k++){
c[i][j]+=a[i][k]*b[k][j];c[i][j]%=12357;
}
}
}
memcpy(a,c,sizeof(c));
}
int a[1<<7][1<<7];
int main(){
while(~scanf("%d%d",&k,&n)){
LL N = 1LL*k*n;
if(N%2!=0){
puts("0");continue;
}
m=1<<k;
cl(M,0);
dfs(0,0,0);
for(int i=0;i<m;i++)a[i][i]=1;
while(n){
if(n&1)mul(a,M);
mul(M,M);
n>>=1;
}
printf("%d\n",a[m-1][m-1]);
}
return 0;
}