USACO Sweet Butter, 难题，Bellman-Ford算法 + SPFA优化

这题碰到我头晕的时候，搞了几天没搞懂，参考了网上的代码
nocow的解释：http://www.nocow.cn/index.php/USACO/butter

Bellman-Ford算法：http://zh.wikipedia.org/wiki/Bellman-Ford%E7%AE%97%E6%B3%95
SPFA的讲解：http://www.cnblogs.com/zen_chou/archive/2009/05/15/1457962.html

实在是没弄懂……

代码

/* 
ID: wangxin12 
PROG: butter
LANG: C++ 
*/ 
#include <iostream>
#include <fstream>
#include <cstring>
#include <queue>
using namespace std;

const int MAX = 805;
const int INF = 0x0fffffff;
ifstream fin("butter.in");
ofstream fout("butter.out");

struct node {
	int end, length;
};

int counter[MAX] = { 0 }, location[505] = { 0 }, n, p, c;
node adj[MAX][MAX];
bool in[MAX];
int dist[MAX];

int SPFA(int i) {
	memset(in, 0, sizeof(in));
	for(int k = 1; k <= p; k++)  dist[k] = INF;
	queue<int> Q;
	Q.push(i);

	in[i] = true;
	dist[i] = 0;
	
	while( !Q.empty() ) {
		int x = Q.front(); Q.pop();
		in[x] = false;
		for(int j = 0; j < counter[x]; j++) {
			if(dist[x] + adj[x][j].length < dist[ adj[x][j].end ] ) {
				dist[adj[x][j].end] = dist[x] + adj[x][j].length;
				if( !in[adj[x][j].end] ) {
					Q.push(adj[x][j].end);
					in[adj[x][j].end] = true;
				}
			}
		}
	}

	int ans = 0;
	for(int k = 1; k <= n; k++) {
		if(dist[location[k]] == INF) return -1;
		else ans += dist[location[k]];
	}

	return ans;
}





int main() {
	memset(counter, 0, sizeof(counter));
	fin>>n>>p>>c;
	for(int i = 1; i <= n; i++)
		fin>>location[i];

	for(int i = 1, s, t, value; i <= c; i++) {
		fin>>s>>t>>value;
		adj[s][counter[s] ].end = t; adj[s][counter[s]].length = value; counter[s]++;
		adj[t][counter[t] ].end = s; adj[t][counter[t]].length = value; counter[t]++;
	}

	int tt, min = INF;

	for(int i = 1; i <= p; i++) {
		tt = SPFA(i);
		if(tt < min && tt != -1) min = tt;
	}
	fout<<min<<endl;
	
	fin.close();
	fout.close();
	return 0;
}

题目翻译

描述

农夫John发现做出全威斯康辛州最甜的黄油的方法：糖。把糖放在一片牧场上，他知道N（1<=N<=500）只奶牛会过来舔它，这样就能做出能卖好价钱的超甜黄油。当然，他将付出额外的费用在奶牛上。

农夫John很狡猾。他知道他可以训练这些奶牛，让它们在听到铃声时去一个特定的牧场。他打算将糖放在那里然后下午发出铃声，以至他可以在晚上挤奶。

农夫John知道每只奶牛都在各自喜欢的牧场呆着（一个牧场不一定只有一头牛）。给出各头牛在的牧场和牧场间的路线，找出使所有牛到达的路程和最短的牧场（他将把糖放在那）。

[编辑]格式

PROGRAM NAME: butter

INPUT FORMAT:

(file butter.in)

第一行: 三个数：奶牛数N，牧场数P（2<=P<=800），牧场间道路数C(1<=C<=1450).

第二行到第N+1行: 1到N头奶牛所在的牧场号.

第N+2行到第N+C+1行： 每行有三个数：相连的牧场A、B，两牧场间距（1<=D<=255），当然,连接是双向的.

OUTPUT FORMAT:

(file butter.out)

一行 输出奶牛必须行走的最小的距离和.

[编辑]SAMPLE INPUT

3 4 5
2
3
4
1 2 1
1 3 5
2 3 7
2 4 3
3 4 5

样例图形

         P2  
P1 @--1--@ C1
    \    |\
     \   | \
      5  7  3
       \ |   \
        \|    \ C3
      C2 @--5--@
         P3    P4

[编辑]SAMPLE OUTPUT

8

{说明： 放在4号牧场最优. }