石头游戏的规则是这样的。
石头游戏在一个n行m列的方格阵上进行。每个格子对应了一个编号在0~9之间的操作序列。
操作序列是一个长度不超过6且循环执行、每秒执行一个字符的字符串。它包括：
数字0~9：拿0~9个石头到该格子。
NWSE：把这个格子内所有的石头推到相邻的格子。
D：拿走这个格子的石头。
石头游戏的问题是：当这个石头游戏进行了t秒之后，所有方格中最多的格子有多少个石头。
注意：所有格子的操作同时执行。
标签：多少 log i++ 矩阵 循环 tput 16px 表示 整数
题解：跟沼泽鳄鱼那道题差不多，只不过周期变成了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;
}
标签：多少 log i++ 矩阵 循环 tput 16px 表示 整数
原文地址：http://www.cnblogs.com/CQzhangyu/p/6791205.html
