hdu 2807 The Shortest Path(矩阵优化)

题目链接：hdu 2807 The Shortest Path

题意：

如果矩阵A*B=C，那么就表示A-->B有一条单向路径，距离为1.

给一些矩阵，然后问任意两个矩阵直接的距离。

题解：

把矩阵读进来后处理一下邻接矩阵，然后Folyd一下。

#include<cstdio>
#include<algorithm>
#include<cstring>
#define mst(a,b) memset(a,b,sizeof(a))
#define F(i,a,b) for(int i=a;i<=b;++i)

using namespace std;
typedef long long ll;

const int mat_N=100,INF=1e9;//矩阵阶数，取膜
int n,m,q,mp[100][100];
struct mat{
    int c[mat_N][mat_N];
    void init(){mst(c,0);}
    mat operator*(mat b){
        mat M;int N=m-1;M.init();
        F(i,0,N)F(j,0,N)if(c[i][j])F(k,0,N)M.c[i][k]=M.c[i][k]+c[i][j]*b.c[j][k];
        return M;
    }
}A[91],tmp;

int check(mat &a,mat &b)
{
    F(i,0,m-1)F(j,0,m-1)if(a.c[i][j]!=b.c[i][j])return 0;
    return 1;
}

int main(){
    while(scanf("%d%d",&n,&m),n+m)
    {
        F(i,1,n)F(j,0,m-1)F(k,0,m-1)scanf("%d",&A[i].c[j][k]);
        F(i,1,n)F(j,1,n)mp[i][j]=INF;
        F(i,1,n)F(j,1,n)if(i!=j)
        {
            tmp=A[i]*A[j];
            F(k,1,n)if(j!=k&&k!=i)
            if(check(tmp,A[k]))mp[i][k]=1;
        }
        F(k,1,n)F(i,1,n)if(mp[i][k]!=INF)F(j,1,n)if(mp[k][j]!=INF)mp[i][j]=min(mp[i][j],mp[i][k]+mp[k][j]);
        scanf("%d",&q);
        while(q--)
        {
            int x,y;
            scanf("%d%d",&x,&y);
            mp[x][y]==INF?puts("Sorry"):printf("%d\n",mp[x][y]);
        }
    }
    return 0;
}