HIHO #1162 : 骨牌覆盖问题·三

题目链接

这题和之前的一个 状态压缩二，是很类似的，这道题的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;
}