HDU-2448 Mining Station on the Sea

先根据不同的起点跑最短路，记录距离，从而建立二分图求最小匹配。

一开始我求最短路的时候我把港口直接加到图中，然后发现进了港口就不能出来了，所以连接港口的边就要从双向边改成单向边…………这也搞得我n和m分不清了……

还不如排除掉港口算最短路后再统计各艘船到各个港口的最短距离……

然后我还傻叉地用了Dijkstra来求最短路，当然TLE咯……

事实证明，Floyd就可以了嘛我智商跑哪里去了。。。

#include <cstdlib>
#include <cstdio>
#include <cmath>
#include <cstring>
#include <algorithm>
#include <fstream>
#include <iostream>
#include <vector>
#include <cctype>
#include <queue>
#define rep(i, l, r) for(int i=l; i<=r; i++)
#define clr(x, c) memset(x, c, sizeof(x))
#define N 234
#define MAX 0x3fffffff
#define ll long long
using namespace std;
int read()
{
	int x=0, f=1; char ch=getchar();
	while (!isdigit(ch)) { if (ch=='-') f=-1; ch=getchar(); }
	while (isdigit(ch)) { x=x*10+ch-'0'; ch=getchar(); }
	return x*f;
}

int n, m, en, p, l[N], st[N], lx[N], ly[N], v[N][N], k[N], d[N][N], w[N];
bool vx[N], vy[N];

bool Find(int x)
{
	vx[x]=1;
	rep(y, 1, n) 
	{
		if (vy[y]) continue;
		int a=lx[x]+ly[y]-v[x][y];
		if (!a)
		{
			vy[y]=1; if (!l[y] || Find(l[y])) { l[y]=x; return 1; }
		}
		else st[y]=min(st[y], a);
	}
	return false;
}

inline int km()
{
	clr(l, 0); rep(i, 1, n) lx[i]=-MAX;
	rep(i, 1, n) rep(j, 1, n) if (lx[i]<v[i][j]) lx[i]=v[i][j]; rep(i, 1, n) if (lx[i]==-MAX) return 1;
	clr(ly, 1); rep(i, 1, n) rep(j, 1, n) if (v[i][j]!=-MAX) ly[j]=0; rep(i, 1, n) if (ly[i]) return 1;
	rep(i, 1, n) 
	{
		rep(j, 1, n) st[j]=MAX;
		while (1)
		{
			clr(vx, 0); clr(vy, 0);
			if (Find(i)) break; int a=MAX;
			rep(j, 1, n) if (!vy[j] && st[j]<a) a=st[j]; 
			if (a==MAX) return 1;
			rep(j, 1, n) if (vx[j]) lx[j]-=a;
			rep(j, 1, n) if (vy[j]) ly[j]+=a; else st[j]-=a;
		}
	}
	int a=0; rep(i, 1, n) a+=lx[i]+ly[i];
	return a;
}

int main()
{
	while (~scanf("%d%d%d%d", &n, &m ,&en, &p))
	{
		rep(i, 1, m) rep(j, 1, m) d[i][j] = i==j?0:MAX;
		rep(i, 1, n) w[i]=read();
		rep(i, 1, en) { int x=read(), y=read(), z=read(); if (d[x][y]>z) d[x][y]=d[y][x]=z; }
		rep(o, 1, m) rep(i, 1, m) rep(j, 1, m) if (d[i][o]+d[o][j]<d[i][j]) d[i][j]=d[i][o]+d[o][j];
		
		rep(i, 1, n) rep(j, 1, n) v[i][j]=MAX;
		rep(j, 1, p) 
		{
			int x=read(), y=read(), z=read();
			rep(i, 1, n) if (v[i][x]>d[w[i]][y]+z) 
				v[i][x]=d[w[i]][y]+z;
		}
		rep(i, 1, n) rep(j, 1, n) v[i][j]*=-1;
		printf("%d
", -km());
	}
	return 0;
}