ZOJ 3497 Mistwald

矩阵快速幂。

邻接矩阵的$P$次方就是走$P$步之后的方案数，这里只记录能否走到就可以了。然后再判断一下三种情况即可。

#pragma comment(linker, "/STACK:1024000000,1024000000")
#include<cstdio>
#include<cstring>
#include<cmath>
#include<algorithm>
#include<vector>
#include<map>
#include<set>
#include<queue>
#include<stack>
#include<ctime>
#include<iostream>
using namespace std;
typedef long long LL;
const double pi=acos(-1.0);
void File()
{
    freopen("D:\\in.txt","r",stdin);
    freopen("D:\\out.txt","w",stdout);
}
template <class T>
inline void read(T &x)
{
    char c = getchar();
    x = ;
    while(!isdigit(c)) c = getchar();
    while(isdigit(c))
    {
        x = x *  + c - '';
        c = getchar();
    }
}

int T;
int n,m,a,b,c,d,e,f,g,h;
char s[];
int r[][];
long long p;

struct Matrix
{
    int A[][];
    int R, C;
    Matrix operator*(Matrix b);
};

Matrix X, Y, Z;

Matrix Matrix::operator*(Matrix b)
{
    Matrix c;
    memset(c.A, , sizeof(c.A));
    int i, j, k;
    for (i = ; i <= R; i++)
        for (j = ; j <= b.C; j++)
            for (k = ; k <= C; k++)
                if(c.A[i][j]==) c.A[i][j]=(A[i][k]&b.A[k][j]);
    c.R = R; c.C = b.C;
    return c;
}

void init()
{
    memset(X.A, , sizeof X.A);
    memset(Y.A, , sizeof Y.A);
    memset(Z.A, , sizeof Z.A);

    for(int i=;i<=n*m;i++) Y.A[i][i]=;
    Y.R=Y.C=n*m;

    for(int i=;i<=n*m;i++)
        for(int j=;j<=n*m;j++) X.A[i][j]=r[i][j];
    X.R=X.C=n*m;
}

void work()
{
    while (p)
    {
        if (p %  == ) Y = Y*X;
        p = p >> ;
        X = X*X;
    }
}

int main()
{
    scanf("%d",&T);
    while(T--)
    {
        scanf("%d%d",&n,&m);

        memset(r,,sizeof r);
        for(int i=;i<=n;i++)
        {
            for(int j=;j<=m;j++)
            {
                scanf(" ((%d,%d),(%d,%d),(%d,%d),(%d,%d))",&a,&b,&c,&d,&e,&f,&g,&h);
                if(i==n&&j==m) continue;
                r[(i-)*m+j][(a-)*m+b]=;
                r[(i-)*m+j][(c-)*m+d]=;
                r[(i-)*m+j][(e-)*m+f]=;
                r[(i-)*m+j][(g-)*m+h]=;
            }
        }

        int Q; scanf("%d",&Q);
        while(Q--)
        {
            scanf("%lld",&p);
            init();
            work();

            if(Y.A[][n*m]==) printf("False\n");
            else
            {
                bool Find=;
                for(int i=;i<=n*m-;i++) if(Y.A[][i]==) Find=;
                if(Find==) printf("True\n");
                else printf("Maybe\n");
            }
        }
        printf("\n");
    }
    return ;
}