HDU 2448 Mining Station on the Sea

题意：有N艘船和N个港口，M个station。现在N艘船分别处在一些station上。有K条边连接站，有P条边连接港口和站。现在这N艘船都要回到港口，且每个港口只能容纳一艘船。问使这N艘船回到港口行程之和的最小值是多少。

比较明显的求最小权匹配，用KM。做最短路的时候要注意如果船已经回到港口了，就不能再跑出去了，所以做floyd闭包或其他最短路的时候，松弛的时候不能用港口去松弛一个点对。然后得到各个船所在位置到各个港口的最短路后，把边权赋成负的船到港口的最短路值，然后跑KM就可以了。

#include <cstdio>
#include <vector>
#include <cstring>
#include <algorithm>
#define maxn 105
#define INF 1<<30
using namespace std;

struct KM{
    vector<int> G[maxn];
    int W[maxn][maxn],n;
    int Lx[maxn],Ly[maxn];
    int left[maxn];
    bool S[maxn],T[maxn];

    void init(int n)
    {
        this->n = n;
        for(int i = 0;i < n;i++)    G[i].clear();
        memset(W,0,sizeof(W));
    }

    void add_edge(int u,int v,int w)
    {
        G[u].push_back(v);
        W[u][v] = w;
    }

    bool match(int i)
    {
        S[i] = true;
        for(int k = 0;k < G[i].size();k++)
        {
            int j = G[i][k];
            if(Lx[i] + Ly[j] == W[i][j] && !T[j])
            {
                T[j] = true;
                if(left[j] == -1 || match(left[j]))
                {
                    left[j] = i;
                    return true;
                }
            }
        }
        return false;
    }
    void update()
    {
        int a = INF;
        for(int i = 0;i < n;i++)    if(S[i])
            for(int k = 0;k < G[i].size();k++)
            {
                int j = G[i][k];    if(!T[j])
                a = min(a,Lx[i] + Ly[j] - W[i][j]);
            }
        for(int i = 0;i < n;i++)
        {
            if(S[i])    Lx[i] -= a;
            if(T[i])    Ly[i] += a;
        }
    }

    void solve()
    {
        for(int i = 0;i < n;i++)
        {
            Lx[i] = *max_element(W[i],W[i]+n);
            left[i] = -1;
            Ly[i] = 0;
        }
        for(int i = 0;i < n;i++)
        {
            while(1)
            {
                for(int j = 0;j < n;j++)    S[j] = T[j] = 0;
                if(match(i))    break;
                else            update();
            }
        }
    }
};
KM solver;
int d[310][310];
int pos[105];

int main(){
    int n,m,k,p;
    while(scanf("%d%d%d%d",&n,&m,&k,&p) != EOF){
        for(int i = 1;i <= n;i++)   scanf("%d",&pos[i]);

        for(int i = 1;i <= m+n;i++)
        for(int j = 1;j <= m+n;j++){
            if(i != j)  d[i][j] = -1;
            else    d[i][j] = 0;
        }

        for(int i = 0;i < k;i++){
            int a,b,c;
            scanf("%d%d%d",&a,&b,&c);
            d[a][b] = d[b][a] = d[a][b] == -1 ? c : min(c,d[a][b]);
        }

        for(int i = 0;i < p;i++){
            int a,b,c;
            scanf("%d%d%d",&a,&b,&c);
            d[m+a][b] = d[b][m+a] = d[m+a][b] == -1 ? c : min(c,d[m+a][b]);
        }

        for(int k = 1;k <= m;k++)
        for(int i = 1;i <= m+n;i++) if(d[i][k] != -1)
        for(int j = 1;j <= m+n;j++) if(d[k][j] != -1 && (d[i][j] == -1 || d[i][j] > d[i][k] + d[k][j]))
        d[i][j] = d[i][k] + d[k][j];

        solver.init(n);

        for(int i = 1;i <= n;i++)
        for(int j = 1;j <= n;j++)
        if(d[pos[i]][m+j] != -1)
        solver.add_edge(i-1,j-1,-d[pos[i]][m+j]);

        solver.solve();

        int ans = 0;
        for(int i = 0;i < n;i++)
        ans -= solver.W[solver.left[i]][i];
        printf("%d
",ans);
    }
    return 0;
}