七十分超时代码，具体思路可以仿照acwing蒙德里安的梦想
就是枚举第i层放 a ，第i-1层放b的合法方案其中b和a都是同一块；

	#include <iostream>
#include <algorithm>
#include <cstring>
#include <sstream>
using namespace std;
typedef long long LL;
const int N = 5e5 + 4, mod = 1e9 + 7, M = 1<<7;

LL n;
int  m;
int dp[N][M];


void dfs(int a, int b, int p)
{
	if(p == m) v[a][b]++, c[a][b]++;
	else
	{
		if(p > 0)
		{
			if(!((a>>p&1)||(b>>p&1)||(b>>p-1&1)))
				dfs(a|1<<p, b|1<<p|1<<p-1, p+1);
			if(!((a>>p&1)||(b>>p&1)||(a>>p-1&1)))
				dfs(a|1<<p|1<<p-1, b|1<<p, p+1);
		}
		if(p+1 != m)
		{
			if(!((a>>p&1)||(b>>p&1)||(b>>p+1&1)))
				dfs(a|1<<p, b|1<<p|1<<p+1, p+1);
			if(!((a>>p&1)||(b>>p&1)||(a>>p+1&1)))
				dfs(a|1<<p|1<<p+1, b|1<<p, p+1);
		}
		dfs(a, b, p+1);
	}
}

int main() 
{
	scanf("%lld%d", &n, &m);
	dfs(0, 0, 0);
	for(int i = 0;i < 1<<m;i ++)
		dp[1][i] = v[i][(1<<m)-1];
	for(int i = 2;i < n;i++)
		for(int j = 0;j < 1<<m;j ++)
			for(int k = 0;k < 1<<m;k ++)
				dp[i][j] = (dp[i][j] + (LL)dp[i-1][((1<<m)-1)^k]*v[j][k])%mod;
	
	cout<<dp[n-1][(1<<m)-1];
	return 0;
}

100 分矩阵快速幂

#include <iostream>
#include <algorithm>
#include <cstring>
#include <sstream>
using namespace std;
typedef long long LL;
const int N = 5e5 + 4, mod = 1e9 + 7, M = 1<<7;

LL n;
int  m, c[M][M], tem[M][M], res[M][M];
int dp[3][M];

void multi(int a[][M], int b[][M], int n)
{
	memset(tem, 0, sizeof tem);	
	for(int i = 0;i < n;i ++)
		for(int j = 0;j < n;j ++)
			for(int k = 0;k < n;k ++)
				tem[i][j] = (tem[i][j] + (LL)a[i][k]*b[k][j])%mod;
	for(int i = 0;i < n;i ++)
		for(int j = 0;j < n;j ++)
			a[i][j] = tem[i][j];
}

void pow(int a[][M], LL n)
{
	memset(res, 0, sizeof res);
	for(int i = 0;i < 1<<m;i ++)
		res[i][i] = 1;
	while(n)
	{
		if(n&1) multi(res, a, 1<<m);
		multi(a, a, 1<<m);
		n>>=1;
	}
}

void dfs(int a, int b, int p)
{
	if(p == m) c[a][((1<<m)-1)^b]++;   // 方便 矩阵快速幂 这题c[i][j]代表第二层是i 第一层是i的合法方案的其他空缺； 
	else
	{
		if(p > 0)
		{
			if(!((a>>p&1)||(b>>p&1)||(b>>p-1&1)))
				dfs(a|1<<p, b|1<<p|1<<p-1, p+1);
			if(!((a>>p&1)||(b>>p&1)||(a>>p-1&1)))
				dfs(a|1<<p|1<<p-1, b|1<<p, p+1);
		}
		if(p+1 != m)
		{
			if(!((a>>p&1)||(b>>p&1)||(b>>p+1&1)))
				dfs(a|1<<p, b|1<<p|1<<p+1, p+1);
			if(!((a>>p&1)||(b>>p&1)||(a>>p+1&1)))
				dfs(a|1<<p|1<<p+1, b|1<<p, p+1);
		}
		dfs(a, b, p+1);
	}
}


int main() 
{
	scanf("%lld%d", &n, &m);
	dfs(0, 0, 0);
	for(int i = 0;i < 1<<m;i ++)
		dp[1][i] = c[i][0];
		
	if(n >= 2)pow(c, n-2); 
	
	for(int i = 0;i < 1<<m;i ++)
		dp[2][(1<<m)-1] = (dp[2][(1<<m)-1] + (LL)res[(1<<m)-1][i]*dp[1][i])%mod;
		
	cout<<dp[2][(1<<m)-1];
	return 0;
}