hdu - 4370 神建图最短路 spfa改进 邻接矩阵存图

题意：有一个n*n的01矩阵，满足以下条件

            1.X₁₂+X₁₃+...X_1n=1
            2.X_1n+X_2n+...X_n-1n=1
            3.for each i (1<i<n), satisfies ∑X_ki (1<=k<=n)=∑X_ij (1<=j<=n).

            另给出一个矩阵C，求∑C_ij*X_ij(1<=i,j<=n)的最小值。

神一样的题，刚开始看就没有想到会是最短路里的题。。。。。。（虽然是在最短路专题里。。。）

难就难在思维转换吧。

将矩阵C看做描述N个点花费的邻接矩阵
再来看三个条件：
条件一：表示1号点出度为1
条件二：表示n号点入度为1
条件三：表示k( 1 < k < n )号点出度等于入度
最后再来看看题目要求，∑Cij*Xij(1<=i,j<=n)，很明显，这是某个路径的花费，而路径的含义可以有以下两种：
一：1号点到n号点的花费

二：1号点经过其它点成环，n号点经过其它点成环，这两个环的花费之和

于是，就变成了一道简单的最短路问题
关于环花费的算法，可以改进spfa算法，初始化dis[start] = INF，且一开始让源点之外的点入队

链接：hdu 4370

#include<iostream>
#include<cstdio>
#include<algorithm>
#include<queue>
#include<cstring>
using namespace std;

const int maxn = 300 + 10;
const int inf = 0x3f3f3f3f;
int n;
int cost[maxn][maxn];
int d[maxn], inq[maxn];

void spfa(int s)      //改进的spfa
{
    memset(inq, 0, sizeof(inq));
    queue<int> q;
    for(int i = 1; i <= n; i++) {
        if(i == s) {
            d[s] = inf;
            continue;
        }
        d[i] = cost[s][i];
        inq[i] = 1;
        q.push(i);
    }
    while(!q.empty()) {
        int u = q.front();
        q.pop();
        inq[u] = 0;
        for(int i = 1; i <= n; i++) {
            if(d[i] > d[u] + cost[u][i]) {
                d[i] = d[u] + cost[u][i];
                if(inq[i]) continue;
                inq[i] = 1;
                q.push(i);
            }
        }
    }
}

int main()
{
    while(~scanf("%d", &n)) {
        for(int i = 1; i <= n; i++)
         for(int j = 1; j <= n; j++)
         scanf("%d", &cost[i][j]);
        int ans = inf;
        spfa(1);
        int loop1 = d[1];
        ans = d[n];
        spfa(n);
        int loopn = d[n];
        ans = min(ans, loop1 + loopn);
        printf("%d\n", ans);
    }
    return 0;
}