hdu 2807 矩阵比较+floyd

题意还是比较好理解的，即如果矩阵A、B、C满足A*B=C，则代表矩阵A的第i个城市与代表矩阵B的第j个城市之间存在通路。。。

orz。。。然后我就按一般的矩阵相等的条件去判断。。。结果TLE了。。。2000MS+了。。。

然后就学了一下矩阵的优化。。。优化后。。。60MS。。orz。。。

一开始的代码：

#include<iostream>
const int N=100;
const int inf=10000000;
using namespace std;

int edge[N][N];
int n,m;
struct Matrix{
    int map[N][N];
};
Matrix matrix[N],M;

Matrix calculate(const int &p,const int &q){
    Matrix Mat;
    for(int i=1;i<=m;i++){
        for(int j=1;j<=m;j++){
            int ans=0;
            for(int k=1;k<=m;k++){
                ans+=matrix[p].map[i][k]*matrix[q].map[k][j];
            }
            Mat.map[i][j]=ans;
        }
    }
    return Mat;
}

int cmp(const Matrix &p,const Matrix &q){
    for(int i=1;i<=m;i++){
        for(int j=1;j<=m;j++){
            if(p.map[i][j]!=q.map[i][j])
                return 0;
        }
    }
    return 1;
}

void floyd(){
    for(int k=1;k<=n;k++){
        for(int i=1;i<=n;i++){
            //if(i==k||edge[i][k]==inf)continue;
            for(int j=1;j<=n;j++){
                //if(i==k||j==k||edge[j][k]==inf)continue;
                if(edge[i][j]>edge[i][k]+edge[k][j]){
                    edge[i][j]=edge[i][k]+edge[k][j];
                }
            }
        }
    }
}


int main(){
    while(scanf("%d%d",&n,&m)!=EOF&&n&&m){
        for(int i=1;i<=n;i++){
            edge[i][i]=0;
            for(int j=i+1;j<=n;j++){
                edge[i][j]=edge[j][i]=inf;
            }
        }
        for(int i=1;i<=n;i++){
            for(int j=1;j<=m;j++){
                for(int k=1;k<=m;k++){
                    scanf("%d",&matrix[i].map[j][k]);
                }
            }
        }
        for(int i=1;i<=n;i++){
            for(int j=1;j<=n;j++){
                if(i==j)continue;
                M=calculate(i,j);
                for(int k=1;k<=n;k++){
                    if(i==k||j==k)continue;
                    if(cmp(M,matrix[k])){
                        edge[i][k]=1;
                    }
                }
            }
        }
        floyd();
        int t;
        scanf("%d",&t);
        while(t--){
            int v0,v;
            scanf("%d%d",&v0,&v);
            if(edge[v0][v]<inf){
                printf("%d\n",edge[v0][v]);
            }else
                printf("Sorry\n");
        }
    }
    return 0;
}

优化后的代码。。

#include<iostream>
const int N=100;
const int inf=10000000;
using namespace std;

int edge[N][N];
int n,m;
struct Matrix{
    int map[N][N];
};
struct ReMatrix{
    int Map[N];
};

Matrix matrix[N];
ReMatrix rematrix[N];

void cmp(const int &a,const int &b){
　//判断矩阵是否相等
    for(int i=1;i<=m;i++){
        if(rematrix[0].Map[i]!=rematrix[b].Map[i])
            return ;
    }
    edge[a][b]=1;
}

void CreateMap(const int &p,const int &q){
    for(int i=1;i<=m;i++){
        rematrix[0].Map[i]=0;
        for(int j=1;j<=m;j++){
            rematrix[0].Map[i]+=matrix[p].map[i][j]*rematrix[q].Map[j];  //就是两个矩阵相乘。。
        }
    }
    for(int i=1;i<=n;i++){
        if(i!=p&&i!=q)cmp(p,i);
    }
}


void floyd(){
    for(int k=1;k<=n;k++){
        for(int i=1;i<=n;i++){
            //if(i==k||edge[i][k]==inf)continue; //加了就ＷＡ了
            for(int j=1;j<=n;j++){
                //if(j==k||edge[j][k]==inf)continue;
                if(edge[i][j]>edge[i][k]+edge[k][j]){
                    edge[i][j]=edge[i][k]+edge[k][j];
                }
            }
        }
    }
}

int main(){
    while(scanf("%d%d",&n,&m)!=EOF&&n&&m){
        for(int i=1;i<=n;i++){
            edge[i][i]=0;
            for(int j=i+1;j<=n;j++){
                edge[i][j]=edge[j][i]=inf;
            }
        }
        for(int i=1;i<=n;i++){
            for(int j=1;j<=m;j++){
                rematrix[i].Map[j]=0;
                for(int k=1;k<=m;k++){
                    scanf("%d",&matrix[i].map[j][k]);
                    rematrix[i].Map[j]+=matrix[i].map[j][k]*k;  //化二维数组为一维数组，矩阵优化
                }
            }
        }
        for(int i=1;i<=n;i++){
            for(int j=1;j<=n;j++){
                if(i==j)continue;
                else CreateMap(i,j);
            }
        }
        floyd();
        int t;
        scanf("%d",&t);
        while(t--){
            int v0,v;
            scanf("%d%d",&v0,&v);
            if(edge[v0][v]<inf){
                printf("%d\n",edge[v0][v]);
            }else
                printf("Sorry\n");
        }
    }
    return 0;
}