传送门：http://oj.cnuschool.org.cn/oj/home/problem.htm?problemID=302
| The Flash |
| 难度级别：B； 运行时间限制：1000ms； 运行空间限制：51200KB； 代码长度限制：2000000B |
|
试题描述
|
|
"My name is Barry Allen and I'm the fastest man alive. When I was a child, I saw my mom killed by something impossible, and my father went to prison for her murder. Then an accident made me the impossible. To the outside world, I'm an ordinary forensic scientist. But secretly, I used my speed to fight crim and find others like me. And one day, I'll find who killed my mother and get justice for my father. I am the flash." 以上是xjr背诵的闪电侠片头的片段。
B.I.T.C.H.们玩得正High时，闪电侠来得瑟他的速度了。他邀请xjr和Avengers观看他夜晚在城市跑步的场景。 他的路线是这样的：（一共有n个城市，城市由1到n编号，城市之间有m条边） 1 → 2 → 1 → 3 → 1 → 4 → 1 → 5 → 1 → … → n-1 → 1 → n → 1 但由于闪电侠速度超过闪电，他也有了一种和电流一样的心理——走最短路，现在你需要回答闪电侠最少需要走多少千米的路程。 |
|
输入
|
|
第一行，两个整数n和m（见试题描述）
接下来m行，每行三个整数a，b和c，表示从a城市到b城市有一条有向边，长度为c |
|
输出
|
|
闪电侠走的最短路程
|
|
输入示例
|
|
5 10
2 3 8 1 5 90 3 5 82 1 2 76 1 3 8 5 3 4 4 1 23 4 5 6 3 5 6 5 4 2 |
|
输出示例
|
|
232
|
|
其他说明
|
|
1 < n, c < 1001
1 < m < 100001 数据保证城市1能到达任何城市，并且任何城市能到达城市1 |
题解：最短路瞎跑大水题（出题人竟然没有卡SPFA！！！）
更新邻接表：
1 #include<iostream> 2 #include<cstdio> 3 #include<cmath> 4 #include<algorithm> 5 #include<queue> 6 #include<cstring> 7 #define PAU putchar(' ') 8 #define ENT putchar(' ') 9 using namespace std; 10 const int maxn=1000+10,maxm=100000+10,inf=1e9; 11 struct SPFA{ 12 struct ted{int x,y,w;ted*nxt;}adj[maxm],*fch[maxn],*ms; 13 bool inq[maxn];int d[maxn]; 14 void init(int n){ 15 ms=adj;fill(d,d+n+1,inf);memset(inq,false,sizeof(inq));return; 16 } 17 void ade(int u,int v,int w){ 18 *ms=(ted){u,v,w,fch[u]};fch[u]=ms++; 19 return; 20 } 21 void spfa(int S){ 22 queue<int>Q;Q.push(S);d[S]=0; 23 while(!Q.empty()){ 24 int u=Q.front();Q.pop();inq[u]=false; 25 for(ted*e=fch[u];e;e=e->nxt){ 26 int v=e->y; 27 if(d[v]>d[u]+e->w){ 28 d[v]=d[u]+e->w; 29 if(!inq[v]) Q.push(v),inq[v]=true; 30 } 31 } 32 } return; 33 } 34 }p1,p2; 35 inline int read(){ 36 int x=0,sig=1;char ch=getchar(); 37 while(!isdigit(ch)){if(ch=='-')sig=-1;ch=getchar();} 38 while(isdigit(ch))x=10*x+ch-'0',ch=getchar(); 39 return x*=sig; 40 } 41 inline void write(int x){ 42 if(x==0){putchar('0');return;}if(x<0)putchar('-'),x=-x; 43 int len=0,buf[20];while(x)buf[len++]=x%10,x/=10; 44 for(int i=len-1;i>=0;i--)putchar(buf[i]+'0');return; 45 } 46 int n,m; 47 void init(){ 48 n=read();m=read();p1.init(n);p2.init(n); 49 int x,y,w; 50 for(int i=1;i<=m;i++){ 51 x=read();y=read();w=read(); 52 p1.ade(x,y,w);p2.ade(y,x,w); 53 } 54 p1.spfa(1);p2.spfa(1); 55 int ans=0; 56 for(int i=1;i<=n;i++){ 57 ans+=p1.d[i]+p2.d[i]; 58 } write(ans); 59 return; 60 } 61 void work(){ 62 return; 63 } 64 void print(){ 65 return; 66 } 67 int main(){init();work();print();return 0;}
数组版：
1 #include<iostream> 2 #include<cstdio> 3 #include<cmath> 4 #include<algorithm> 5 #include<queue> 6 #include<cstring> 7 #define inf 10000000 8 #define PAU putchar(' ') 9 #define ENT putchar(' ') 10 using namespace std; 11 const int maxn=1000+10,maxm=100000+10; 12 struct ShortiestPath{ 13 struct Tedge{int x,y,w,next;}adj[maxm];int fch[maxn],ms; 14 int d[maxn],S,n;bool inque[maxn]; 15 void AddEdge(int u,int v,int w){adj[++ms]=(Tedge){u,v,w,fch[u]};fch[u]=ms;return;} 16 void init(int S,int n){ 17 this->S=S;this->n=n;ms=0; 18 memset(inque,false,sizeof(inque)); 19 memset(fch,0,sizeof(fch)); 20 for(int i=1;i<=n;i++) d[i]=inf; 21 return; 22 } 23 void SPFA(){ 24 queue<int>Q;d[S]=0;Q.push(S); 25 while(!Q.empty()){ 26 int u=Q.front();Q.pop();inque[u]=false; 27 for(int i=fch[u];i;i=adj[i].next){ 28 int v=adj[i].y; 29 if(d[v]>d[u]+adj[i].w){ 30 d[v]=d[u]+adj[i].w; 31 if(!inque[v]){ 32 inque[v]=true; 33 Q.push(v); 34 } 35 } 36 } 37 } return; 38 } 39 }p1,p2; 40 inline int read(){ 41 int x=0,sig=1;char ch=getchar(); 42 while(!isdigit(ch)){if(ch=='-')sig=-1;ch=getchar();} 43 while(isdigit(ch))x=10*x+ch-'0',ch=getchar(); 44 return x*=sig; 45 } 46 inline void write(int x){ 47 if(x==0){putchar('0');return;}if(x<0)putchar('-'),x=-x; 48 int len=0,buf[15];while(x)buf[len++]=x%10,x/=10; 49 for(int i=len-1;i>=0;i--)putchar(buf[i]+'0');return; 50 } 51 int n,m; 52 void init(){ 53 n=read();m=read(); 54 p1.init(1,n);p2.init(1,n); 55 for(int i=1;i<=m;i++){ 56 int a=read(),b=read(),c=read(); 57 p1.AddEdge(a,b,c); 58 p2.AddEdge(b,a,c); 59 } 60 return; 61 } 62 void work(){ 63 p1.SPFA();p2.SPFA(); 64 return; 65 } 66 void print(){ 67 int ans=0; 68 for(int i=1;i<=n;i++){ 69 ans+=p1.d[i]+p2.d[i]; 70 } 71 write(ans); 72 return; 73 } 74 int main(){init();work();print();return 0;}
当然啦还是写Dijkstra比较稳妥，小健建のDijkstra：（原谅我懒到一定境界了。。。）
1 #include<cstdio> 2 #include<cstring> 3 #include<queue> 4 #include<algorithm> 5 using namespace std; 6 const int maxn=1010; 7 const int maxm=100010; 8 const int INF=100000000; 9 struct Edge 10 { 11 int from,to,dist; 12 }; 13 struct Heapnode 14 { 15 int d,u; 16 bool operator < (const Heapnode& rhs) const 17 { 18 return d>rhs.d; 19 } 20 }; 21 struct Dijkstra 22 { 23 int n,m; 24 int p[maxn]; 25 Edge edges[maxm]; 26 int first[maxn],next[maxm]; 27 int d[maxn]; 28 bool done[maxn]; 29 void init(int n) 30 { 31 this->n=n; 32 m=0; 33 memset(first,0,sizeof(first)); 34 } 35 void AddEdge(int from,int to,int dist) 36 { 37 edges[++m]=(Edge){from,to,dist}; 38 next[m]=first[from]; 39 first[from]=m; 40 } 41 void dijkstra(int s) 42 { 43 memset(done,0,sizeof(done)); 44 priority_queue<Heapnode> Q; 45 for(int i=1;i<=n;i++) d[i]=INF; 46 d[s]=0; 47 Q.push((Heapnode){0,s}); 48 while(!Q.empty()) 49 { 50 Heapnode x=Q.top(); 51 Q.pop(); 52 if(done[x.u]) continue; 53 done[x.u]=1; 54 for(int i=first[x.u];i;i=next[i]) 55 { 56 Edge& v=edges[i]; 57 if(d[v.to]>d[x.u]+v.dist) 58 { 59 d[v.to]=d[x.u]+v.dist; 60 Q.push((Heapnode){d[v.to],v.to}); 61 } 62 } 63 } 64 } 65 }sol1,sol2; 66 int main() 67 { 68 int n,m,a,b,c; 69 scanf("%d%d",&n,&m); 70 sol1.init(n); sol2.init(n); 71 while(m--) 72 { 73 scanf("%d%d%d",&a,&b,&c); 74 sol1.AddEdge(a,b,c); 75 sol2.AddEdge(b,a,c); 76 } 77 sol1.dijkstra(1); sol2.dijkstra(1); 78 long long ans=0; 79 for(int i=1;i<=n;i++) ans+=sol1.d[i]+sol2.d[i]; 80 printf("%lld ",ans); 81 return 0; 82 }