NYOJ-115 城市平乱 SPFA

城市平乱

时间限制：1000 ms  |  内存限制：65535 KB

难度：4

描述

南将军统领着N个部队，这N个部队分别驻扎在N个不同的城市。

他在用这N个部队维护着M个城市的治安，这M个城市分别编号从1到M。

现在，小工军师告诉南将军，第K号城市发生了暴乱，南将军从各个部队都派遣了一个分队沿最近路去往暴乱城市平乱。

现在已知在任意两个城市之间的路行军所需的时间，你作为南将军麾下最厉害的程序员，请你编写一个程序来告诉南将军第一个分队到达叛乱城市所需的时间。

注意，两个城市之间可能不只一条路。

输入

第一行输入一个整数T，表示测试数据的组数。(T<20)
每组测试数据的第一行是四个整数N,M,P,Q(1<=N<=100,N<=M<=1000,M-1<=P<=100000)其中N表示部队数，M表示城市数，P表示城市之间的路的条数,Q表示发生暴乱的城市编号。
随后的一行是N个整数，表示部队所在城市的编号。
再之后的P行，每行有三个正整数，a,b,t(1<=a,b<=M,1<=t<=100)，表示a,b之间的路如果行军需要用时为t

数据保证暴乱的城市是可达的。

输出

对于每组测试数据，输出第一支部队到达叛乱城市时的时间。每组输出占一行

样例输入

13 8 9 81 2 31 2 12 3 21 4 22 5 33 6 24 7 15 7 35 8 26 8 2 

样例输出

4

邻接表，实现SPFA算法。（以前只知道Dijkstra算法，菜鸟就是孤陋寡闻啊）

01.#include <iostream>

02.#include <vector>

03.#include <utility>

04.#include <fstream>

05.#include <queue>

06.using namespace std;

07.#define inf 100000000

08.#define maxM 1002

09.//#define DEBUG

10.int N,M,P,Q,T,city[maxM],dist[maxM],InQueue[maxM];

11.vector< pair<int,int> > c[maxM];

12. 

13.void FileIn(istream &in)

14.{

15.in>>N>>M>>P>>Q;

16.for(int i=0;i<N;i++)in>>city[i];

17.for(int i=0;i<maxM;i++)

18.{

19.c[i].clear();

20.}

21.int a,b,t;

22.for(int i=0;i<P;i++)

23.{

24.in>>a>>b>>t;

25.c[a].push_back(make_pair(b,t));

26.c[b].push_back(make_pair(a,t));

27.}

28.}

29. 

30.void SPFA()

31.{

32.for(int i=0;i<maxM;i++)

33.{

34.dist[i]=inf;

35.InQueue[i]=0;

36.}

37.dist[Q]=0;

38.queue<int> q;

39.q.push(Q);

40.InQueue[Q]=1;

41.int u,v,u2v;

42.while(!q.empty())

43.{

44.u=q.front();

45.q.pop();

46.InQueue[u]=0;

47.for(int i=0;i<c[u].size();i++)

48.{

49.v=c[u][i].first;

50.u2v=c[u][i].second;

51.if(dist[u]+u2v<dist[v])

52.{

53.dist[v]=dist[u]+u2v;

54.if(!InQueue[v])

55.{

56.q.push(v);

57.InQueue[v]=1;

58.}

59.}

60.}

61.}

62.}

63. 

64.int main()

65.{

66.#ifdef DEBUG

67.ifstream in("data.txt");

68.#else

69.istream &in=cin;

70.#endif

71.in>>T;

72.while(T--)

73.{

74.FileIn(in);

75.SPFA();

76.int min=inf;

77.for(int i=0;i<N;i++)

78.{

79.if(dist[city[i]]<min)

80.min=dist[city[i]];

81.}

82.cout<<min<<endl;

83.}

84.#ifdef DEBUG

85.system("PAUSE");

86.#endif

87.return 0;

88.}