标签：tem 拼图 201409 ++ dfs int long include

七十分超时代码，具体思路可以仿照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;
}

标签：tem,拼图,201409,++,dfs,int,long,include
来源： https://blog.csdn.net/weixin_51942413/article/details/121196891

本站声明： 1. iCode9 技术分享网（下文简称本站）提供的所有内容，仅供技术学习、探讨和分享；
2. 关于本站的所有留言、评论、转载及引用，纯属内容发起人的个人观点，与本站观点和立场无关；
3. 关于本站的所有言论和文字，纯属内容发起人的个人观点，与本站观点和立场无关；
4. 本站文章均是网友提供，不完全保证技术分享内容的完整性、准确性、时效性、风险性和版权归属；如您发现该文章侵犯了您的权益，可联系我们第一时间进行删除；
5. 本站为非盈利性的个人网站，所有内容不会用来进行牟利，也不会利用任何形式的广告来间接获益，纯粹是为了广大技术爱好者提供技术内容和技术思想的分享性交流网站。