【洛谷P1730】最小密度路径

题目大意：给定一个 N 个点，M 条边的有向图，现有 Q 个询问，每次询问 X 到 Y 的最小密度路径是多少。最小密度路径的定义是路径长度除以路径边数。

题解：利用矩阵乘法，可以预处理出从 X 到 Y 恰好经过 K 条边的最短路是多少。对于每次询问，直接处理处理即可，时间复杂度为 \(O(n^4)\)。

注意：恰好经过 K 条边的最短路不能将 G[i][i] 初始化成 0，因为边数有实际意义，若这样初始化意味着有自环出现。至少经过 K 条边的同理，也不能这样初始化。

代码如下

#include <bits/stdc++.h>
using namespace std;
const int maxn=55;
const int inf=0x3f3f3f3f;

int n,m,q;
struct mat{
    int d[maxn][maxn];
    mat(){memset(d,0x3f,sizeof(d));}
    int *operator[](int i){return d[i];}
    friend mat operator*(mat &x,mat &y){
        mat z;
        for(int i=1;i<=n;i++)
            for(int j=1;j<=n;j++)
                for(int k=1;k<=n;k++)
                    z[i][j]=min(z[i][j],x[i][k]+y[k][j]);
        return z;
    }
}d[maxn];

void read_and_parse(){
    scanf("%d%d",&n,&m);
    for(int i=1;i<=m;i++){
        int x,y,z;scanf("%d%d%d",&x,&y,&z);
        d[1][x][y]=min(d[1][x][y],z);
    }
}
void solve(){
    for(int i=2;i<=n;i++)d[i]=d[i-1]*d[1];
    scanf("%d",&q);
    while(q--){
        int x,y;scanf("%d%d",&x,&y);
        bool flag=0;double ans=1e18;
        for(int i=n;i;i--){
            if(d[i][x][y]!=inf){
                flag=1;
                ans=min(ans,(double)d[i][x][y]/(double)i);
            }
        }
        if(flag)printf("%.3lf\n",ans);
        else puts("OMG!");
    }
}
int main(){
    read_and_parse();
    solve();
    return 0;
}