hdoj 1874 畅通工程续

题目链接:http://acm.hdu.edu.cn/showproblem.php?pid=1874

解题思路:SPFA 算法(邻接矩阵)

  1 /**************************************************************************
  2 user_id: SCNU20102200088
  3 problem_id: hdoj 1874
  4 problem_name: 畅通工程续
  5 **************************************************************************/
  6 
  7 #include <algorithm>
  8 #include <iostream>
  9 #include <iterator>
 10 #include <iomanip>
 11 #include <sstream>
 12 #include <fstream>
 13 #include <cstring>
 14 #include <cstdlib>
 15 #include <climits>
 16 #include <bitset>
 17 #include <string>
 18 #include <vector>
 19 #include <cstdio>
 20 #include <cctype>
 21 #include <ctime>
 22 #include <cmath>
 23 #include <queue>
 24 #include <stack>
 25 #include <list>
 26 #include <set>
 27 #include <map>
 28 using namespace std;
 29 
 30 //线段树
 31 #define lson l,m,rt<<1
 32 #define rson m+1,r,rt<<1|1
 33 
 34 //手工扩展栈
 35 #pragma comment(linker,"/STACK:102400000,102400000")
 36 
 37 const double EPS=1e-9;
 38 const double PI=acos(-1.0);
 39 const double E=2.7182818284590452353602874713526;  //自然对数底数
 40 const double R=0.5772156649015328606065120900824;  //欧拉常数:(1+1/2+...+1/n)-ln(n)
 41 
 42 const int x4[]={-1,0,1,0};
 43 const int y4[]={0,1,0,-1};
 44 const int x8[]={-1,-1,0,1,1,1,0,-1};
 45 const int y8[]={0,1,1,1,0,-1,-1,-1};
 46 
 47 typedef long long LL;
 48 
 49 typedef int T;
 50 T max(T a,T b){ return a>b? a:b; }
 51 T min(T a,T b){ return a<b? a:b; }
 52 T gcd(T a,T b){ return b==0? a:gcd(b,a%b); }
 53 T lcm(T a,T b){ return a/gcd(a,b)*b; }
 54 
 55 ///////////////////////////////////////////////////////////////////////////
 56 //Add Code:
 57 const int maxn=205;
 58 const int INF=2000000;
 59 int n,dist[maxn],f[maxn][maxn];
 60 bool visit[maxn];
 61 
 62 void SPFA(int k){
 63     int i;
 64     for(i=0;i<n;i++){
 65         dist[i]=INF;
 66         visit[i]=0;
 67     }
 68     queue<int> q;
 69     dist[k]=0;
 70     q.push(k);
 71     visit[k]=1;
 72     while(!q.empty()){
 73         int now=q.front();
 74         q.pop();
 75         visit[now]=0;
 76         for(i=0;i<n;i++){
 77             if(dist[now]+f[now][i]<dist[i]){
 78                 dist[i]=dist[now]+f[now][i];
 79                 if(!visit[i]){
 80                     q.push(i);
 81                     visit[i]=1;
 82                 }
 83             }
 84         }
 85     }
 86 }
 87 ///////////////////////////////////////////////////////////////////////////
 88 
 89 int main(){
 90     std::ios::sync_with_stdio(false);
 91     //freopen("in.txt","r",stdin);
 92     //freopen("out.txt","w",stdout);
 93     ///////////////////////////////////////////////////////////////////////
 94     //Add Code:
 95     int m,a,b,c,start,end,i,j;
 96     while(scanf("%d%d",&n,&m)!=EOF){
 97         for(i=0;i<n;i++) f[i][i]=0;
 98         for(i=0;i<n;i++){
 99             for(j=i+1;j<n;j++) f[i][j]=f[j][i]=INF;
100         }
101         while(m--){
102             scanf("%d%d%d",&a,&b,&c);
103             if(c<f[a][b]) f[a][b]=f[b][a]=c;
104         }
105         scanf("%d%d",&start,&end);
106         SPFA(start);
107         if(dist[end]==INF) printf("-1
");
108         else printf("%d
",dist[end]);
109     }
110     ///////////////////////////////////////////////////////////////////////
111     return 0;
112 }
113 
114 /**************************************************************************
115 Testcase:
116 Input:
117 3 3
118 0 1 1
119 0 2 3
120 1 2 1
121 0 2
122 3 1
123 0 1 1
124 1 2
125 Output:
126 2
127 -1
128 **************************************************************************/
原文地址:https://www.cnblogs.com/linqiuwei/p/3339298.html