ZOJ - 3497 Mistwald（离散数学）

题意：一张n*m的图，告诉你每一个点可以到达哪4个点，然后给出一个p，问从（1,1）出发P步之后，是否能到达（n，m），如果一定能达，输出true，如果可能到达输出maybe，如果不能到达，输出false

思路：我们可以给每一个点编号0 - n * m - 1，再建立一个邻接矩阵A，这样就是求0是否能到达n * m - 1，离散数学曾学过这样一个性质，A的p方得到的邻接矩阵，aij表示从i到j有多少长度为p的通路。所以我们可以通过矩阵快速幂求A的p次方来判断，注意题目说到达终点后就不能再移动，所以邻接矩阵的最后一行都为0。

#include <iostream>
#include <stdio.h>
#include <vector>

using namespace std;
typedef vector<int> vec;  
typedef vector<vec> mat; 

mat mul(mat &A,mat &B){
	mat C(A.size(),vec(B[0].size()));
	for(int i = 0; i < A.size(); i++){
		for(int j = 0; j < B[0].size(); j++){
			for(int k = 0; k < B.size(); k++){
				C[i][j] = C[i][j] + A[i][k] * B[k][j];
			}
		}
	}
	return C;
}

mat pow(mat A,int n){
	mat B(A.size(),vec(A.size()));
	for(int i = 0; i < A.size(); i++){
		B[i][i] = 1;
	}
	
	while(n){
		if(n & 1){
			B = mul(B,A);
		}
		A = mul(A,A);
		n >>= 1;
	}
	return B;
}

int main(void){
	int t;
	scanf("%d",&t);
	
	while(t--){
		int n,m;
		scanf("%d%d",&n,&m);
		getchar();
		mat C(n * m,vec(n * m));
		int a,b,c,d,e,f,g,h;
		for(int i = 0; i < n; i++){
			for(int j = 0; j < m; j++){
				scanf("((%d,%d),(%d,%d),(%d,%d),(%d,%d))",&a,&b,&c,&d,&e,&f,&g,&h);
				a--;b--,c--;d--;e--;f--;g--;h--;
				getchar();
			    int temp = i * m + j;
				if(temp == n * m - 1){
					continue;
				}
				C[temp][a * m + b] = 1;
				C[temp][c * m + d] = 1;
				C[temp][e * m + f] = 1;
				C[temp][g * m + h] = 1;
			}
		}
		int q;
		scanf("%d",&q);
		while(q--){
			int p;
			scanf("%d",&p);
			mat ans = pow(C,p);
			if(ans[0][n * m - 1] == 0){
				printf("False\n");
			}
			else{
				int flag = 0;
				for(int i = 0; i < n * m - 1; i++){
					if(ans[0][i]){
						flag = 1;
						break;
					}
				}
				if(flag){
					printf("Maybe\n");
				}
				else{
					printf("True\n");
				}
			}
		}
		printf("\n");
	}
	return 0;
}