zoj 3497 Mistwald (矩阵应用)

/*
邻接矩阵A的k次幂得到的新矩阵B中，bij表示ai到aj长度为k的通路数。
上面的方法计算的是任意两点通路数，
而本题目要求终点不能有出边，
只要把终点那一行全部置为0就ok，
这里不求通路数，改为求能否k步到达终点。
*/
# include <stdio.h>
# include <algorithm>
# include <string.h>
# include <math.h>
using namespace std;
struct node
{
int m1[50][50];
};
int n,m;
struct node a,I;
node f(node a1,node a2)
{
int i,j,k;
node a3;
for(i=1; i<=n*m; i++)
{
for(j=1; j<=n*m; j++)
{
a3.m1[i][j]=0;
for(k=1; k<=n*m; k++)
a3.m1[i][j]+=(a1.m1[i][k]*a2.m1[k][j]);
}
}
return a3;
}

node quick(int x)
{
node mm=a;
node b=I;
while(x)
{
if(x&1)
b=f(b,mm);
mm=f(mm,mm);
x=x>>1;
}
return b;
}
int main()
{
int i,j,t,x1,x2,x3,x4,y1,y2,y3,y4;
while(~scanf("%d",&t))
{
while(t--)
{
scanf("%d%d",&n,&m);
getchar();
memset(a.m1,0,sizeof(a.m1));
memset(I.m1,0,sizeof(I.m1));//单位矩阵
for(i=1; i<=n*m; i++)
I.m1[i][i]=1;
for(i=1; i<=n; i++)
{
for(j=1; j<=m; j++)
{
scanf("((%d,%d),(%d,%d),(%d,%d),(%d,%d))",&x1,&y1,&x2,&y2,&x3,&y3,&x4,&y4);
getchar();
if((i-1)*m+j==n*m)
continue;
a.m1[(i-1)*m+j][(x1-1)*m+y1]=1;
a.m1[(i-1)*m+j][(x2-1)*m+y2]=1;
a.m1[(i-1)*m+j][(x3-1)*m+y3]=1;
a.m1[(i-1)*m+j][(x4-1)*m+y4]=1;
}
}
int q,k;
scanf("%d",&q);
while(q--)
{
scanf("%d",&k);
node b=quick(k);
if(b.m1[1][n*m]==0)
printf("False\n");
else
{

for(i=1; i<=n*m; i++)
{
if(b.m1[1][i])
break;
}
if(i==n*m)///只能到达终点
printf("True\n");
else///还可以到达其他点
printf("Maybe\n");
}
}
printf("\n");

}
}
return 0;
}