【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;
}