首页 > 代码库 > 【BZOJ2973】石头游戏 矩阵乘法

【BZOJ2973】石头游戏 矩阵乘法

2024-09-18 15:14:50   218人阅读

【BZOJ2973】石头游戏

Description

石头游戏的规则是这样的。
石头游戏在一个n行m列的方格阵上进行。每个格子对应了一个编号在0~9之间的操作序列。
操作序列是一个长度不超过6且循环执行、每秒执行一个字符的字符串。它包括：
数字0~9：拿0~9个石头到该格子。
NWSE：把这个格子内所有的石头推到相邻的格子。
D：拿走这个格子的石头。
石头游戏的问题是：当这个石头游戏进行了t秒之后，所有方格中最多的格子有多少个石头。
注意：所有格子的操作同时执行。

Input

第一行三个整数n, m, t, act。
接下来n行，每行m个字符，表示每个格子对应的操作序列。
最后act行，每行一个字符串，表示从0开始的每个操作序列。

Output

一个整数：游戏进行了t秒之后，所有方格中最多的格子有多少个石头。

Sample Input

1 6 10 3
011112
1E
E
0

Sample Output

3
【样例解释】
这是另一个类似于传送带的结构。左边的设备0间隔地产生石头并向东传送。设备1向右传送，直到设备2。10秒后，总共产生了5个石头，2个在传送带上，3个在最右边。
【数据约定】
0≤n, m≤8。
1≤act≤10。
1≤t≤10^9。

题解：跟沼泽鳄鱼那道题差不多，只不过周期变成了60。那我们构造60个转移矩阵f[0..59]，它们的乘积为x。那么ans=x^(T/12)*f[0..T%12]。至于怎么构造矩阵。。自己yy一下就好了吧？

又一次get了惨痛的历史教训，矩乘不满足交换律，所以一定要先ans*=x^(T/12)再ans*=f[0..T%12]。

#include <cstdio>#include <cstring>#include <iostream>#define P(A,B) ((A-1)*m+B)using namespace std;typedef long long ll;typedef struct matrix{	ll v[100][100];}M;M ans,x,emp,f[65];int n,m,tim,act;int len[20];ll maxx;char A[20][20],B[20][20],ch;M mmul(M a,M b){	M c=emp;	int i,j,k;	for(i=0;i<=n*m;i++)		for(j=0;j<=n*m;j++)			for(k=0;k<=n*m;k++)				c.v[i][j]+=a.v[i][k]*b.v[k][j];	return c;}void pm(int y){	while(y)	{		if(y&1)	ans=mmul(ans,x);		x=mmul(x,x),y>>=1;	}}int main(){	scanf("%d%d%d%d",&n,&m,&tim,&act);	int i,j,k,l,a,b;	for(i=1;i<=n;i++)	scanf("%s",A[i]);	for(i=0;i<act;i++)	scanf("%s",B[i]),len[i]=strlen(B[i]);	for(i=0;i<=n*m;i++)	x.v[i][i]=1;	ans.v[0][0]=1;	for(l=0;l<60;l++)	{		f[l].v[0][0]=1;		for(i=1;i<=n;i++)		{			for(j=1;j<=m;j++)			{				k=A[i][j-1]-‘0‘,ch=B[k][l%len[k]];				if(ch>=‘0‘&&ch<=‘9‘)	f[l].v[P(i,j)][P(i,j)]=1,f[l].v[0][P(i,j)]=ch-‘0‘;				if(ch==‘N‘&&i>1)	f[l].v[P(i,j)][P(i-1,j)]=1;				if(ch==‘W‘&&j>1)	f[l].v[P(i,j)][P(i,j-1)]=1;				if(ch==‘S‘&&i<n)	f[l].v[P(i,j)][P(i+1,j)]=1;				if(ch==‘E‘&&j<m)	f[l].v[P(i,j)][P(i,j+1)]=1;			}		}		x=mmul(x,f[l]);	}	pm(tim/60);	for(i=0;i<tim%60;i++)	ans=mmul(ans,f[i]);	for(i=1;i<=n*m;i++)	maxx=max(maxx,ans.v[0][i]);	printf("%lld",maxx);	return 0;}

【BZOJ2973】石头游戏 矩阵乘法


联系
我们