hbcpc problem J 骨牌游戏（状压dp+矩阵快速幂）

已于 2023-05-24 23:54:08 修改
阅读量213 收藏 4
点赞数 3
文章标签： 游戏 算法
于 2023-05-24 23:51:29 首次发布
版权声明：本文为博主原创文章，遵循 CC 4.0 BY-SA 版权协议，转载请附上原文出处链接和本声明。
本文链接：https://blog.csdn.net/m0_74315028/article/details/130819803
版权
文章讲述了如何使用矩阵快速幂和状态压缩动态规划方法解决覆盖特定大小棋盘的问题。作者最初尝试了简单的状态压缩dp，但因数据范围过大导致超时，然后通过观察斐波那契数列的规律，结合矩阵快速幂优化算法，最终成功解决了当n为奇数时，覆盖n*m棋盘的方案数问题。
摘要由CSDN通过智能技术生成

题目要求：

用1*2的骨牌完全覆盖n*m棋盘，问有多少种方案

1<=n<=5,1<=m<=10^18

答案998244353取模

思路

你可以当我虚线下面是废话，我想了两天终于想出来了，我是菜鸡

总的来说就是矩阵快速幂加状态压缩dp

#include<bits/stdc++.h> 
#define int long long
using namespace std;
const int N=5,M=1<<N,mod=998244353;
int n,m;
bool st[M];
int A[M][M],tmp[M][M],res[M][M];
 void MXMP(int a[][M], int b[][M])
{
	memset(tmp, 0, sizeof tmp);//tmp用来存放相乘的结果

	for (int i = 0; i < (1<<n); i++) {
		for (int j = 0; j < (1<<n); j++) {
			for (int k = 0; k < (1<<n); k++) {
				tmp[i][j] = (tmp[i][j] + (a[i][k] * b[k][j]) % mod) % mod;//记得摸mod
			}
		}
	}

	for (int i = 0; i < (1<<n); i++) {
		for (int j = 0; j < (1<<n); j++) {
			a[i][j] = tmp[i][j];//记得把tmp赋值给a
		}
	}
}

void powermod(int A[][M], int x)
{
	memset(res, 0, sizeof res);//res存放结果
	for (int i = 0; i < (1<<n); i++) {
		for (int j = 0; j < (1<<n); j++) {
			res[i][j] = A[i][j];
		}
	}
	while (x) {
		if (x & 1) MXMP(res, A);//x为奇数就将A与res相乘
		MXMP(A, A);//A自身相乘
		x >>= 1;//x除以二
	}

	for (int i = 0; i < (1<<n); i++) {//将结果赋值到A
		for (int j = 0; j < (1<<n); j++) {
			A[i][j] = res[i][j];
		}
	}
}
signed main()
{
	ios_base::sync_with_stdio(0); cin.tie(0),cout.tie(0);
	cin>>n>>m;
	if((n&1)&&(m&1)){//判断是否都是奇数
		cout<<0;
		return 0;
	}
	memset(st,true,sizeof st);//记录是否可以放
	for(int i=0;i<(1<<n);i++){
		int cnt=0;
		for(int j=0;j<n;j++){
			if((i>>j)&1){
				if(cnt&1){
					st[i]=false;
					break;
				}
			}else cnt++;
		}
		if(cnt&1) st[i]=false;
	}
	memset(A,0,sizeof A);
	for(int j=0;j<(1<<n);j++){
		for(int k=0;k<(1<<n);k++){
			if((j&k)==0&&st[j|k]){
				A[j][k]=1; //记录是由哪个状态可以继承过来
			}
		}
	}
	m-=1;//因为res是由A赋值过去的所以res*A会多算一次得减一
	powermod(A,m);//矩阵快速幂
	cout<<A[0][0];
	/*for(int i=0;i<(1<<n);i++){
		for(int j=0;j<(1<<n);j++){
			cout<<A[i][j]<<' ';
		}
		cout<<'\n';
	}*/
	
}

---------------------------------------------------------------------------------------------------------------------------------

原本没看范围以为就是一道简单的状态压缩dp很开心的按版子打上去交了一发就超时了

后面一看范围10的18次方,比赛时也是没有做出来

这是状态压缩dp的代码（超时）

#include<bits/stdc++.h> 
#define int long long
using namespace std;
const int N=5,M=1<<N;
int n,m;
bool st[M];
int f[M][M];
signed main()
{
	ios_base::sync_with_stdio(0); cin.tie(0),cout.tie(0);
	cin>>n>>m;
	memset(st,true,sizeof st);
	for(int i=0;i<(1<<n);i++){
		int cnt=0;
		for(int j=0;j<n;j++){
			if((i>>j)&1){
				if(cnt&1){
					st[i]=false;
					break;
				}
			}else cnt++;
		}
		if(cnt&1) st[i]=false;
	}
	memset(f,0,sizeof f);
	f[0][0]=1;
	for(int i=1;i<=m;i++){
		for(int j=0;j<(1<<n);j++){
			for(int k=0;k<(1<<n);k++){
				if((j&k)==0&&st[j|k]){
					f[i][j]+=f[i-1][k];
				}
			}
		}
	}
	cout<<f[m][0];
}

后面我们还是找了一下规律当n=1 就是0,1,0,1,0,1

n=2就是1,2,3,5,8,13

发现是一个斐波那契数列,这时候m可以取到10^18正常求肯定会超时，但是斐波那契还可以用矩阵快速幂来求，这时候我就有大致思路，是不是整体都可以用矩阵快速幂来优化

矩阵快速幂

下面是求斐波那契的代码

#include<bits/stdc++.h>
#define int long long
using namespace std;
const int mod = 998244353;
int A[2][2] = { {1,1},{1,0} };
int res[2][2], n;
int tmp[2][2];

void MXMP(int a[][2], int b[][2])
{
	memset(tmp, 0, sizeof tmp);

	for (int i = 0; i < 2; i++) {
		for (int j = 0; j < 2; j++) {
			for (int k = 0; k < 2; k++) {
				tmp[i][j] = (tmp[i][j]%mod + (a[i][k] * b[k][j]) % mod) % mod;
			}
		}
	}

	for (int i = 0; i < 2; i++) {
		for (int j = 0; j < 2; j++) {
			a[i][j] = tmp[i][j]%mod;
		}
	}
}

void powermod(int A[][2], int x)
{
	memset(res, 0, sizeof res);
	for (int i = 0; i < 2; i++) {
		for (int j = 0; j < 2; j++) {
			res[i][j] = 1;
		}
	}
	while (x) {
		if (x & 1) MXMP(res, A);
		MXMP(A, A);
		x >>= 1;
	}

	for (int i = 0; i < 2; i++) {
		for (int j = 0; j < 2; j++) {
			A[i][j] = res[i][j];
		}
	}
}
signed main()
{
	ios_base::sync_with_stdio(0); cin.tie(0), cout.tie(0);
	cin >> n;
	n-=1;
	if(n==0||n==1){
		cout<<1;
		return 0;
	}
	powermod(A, n - 2);
	int q = (A[0][1] + A[0][0])%mod;
	cout << q;
}

然后我n=3是1，0,3,0,11,0

我想这先找一下偶数的规律

然后去找了n=4  1,1,5,11,36,95,281

 看他数出现的规律然后去推断是由哪个数继承过来的

我把出现过数字的号码记下来了

 看上图1可以加到12456进去，2可以加到12进去

后面想到可以用矩阵就是右边那图

可以写成n==4的代码

#include<bits/stdc++.h>
#define int long long
using namespace std;
const int mod = 998244353;
int A[6][6] = { {1,1,0,1,1,1},{1,1,0,0,0,0},{0,0,0,1,0,0},{1,0,1,0,0,0},{1,0,0,0,1,0},{1,0,0,0,0,0}};
int res[6][6], n;
int tmp[6][6];

void MXMP(int a[][6], int b[][6])
{
	memset(tmp, 0, sizeof tmp);

	for (int i = 0; i < 6; i++) {
		for (int j = 0; j < 6; j++) {
			for (int k = 0; k < 6; k++) {
				tmp[i][j] = (tmp[i][j] + (a[i][k] * b[k][j]) % mod) % mod;
			}
		}
	}

	for (int i = 0; i < 6; i++) {
		for (int j = 0; j < 6; j++) {
			a[i][j] = tmp[i][j];
		}
	}
}

void powermod(int A[][6], int x)
{
	memset(res, 0, sizeof res);
	for (int i = 0; i < 6; i++) {
		for (int j = 0; j < 6; j++) {
			res[i][j] = A[i][j];
		}
	}
	while (x) {
		if (x & 1) MXMP(res, A);
		MXMP(A, A);
		x >>= 1;
	}

	for (int i = 0; i < 6; i++) {
		for (int j = 0; j < 6; j++) {
			A[i][j] = res[i][j];
		}
	}
}
signed main()
{
	ios_base::sync_with_stdio(0); cin.tie(0), cout.tie(0);
	cin >> n;
	n-=2;
	if(n==-1){
		cout<<1;
		return 0;
	} 
	powermod(A, n);
	int q = A[0][0];
	cout << q;
}

做到这里突然就悟了，就是打一张表看能不能继承在用一下快速幂就行了

卧槽我太聪明了

#include<bits/stdc++.h> 
#define int long long
using namespace std;
const int N=5,M=1<<N,mod=998244353;
int n,m;
bool st[M];
int A[M][M],tmp[M][M],res[M][M];
 void MXMP(int a[][M], int b[][M])
{
	memset(tmp, 0, sizeof tmp);//tmp用来存放相乘的结果

	for (int i = 0; i < (1<<n); i++) {
		for (int j = 0; j < (1<<n); j++) {
			for (int k = 0; k < (1<<n); k++) {
				tmp[i][j] = (tmp[i][j] + (a[i][k] * b[k][j]) % mod) % mod;//记得摸mod
			}
		}
	}

	for (int i = 0; i < (1<<n); i++) {
		for (int j = 0; j < (1<<n); j++) {
			a[i][j] = tmp[i][j];//记得把tmp赋值给a
		}
	}
}

void powermod(int A[][M], int x)
{
	memset(res, 0, sizeof res);//res存放结果
	for (int i = 0; i < (1<<n); i++) {
		for (int j = 0; j < (1<<n); j++) {
			res[i][j] = A[i][j];
		}
	}
	while (x) {
		if (x & 1) MXMP(res, A);//x为奇数就将A与res相乘
		MXMP(A, A);//A自身相乘
		x >>= 1;//x除以二
	}

	for (int i = 0; i < (1<<n); i++) {//将结果赋值到A
		for (int j = 0; j < (1<<n); j++) {
			A[i][j] = res[i][j];
		}
	}
}
signed main()
{
	ios_base::sync_with_stdio(0); cin.tie(0),cout.tie(0);
	cin>>n>>m;
	if((n&1)&&(m&1)){//判断是否都是奇数
		cout<<0;
		return 0;
	}
	memset(st,true,sizeof st);//记录是否可以放的了
	for(int i=0;i<(1<<n);i++){
		int cnt=0;
		for(int j=0;j<n;j++){
			if((i>>j)&1){
				if(cnt&1){
					st[i]=false;
					break;
				}
			}else cnt++;
		}
		if(cnt&1) st[i]=false;
	}
	memset(A,0,sizeof A);
	for(int j=0;j<(1<<n);j++){
		for(int k=0;k<(1<<n);k++){
			if((j&k)==0&&st[j|k]){
				A[j][k]=1; //记录是由哪个状态可以继承过来
			}
		}
	}
	m-=1;//因为res是由A赋值过去的所以res*A会多算一次得减一
	powermod(A,m);//矩阵快速幂
	cout<<A[0][0];
	/*for(int i=0;i<(1<<n);i++){
		for(int j=0;j<(1<<n);j++){
			cout<<A[i][j]<<' ';
		}
		cout<<'\n';
	}*/
	
}

阿根廷必胜
关注
  • 3
    点赞
  • 踩
  • 4
    收藏
    觉得还不错? 一键收藏
  • 0
    评论
棋盘覆盖 状压DP+矩阵快速幂
Ark
08-23 616
题意：有一个m 行n 列的矩形方格棋盘，1 &lt; = m&lt; = 5,1=&lt; n&lt; =10^9，用1*2 的骨牌(可横放或竖放)完全覆盖，骨牌不能重叠，有多少种不同的覆盖的方法。你只需要求出覆盖方法总数 mod p 的值即可。 看到1e9立马知道快速幂DP或者数学方法，然后m&lt;=5就状压吧 定义f[s][t]表示从s到t有多少种方案转移:则有f[s][t] = sig...
POJ 3420 Quad Tiling 状压DP+矩阵快速幂
看曙光
07-26 1292
链接：http://poj.org/problem?id=3420 题意：给一个4*N(1 ≤ N ≤ 1e9)的矩形空间，并且给不限块数的1*2的多米诺骨牌，问是由多少种方式能把这个矩形空间填满。 思路：看到这种问题果断想到状压，虽然是在看矩阵的时候看到的这道题。dp[i][j]表示在第i行状态为j的情况下的填满方式数，j的二进制表示中0表示对应位置上一行的骨牌是竖放，或者对应位置的骨牌是横
Jzoj P4787 数格子___状压dp+矩阵快速幂
生如夏花之绚烂
05-29 137
题目大意： 用2∗12*12∗1的骨牌覆盖4∗N4*N4∗N的矩形。 问方案数对MMM取模后的结果。 多组数据，不超过101010组 1≤N≤109，0≤M≤1051≤N≤10^9，0≤M≤10^51≤N≤109，0≤M≤105 分析： 令fi,sf_{i,s}fi,s​，表示第iii列放了多少个横着的骨牌时覆盖完前iii列的方案总数， sss为这些横着的骨牌的位置的状态 [i,1][i,1][i...
骨牌问题（DP+矩阵快速幂）
feng_zhiyu的博客
07-30 1361
在2*N的一个长方形方格中，用一个1*2的骨牌排满方格。 问有多少种不同的排列方法。例如：2 * 3的方格，共有3种不同的排法。（由于方案的数量巨大，只输出 Mod 10^9 + 7 的结果） Input 输入N(N <= 1000) Output 输出数量 Mod 10^9 + 7 Sample Input 3 Sample Output 3分析：从N小的时候入手，比如N=2
基于C#和Unity开发的米诺骨牌的小游戏源码(含项目说明+sln解决方案).zip
12-28
基于C#和Unity开发的米诺骨牌的小游戏源码(含项目说明+sln解决方案).zip 需求分析 1. 老师上课演示关于多米诺骨牌的那些关于重量、大小、重心位置的实验过于浪费时间，直接放视频的话，又体现不出操作性，这个时候...
Cut Throat v1.0 - 多米诺骨牌多人游戏 Unity项目源码
01-26
Cut Throat v1.0 - 多米诺骨牌多人游戏 Unity 支持Unity2020.3.20f1及以上 Cut Throat – Dominoes 是一款使用 Unity 3D 制作的离线和在线多人游戏。它包括apk、源文件和详细文档 设置： 在设置中您可以打开或关闭...
chickenfoot_vis:多米诺骨牌游戏模拟的网络可视化
07-11
鸡脚_vis 多米诺骨牌游戏模拟的网络可视化。入门确保您安装了最新的软件包 npm installbower install注意：如果你没有安装npm ，请确保你已经安装了 。 如果你没有凉亭， npm install -g bower 。 上述步骤将下载...
domino-mofo:Ruby 命令行多米诺骨牌游戏
06-11
Ruby 命令行多米诺骨牌游戏 这是我使用 Cucumber 和 Rspec 以 BDD 方式开发的一个有趣的命令板多米诺骨牌游戏。 玩 下载并运行： bin / domino_mofo 命令 铅 <多米诺骨牌> 用 <suit> 在 <domino> 上玩 ...
MyDomino42:我的多米诺骨牌42游戏
03-21
《我的多米诺骨牌42游戏》是一款基于C#编程语言开发的桌面游戏，旨在为玩家提供一个体验经典多米诺骨牌游戏42的平台。多米诺骨牌42是一种策略性极强的双人或团队游戏，源自美国，广受欢迎。通过这个项目，我们可以...
unity2D游戏开发16弹弓动画
最新发布
m0_58498702的博客
07-30 172
设置属性为2D Simple Directional，再点击加号选择Add Motion Field。创建这三个参数isWalking:Bool ,xDir:Float，yDir:Float。选择player-idle，右键MakeTransition,链接到Walk Tree。选中PlayerObject，打开Animator,删除原来的四个状态。双击Walk Tree查看Blend Tree Graph。将Any State到idle之间的转换，设置false。点击四个节点，勾选Loop Time。
unity2D游戏开发12单例
m0_58498702的博客
07-29 31
当一个类提供了游戏中其他几个类使用的功能时，单例会很有用，例如，在Game Manager 类中协调游戏逻辑，单例可以提供对该类及其功能的公共统一访问入口。因为现在是动态生成玩家，所以，Virtual Camera的follow属性没有实例，为了让摄像机跟随玩家走动，我们将创建一个类来设置Virtual Camera。给SpawnPoint添加SpawnPoint脚本，再将SpawnPoint拖到Prefabs文件夹生成预购件，再删除Hierarhy面板的SpawnPoint。
cocos creator 3学习记录01——如何替换图片
SilenceJude的博客
07-26 330
指的是 loader.loaderRes、loader.loadResDir、loader.release等loader开头的API，并不是 resources.load，因此无需担心。3、对于图片类型，直接访问路径所对应的图片地址，得到的图片类型会是ImageAsset，此时必须指定具体的子资源。2、调用resources.load API，路径是默认在resources文件夹下，所以直接写后面的路径即可。不过吧，这种方法我并不推荐，因为需要提前声明好属性，甚至需要提前关联好。
「Unity3D」场景中的距离单位Unit与相关设置PixelsToUnits、PixelsPerUnit
scott.cgi
07-25 435
这就涉及到一个问题，即如何按照比例显示像素？或说如何显示2D图片的大小？
Unity Yaml
记录
07-29 158
选择用于存储序列化资产的格式。
ScriptableObject使用
记录
07-26 613
结合两者特点，可将SO数据存储在JSON中，修改时，直接修改JSON文件;JSON和XML等文件格式可能很难在编辑器中修改，在Unity外，使用文本编辑都可轻松修改。多个对象共享一份So实例，多个对象共享的数据可提取出来作为so实例，减少重复，节约内存。So脚本可正常使用，但显示不正确，集合元素显示为“Type mismatch”。需要监听或触发该委托的对象可获取该So实例，监听者注册方法，广播者触发方法。，有点类似静态数据，同一个实例文件可被多个对象引用；此外可拓展So脚本，直接触发方法，方便测试。
Unity GameObject学习笔记
2402_83809362的博客
07-29 228
GameObject成员变量。GameObject成员变量。GameObject静态方法。
Unity--UI事件触发器--EventTrigger
qq_40669895的博客
07-29 599
Unity--UI事件触发器--EventTrigger
右手系转Unity左手系坐标系
薛文旺
07-26 374
【代码】右手系转Unity左手系坐标系。
python骨牌铺方格
09-28
然后我们可以使用一个循环从2到n，每次更新dp[i]的值为dp[i-1]+dp[i-2]，即长度为i的方格可以通过在长度为i-1的方格的右边加一个骨牌或者在长度为i-2的方格的右边加两个骨牌得到。最后，结果就是dp[n]。 以下是用...
评论
添加红包

请填写红包祝福语或标题

个

红包个数最小为10个

元

红包金额最低5元

当前余额3.43元 前往充值 >
需支付：10.00元
成就一亿技术人!
领取后你会自动成为博主和红包主的粉丝 规则
hope_wisdom
发出的红包
实付元
使用余额支付
点击重新获取
扫码支付
钱包余额 0

抵扣说明：

1.余额是钱包充值的虚拟货币，按照1:1的比例进行支付金额的抵扣。
2.余额无法直接购买下载，可以购买VIP、付费专栏及课程。

余额充值