题意:给出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;
}