Mistwald POJ

一开始看不出来是快速幂矩阵的题目

先要把整个地图离散化为1，2，3，4，。。。。

连成一个有向图

邻接矩阵的平方意为：假如a->b  且b->c     那么一次平方后   a->c    相当于floyd路径的连通

所以p次方就是   该矩阵经过p次幂     如果路径为1 则代表可以走

离散化   i*m+j  ； x*m+y        i，j，x，y必须是基于0~n-1 标准的！

#include <iostream>
#include <cstdio>
#include<cstring>

using namespace std;

typedef long long ll;

 struct Matrix{
    ll m[][];
};

int n,m;
int siz;

Matrix Mul(Matrix a, Matrix b)
{
    Matrix c;
    memset(c.m, , sizeof(c.m));
    for (int i = ; i < siz; i++)
    {
        for (int j = ; j < siz; j++)
        {
            for (int k = ; k < siz; k++)
            {
                c.m[i][j] = (c.m[i][j] + (a.m[i][k] * b.m[k][j])  ) ;
            }
        }
    }
    return c;
}

Matrix fastm(Matrix a, ll num)
{
    Matrix res;
    memset(res.m, , sizeof(res.m));
    for(int i=;i<siz;i++)
        res.m[i][i]=;
    while (num)
    {
        if (num & )
            res = Mul(res, a);
        num >>= ;
        a = Mul(a, a);
    }
    return res;
}

int main()
{
    int T;
    scanf("%d",&T);
    while(T--)
    {
        scanf("%d%d",&n,&m);
        getchar();
        Matrix a;
        memset(a.m,,sizeof(a.m));
        for(int i=;i<n;i++)
        {
            for(int j=;j<m;j++)
            {
                int x1,y1,x2,y2,x3,y3,x4,y4;
                scanf("((%d,%d),(%d,%d),(%d,%d),(%d,%d))",&x1,&y1,&x2,&y2,&x3,&y3,&x4,&y4);
                getchar();
                if(i==n-&&j==m-)  continue;//到达结尾就结束  所以不处理
                int now=i*m+j;
                a.m[now][(x1-)*m+y1-]=;
                a.m[now][(x2-)*m+y2-]=;
                a.m[now][(x3-)*m+y3-]=;
                a.m[now][(x4-)*m+y4-]=;
            }
        }
        siz=n*m;
        int Q;
        scanf("%d",&Q);
        while(Q--)
        {
            int p;
            scanf("%d",&p);
            Matrix res=fastm(a,p);
            int flag=;
            if(res.m[][siz-]==)
                printf("False\n");
            else
            {
                for(int i=;i<siz-;i++)
                {
                    if(res.m[][i]){
                        flag=;
                        break;
                    }
                }
                if(flag==) printf("Maybe\n");
                else    printf("True\n");
            }
        }
        printf("\n");
    }
    return ;
}