hdu2807（floyd +快速矩阵计算）

题意：给出n个m×m的矩阵，若矩阵a×b=c，则a与c之间距离为1，给出k个问题，求出最小步数，若没有，输出sorry。

思路：本身很简单的Floyd可以解决，但是在矩阵计算的时候将二维矩阵化为一位矩阵，具体操作方法是给不同列（或者行）的数一个权值，将行的总权值求出来，存为一位数组，之后比较数组就只要一个for循环就行了。

#include <iostream>
#include <cstdio>
#include <cstdlib>
#include <cstring>
#include <algorithm>
#include <cmath>
#define N 85
#define inf 0xffffff
using namespace std;
int n,m,map[N][N];
struct node
{
    int p[N][N];
}Node[N];
struct nope
{
    int d[N];
}Nope[N];
void calc(int a,int b)
{
    for(int i=0;i<m;i++)
    {
        Nope[82].d[i]=0;
        for(int j=0;j<m;j++)
        {
            Nope[82].d[i]+=Node[a].p[i][j]*Nope[b].d[j];
        }
    }
    for(int i=0;i<n;i++)
    {
        if(i==a ||i==b) continue;
        int j;
        for(j=0;j<m;j++)
            if(Nope[82].d[j]!=Nope[i].d[j]) break;
        if(j==m) map[a][i]=1;
    }
}
int main()
{
    int p,q;
    while(scanf("%d %d",&n,&m))
    {
        if(n==0 && m==0) break;
        for(int i=0;i<n;i++)
        {
            for(int j=0;j<n;j++)
            {
                if(i==j) map[i][j]=0;
                map[i][j]=inf;
            }
        }
        for(int i=0;i<n;i++)
        {
            for(int j=0;j<m;j++)
            {
                Nope[i].d[j]=0;
                for(int t=0;t<m;t++)
                {
                    scanf("%d",&Node[i].p[j][t]);
                    Nope[i].d[j]+=Node[i].p[j][t]*(t+1);
                }
            }
        }
        for(int i=0;i<n;i++)
            for(int j=0;j<n;j++)
        {
            if(i==j) continue;
            calc(i,j);
        }
        for(int t=0;t<n;t++)
            for(int i=0;i<n;i++)
                for(int j=0;j<n;j++)
        {
            if(map[i][t]+map[t][j]<map[i][j]) map[i][j]=map[i][t]+map[t][j];
        }
        int k;
        scanf("%d",&k);
        while(k--)
        {
            scanf("%d %d",&p,&q);
            if(map[p-1][q-1]!=inf) printf("%d\n",map[p-1][q-1]);
            else printf("Sorry\n");
        }
    }

    return 0;
}