hdu2807（floyd +快速矩阵计算）

来源：互联网发布：知乎南京装修工作室编辑：程序博客网时间：2024/06/05 09:28

题意：给出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 0xffffffusing 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;}

0 0