[USACO Mar08] 牛跑步
Bessie准备用从牛棚跑到池塘的方法来锻炼. 但是因为她懒,她只准备沿着下坡的路跑到池塘,然后走回牛棚.
Bessie也不想跑得太远,所以她想走最短的路经. 农场上一共有M(1<=M<=10,000)条路,每条路连接两个用1..N(1<=N<=1000)标号的地点. 更方便的是,如果X>Y,则地点X的高度大于地点Y的高度. 地点N是Bessie的牛棚;地点1是池塘.
很快, Bessie厌倦了一直走同一条路.所以她想走不同的路,更明确地讲,她想找出K(1<=K<=100)条不同的路经.为了避免过度劳累,她想使这K条路径为最短的K条路径.
请帮助Bessie找出这K条最短路经的长度.你的程序需要读入农场的地图, 一些从Xi到Yi的路径和它们的长度(Xi,Yi,Di). 所有(Xi,Yi,Di)满足(1<=Yi<Xi;Yi<Xi<=N,1<=Di<=1,000,000).
题目名称: cowjog
输入格式:
- 第1行: 3个数: N,M,K
- 第2..M+1行: 第 i+1行包含3个数 Xi,Yi,Di, 表示一条下坡的路.
样例输入 (cowjog.in):
5 8 7 5 4 1 5 3 1 5 2 1 5 1 1 4 3 4 3 1 1 3 2 1 2 1 1
输出格式:
- 第1..K行: 第i行包含第i最短路径的长度,或−1如果这样的路径不存在.如果多条路径有同样的长度,请注意将这些长度逐一列出.
样例输出 (cowjog.out):
1 2 2 3 6 7 -1
输出解释:
路径分别为(5−1),(5−3−1),(5−2−1),(5−3−2−1),(5−4−3−1),(5−4−3−2−1)
1 //0.002s 2 3 #include <queue> 4 #include <cstdio> 5 #include <ctype.h> 6 #include <cstring> 7 #include <algorithm> 8 9 using namespace std; 10 11 const int MAXM=10010; 12 const int MAXN=1010; 13 14 int n,m,k; 15 16 int dis[MAXN]; 17 18 bool vis[MAXN]; 19 20 struct SKT { 21 int v; 22 int dist; 23 bool operator < (const SKT&p) const { 24 return dist+dis[v]>p.dist+dis[p.v]; 25 } 26 }; 27 SKT s; 28 29 struct edge { 30 int to; 31 int next; 32 int val; 33 edge() {} 34 edge(int to,int val,int next):to(to),val(val),next(next){} 35 }; 36 edge e[MAXM<<1],r[MAXM<<1]; 37 38 int head[MAXN],Head[MAXN],tot; 39 40 inline void add(int x,int y,int v) { 41 e[++tot]=edge(y,v,head[x]); 42 r[tot]=edge(x,v,Head[y]); 43 head[x]=Head[y]=tot; 44 } 45 46 inline void read(int&x) { 47 int f=1;register char c=getchar(); 48 for(x=0;!isdigit(c);c=='-'&&(f=-1),c=getchar()); 49 for(;isdigit(c);x=x*10+c-48,c=getchar()); 50 x=x*f; 51 } 52 53 void SPFA() { 54 queue<int> Q; 55 for(int i=1;i<=n;++i) dis[i]=0x7fffffff; 56 dis[1]=0; 57 vis[1]=true; 58 Q.push(1); 59 while(!Q.empty()) { 60 int now=Q.front(); 61 Q.pop(); 62 vis[now]=false; 63 for(int i=Head[now];i!=-1;i=r[i].next) { 64 int to=r[i].to; 65 if(dis[to]>dis[now]+r[i].val) { 66 dis[to]=dis[now]+r[i].val; 67 if(!vis[to]) Q.push(to),vis[to]=true; 68 } 69 } 70 } 71 return; 72 } 73 74 void Astar() { 75 int ans=0; 76 priority_queue<SKT> Q; 77 s.v=n,s.dist=0; 78 Q.push(s); 79 while(!Q.empty()) { 80 SKT now=Q.top(); 81 Q.pop(); 82 if(now.v==1) { 83 printf("%d ",now.dist); 84 ++ans; 85 } 86 if(ans==k) return; 87 for(int i=head[now.v];i!=-1;i=e[i].next) { 88 int to=e[i].to; 89 if(to<now.v) { 90 SKT FAKER=now; 91 FAKER.v=to;FAKER.dist=now.dist+e[i].val; 92 Q.push(FAKER); 93 } 94 } 95 } 96 while(ans<k) printf("-1 "),++ans; 97 return; 98 } 99 100 int hh() { 101 freopen("cowjog.in","r",stdin); 102 freopen("cowjog.out","w",stdout); 103 memset(Head,-1,sizeof Head); 104 memset(head,-1,sizeof head); 105 read(n);read(m);read(k); 106 for(int x,y,z;m--;) { 107 read(x);read(y);read(z); 108 if(x>y) add(x,y,z); 109 } 110 SPFA(); 111 Astar(); 112 return 0; 113 } 114 115 int sb=hh(); 116 int main() {;}
1 #include <queue> 2 #include <cstdio> 3 #include <ctype.h> 4 #include <cstring> 5 #include <algorithm> 6 7 using namespace std; 8 9 const int MAXM=10010; 10 const int MAXN=1010; 11 12 int n,m,k; 13 14 int dis[MAXN]; 15 16 bool vis[MAXN]; 17 18 struct SKT { 19 int v,val,h; 20 int dist; 21 SKT() {} 22 SKT(int v,int dist):v(v),dist(dist) {} 23 SKT(int x,int y,int z) { 24 v=x;val=y;h=z; 25 dist=val+h; 26 } 27 bool operator < (const SKT&p) const { 28 return dist>p.dist; 29 } 30 }; 31 SKT s; 32 33 struct edge { 34 int to; 35 int next; 36 int val; 37 edge() {} 38 edge(int to,int val,int next):to(to),val(val),next(next){} 39 }; 40 edge e[MAXM<<1],r[MAXM<<1]; 41 42 int head[MAXN],Head[MAXN],tot; 43 44 inline void add(int x,int y,int v) { 45 e[++tot]=edge(y,v,head[x]); 46 r[tot]=edge(x,v,Head[y]); 47 head[x]=Head[y]=tot; 48 } 49 50 inline void read(int&x) { 51 int f=1;register char c=getchar(); 52 for(x=0;!isdigit(c);c=='-'&&(f=-1),c=getchar()); 53 for(;isdigit(c);x=x*10+c-48,c=getchar()); 54 x=x*f; 55 } 56 57 void Dijs() { 58 priority_queue<SKT> Q; 59 for(int i=1;i<=n;++i) dis[i]=0x3f3f3f3f; 60 Q.push(SKT(1,0)); 61 dis[1]=0; 62 while(!Q.empty()) { 63 SKT now=Q.top(); 64 Q.pop(); 65 for(int i=Head[now.v];i!=-1;i=r[i].next) { 66 int to=r[i].to; 67 if(dis[to]>dis[now.v]+r[i].val) { 68 dis[to]=dis[now.v]+r[i].val; 69 Q.push(SKT(to,dis[to])); 70 } 71 } 72 } 73 return; 74 } 75 76 void Astar() { 77 int ans=0; 78 priority_queue<SKT> Q; 79 s.v=n,s.dist=0; 80 Q.push(s); 81 while(!Q.empty()) { 82 SKT now=Q.top(); 83 Q.pop(); 84 if(now.v==1) {printf("%d ",now.dist);++ans;} 85 if(ans==k) return; 86 for(int i=head[now.v];i!=-1;i=e[i].next) 87 Q.push(SKT(e[i].to,now.val+e[i].val,dis[e[i].to])); 88 } 89 while(ans<k) printf("-1 "),++ans; 90 return; 91 } 92 93 int hh() { 94 freopen("cowjog.in","r",stdin); 95 freopen("cowjog.out","w",stdout); 96 memset(Head,-1,sizeof Head); 97 memset(head,-1,sizeof head); 98 read(n);read(m);read(k); 99 for(int x,y,z;m--;) { 100 read(x);read(y);read(z); 101 add(x,y,z); 102 } 103 Dijs(); 104 Astar(); 105 return 0; 106 } 107 108 int sb=hh(); 109 int main() {;}